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THE BOUNDARY LAYER IN THE CONVERGING 
NOZZLE OF A SWIRL ATOMIZER 
By G. 1. TAYLOR (Trinity College, Cambridge) 
[Received 15 November 1949] 


SUMMARY* 


Swirl atomizers are used in agricultural spraying machinery and in oil-fired 
furnaces. The swirl is produced by leading in the liquid tangentially at the outside 
of achamber and allowing it to spray out through a central orifice of small diameter. 
The liquid is foreed towards the centre by high pressure in the supply pipe. When 
aliquid which is not swirling approaches an orifice through a converging chamber, 
it converges from all points of a cross-section, only decreasing in velocity in the 
boundary layer near the wall. When the liquid is swirling, however, the radial 
pressure gradient, which necessarily accompanies the swirling motion and holds the 
liquid particles in their circular paths, acts on the retarded boundary layer, driving 
it along the surface of the chamber towards the orifice. Thus a condition can arise 
in which practically the whole of the outflow from the orifice is fed by a boundary- 
layer current close to the surface of the swirl chamber. The study of boundary layers 
containing transverse as well as longitudinal components of velocity is of interest 
apart from its practical importance, and one case has been investigated completely, 
namely, the boundary layer close to a rotating disk. The case in which the solid 
surface is at rest and the liquid rotating, or swirling in a free vertex, is much more 
complicated, and the present analysis is only an approximate one based on Pohl- 
hausen’s use of the momentum integrals through the boundary layer. The growth 
of the boundary layer from the edge of a cone immersed in a fluid swirling in a free 
vertex is discussed, and the velocity with which the fluid in the layer flows towards 
the vertex is calculated. 

It is found that even when the fluid concerned is water with kinematic viscosity 
v= 0-01 and when it is forced through an orifice 2 mm. diameter by a pressure of 
150 Ib. sq. in., a large part of the efflux reaches the orifice by way of the boundary 
layer. 


Statement of simplified problem 


TE swirling jet here considered is one in which liquid enters tangentially 
a chamber of radius R, and emerges through an orifice of radius R, which 


is considerably less than R,, after passing down a converging cone whose 
vertex subtends an angle 2a. If the entry conditions are suitably designed, 
the longitudinal component of velocity is small compared with the swirl 
component except at points close to the orifice. In order to get some idea 
of the thickness of the boundary layer at the surface of the cone, and the 
distribution of velocity in it, an attempt is made to solve a simplified 
problem in which the longitudinal or axial component of velocity outside 

* This paper was issued as a report to the Aeronautical Research Committee in March 
1945, 

[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1950)] 
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130 G. I. TAYLOR 

the boundary layer is neglected, only the swirl being considered. The main 
body of the liquid is, in fact, taken to be moving with velocity Q/r in circles 
round the axis, where r is the distance from the axis, and Q is the circula- 
tion constant. Within the boundary layer the retarded liquid is slowed 
down by viscosity and consequently has not sufficient centrifugal accelera- 
tion to hold it in a circular path against the inward radial pressure gradient, 
which must exist throughout the fluid in order to balance the swirling 
motion outside the boundary layer. A current directed towards the point 
of the cone must therefore set in over a layer close to the surface of the 
cone. This current must of course be removed before the point of the 
cone is reached if the motion is to remain steady. 


Boundary-layer equations for viscous fluid at the inner surface of 
a cone 
The full equations for the motion of a viscous liquid which flows sym- 
metrically with respect to the axis @ =: 0 are, in spherical polar coordi- 
nates (1), 


oe vey uw v? lap | v2 2u 2 @v = 2qecoté (1) 
iene — a —_ — a P08 cums aaa ame anne oe 
aR ' Reo R R peR- Rk? R* 06 R 
ye 8 Ow 3 we, ew cot @ vn w (2) 
— —_ +- =v “W — - ° pA 
oR’ Reo’ RR’ R RR? sin? 


Here R is the distance of a point from the origin which is taken at the 
vertex of the cone, 

6 is the angle a radius vector makes with the axis, 

u is the radial velocity, 

w is the component of velocity perpendicular to axial planes, 

v is the component in the axial plane perpendicular to R.) Since + 
is positive when the flow is in the direction of increasing @ it is 
directed towards the surface of the cone in the boundary layer, 

is the pressure, 


~ 


I 
p is the density, 
v (= p/p) is the kinematic viscosity, 
and the coordinate system is shown in Fig. 1. 
It may now beassumed that the boundary layer is thin and that variation 
in p through its thickness may be neglected. Since outside the boundary 
layer QO 


= <— (3) 
Rsin @ 


w 
the pressure is given by 


1 QQ? 
Pp _ 3 gp + Constant, 
p a sin* 
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1 ép Q2 
pc R R3 sin26" 


It will be found that if the boundary layer is thin compared with R the 


so that (4) 


following terms may be neglected. In (1) v?/R compared with w?/R, and 
ou 2cv. 2veoté 


A compared with V2u which itself may be represented 
R2' R206 ~— OR? ; 





Limit pf Pohlhausen- 
boundary layer at 7=1 


; l cu 
by its largest term Re : 3° In (2), (ww cot 6)/R may be neglected compared 
2 HA2 


with the other three terms on the left-hand side and (w/ R? sin?@) compared 


“9 
“— ; ' Cw 
with Vw which may be written — —. 
} R? 06? 
The approximate forms of (1) and (2) are now 
ou. v Ou Q)? we v mu (5) 
u + + he er * ee thee eer e 
eR’ Rob Rsin?@ RR? 66?’ 
ow v ew wu v Ow , 
u- : : + , oo ee 2° (6) 
eR’ RebO’ R_~ RO 
The full equation of continuity is 
ou  2u la, v ” 
| __ cot @ = 0, (7) 


oR' R'RO'R 

but the term v cot @/R is small compared with the other three so that the 
approximate form which will be used is 

cu , 2u, 1 dv 


eR' R'RO 


0. (8) 
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It seems unlikely that a full solution of the equations (5), (6), and (8) could 
be obtained, but it has been found in the less complicated problems of 
two-dimensional boundary-layer flow near a cylindrical aerofoil that the 
approximate method introduced by Pohlhausen gives results which agree 
well with measurements, particularly, as in the present case, when the 
flow in the layer is in the direction of the pressure gradient. To apply this 
method the momentum integrals through the boundary layer must first 
be developed in a form suitable for application to a fluid circulating round 
the inside of a conical surface. The Pohlhausen method then consists in 
assuming an arbitrary form for the distribution of velocity in the boundary 
layer which satisfies the boundary conditions at the solid surface and at 
the outside of the boundary layer, and also satisfies the momentum 


integrals. 


Boundary-layer momentum integrals 

To obtain these integrals we may integrate the equations of motion 
(5) and (6) through the boundary layer, which must be assumed to have 
a definite thickness 5. If 2a is the vertex angle, variations in sin @ in the 
boundary layer may be neglected and sin 6 taken as sina. (5) then becomes 


cu cr vw Cu 


F P ()? ue v [cu 
ial dO +. | Td Fe 9 
| "or@+) Rw (ame R) Foil, a; 


R is constant over the path of integration and the integrals extend from 
6 r—86/R to 0 x. The second term in (9) can be transformed using 
(8), giving 


a 


v CU a 0 &. ov ] ‘ Pi du .. 20% 
IR =p - | aso "| dé — ple lant | (w’ at 4 dé. 
(10 


In the present problem w vanishes at 6 = « and at @ = a—8/R, so that 


. X } 
[we], -sie = 0 (38) 
and (9) becomes 


a 


cu a i Q)? w v Jou ; 

2 | u—dé+2 | d@- | ———— a dj = —, . (12) 
J éR JR J \Rsin’a R R?| 6} ,_, 

This may be called the longitudinal momentum integral. In the same way 

the transverse momentum integral can be obtained by integrating 

(6) through the boundary layer and again using (8). The result is 


r uw v low 
») dd+-3 -dé = — , 13) 
(uw) ¢ R , ( 


l ow re 
RY 4a-8R1 | eR R?| 60 


Here a difference between (12) and (13) may be noticed, for though wv = 9 
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it 0 x, at 0 .—5/R this is not the case, in fact 


Q) 


(14) 


(UW) y-3/R ()q—s/r- 


? sin a 
To calculate the value of v at the outer edge of the boundary layer (8) may 
used, thus 


(v) s/n ( eh ten) dé, (15) 


Cc 


3 that (13) becomes 


QO cu -_ uw v [ow 
m on R 1 9u4) dé ») dd+-3 | 10 - . 
R? sin 6 | | oR | | ep) R aha), x 


— (13.1) 
{pplication of the Pohlhausen method 

To apply the Pohlhausen method (2) it is necessary to assume an 

bitrary form for the distribution of velocity through the boundary layer. 


The following forms are the simplest power series expressions which 


satisfy the required boundary conditions, namely, u = w= 0 at 0= « 
; v e 
nd u 0, w Q/(R sin a= 0, <= 0 at 0 .—6/R. 
c6 fale] 
QE QE e : 
u : }() : n—27 n°), (16 
R sin a‘ Rsin « | ) 
Q) Q (17) 
w a"? *__ (2n—n? 7 
Rsina’ ° Rsin « —7) 
: R(xa—6 
nere 2 5 ) (18) 
c 


s0 that » varies from 0 to 1-0 through the boundary layer. With these 
simple assumptions for the functions f(7) and ¢(7) two variable quantities 
E and § determine the boundary layer at any point, so that E and 6 are 
functions of R only and can therefore be determined by using the momen- 
tum integrals (12) and (13.1). In applying (16) and (17) to (12) and (13.1), R 
s constant, and so may be removed outside the integral sign where it 
ccurs as a factor, but the values of éu/é@R and é(uw)/eR which occur in 
12) and (13.1) are differentials keeping 6 constant so that they must first be 
expressed in terms of the variables 5, £, and », thus 
sin a/cu L me f cE E ef (: n\. 
, ), 


E 4 
Q \éeR RI ROR Ren\eR 


from (18) } r— 0 R(a- 6) ds _ n ds 
F R A Oo 82 dR ~ R : 5 dR’ 


(7 dE al nal! t TR) Tae: (19) 
: dy 


0 that sin “{: >| £. . 
Q \cR RdR R R? 5 dR 
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Similarly 


Fs) 


eR 


— 
SIn“a 


2 (uw) 








( 1 dE 2E fb E R ds 
9 \R?dR_ R2/*" "RA 8dR 


ee 
)n (fh). (20) 
dy 

Inserting these expressions in (12), (13.1), and (15), changing the variable 
of integration from dé to —(8/R) dy and noticing that 

: if : ; ] : 

of J dy —} | f?dy and n (fb) dy - (fd) dy 

J ~ dy a J dy* 

0 0 0 0 


since f vanishes at both limits of integration, (12) becomes 


1 1 
§2 dE? RE*d82 _\ fC... iad sin a/df 
Ae — 2 2 dy -t-j — 2 =— —F - 
(® dR 28 dR ) | fray | ¢ a] i (is) ‘ 
0 0 


(21 
and (13.1) becomes after some reduction 
32/ dE? REP Is 
1 6?/_dE? 2 dd? i TP i vsin a/¢ 
R sollte . ase | 1 ii 1 — an E . 
> | dR’ & Al jie | fa dn) Q aa). : 
; 0 0 
(22 


To reduce (21) and (22) to non-dimensional form the length along the 
generator from the vertex to the point at which the boundary layer may 
be supposed to start (i.e. the edge of the base of the cone) may be taken 
as R,. The non-dimensional variables are then 


ae 
1 R, | (23) 
5,2 i Q )J 7 
Ry A v sin a 
and (21) becomes 
(dE? 1e2 ds ENC, = laf 
al (a : a aR, z) Peat Ri _ - 2] me ai, . 
: ' (24) 


(22) becomes 


, 1 1 
sof CH? K? dt ba fee re _ER " 95) 
131 "3 dR, | fay | ay dn) -0 : 

0 0 


So far the only features of the arbitrary functions f and ¢ which have 
been used are the boundary conditions f = ¢ = 0 at yn = 0,f = 0,¢=!1 
at 7 = 1. The particular forms for f and ¢ given in (16) and (17) will 
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) 


These functions are plotted in Fig. 2. 


0-009524, 


1 1 
| fo dy = 0-0500, 


0 0 


1-0. (7) 
dy n=O 


0-5333, 


In this case 


[ f dy = 0-08333, 
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.2. Assumed distribution of velocity in boundary layer. 


Using these values (24) and (25) become 


. ds E82 00-4666 8? 
OF ~- 5 1~ T 7 
‘dR, |” dR, R, — 0-009524 R, 
dE? , dd? , 
oF ae he 1L20E R,, 
d R, d R, 
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These equations are now in a form suitable for numerical computation 
using E£? and E8} as variables. The calculation has been carried out 


starting with E? E83 Oat R, 1-0. The first approximation to the 
solution in the immediate neighbourhood of R, = 1 is 
E* = 19-1(1—R,) |) sn 
E® —80(1— R,) J 


This approximation has been assumed to hold over the interval R, = 1-0 
to R, = 0-995 and the values of dE?/d R,, d( £8?)/dR, calculated from (27) 


and (28) at R, 0-995. These values have been used to calculate F? 


and Hd; at R, 0-990 and the process continued to R, 0-940. At this | 


stage the interval for each step was lengthened from 0-005 to 0-01. The 
range of intervals given below. 

1-00 > R > 0-940, interval in R, was 0-005 

0-94 > R > 0-90, interval in R, was 0-010 

0-90 > R > 0-80, interval in R, was 0-020 

0-80 > R > 0-74, interval in R, was 0-030 

0-74 > R > 0-70, interval in R, was 0-040 

0-70 > R > 0-05, interval in R, was 0-050. 
The results of the calculation are given in Table 1 and are shown graphi- 


cally in Figs. 3 and 4. Fig. 3 shows 6, as a function of R,, Fig. 4 shows the 


TABLE | 


950 0°940 0°93 o*92 o'go0 


09 | 0°442 5 0625 | 0°764 o'884 0°990 1'093 1185 1'273 1°432 


I 5 I 73 1°87 7 2-046 2°115 2 2°274 
- -) Sm°5 75°5 751 72°¢ 690°1 66-2 63°7 61°3 59°4 57°5 
o'77 7 - 0°65 0°60 55 
f 72 1°55 I 7 2°283 2°46 \ 2°960 3°24 3°515 
15 7 2 2°595 2°627 652 2°665 2°658 2°6 2°60 
7 41°2 39°1 6°7 33°9 57 29°7 77 
5 5 or25 2 I o’'10 o'o5 
[ 4°38 69 2 5°35 5°7 6:06 6°45 6°81 
C ) 2 2°2 2°19 Io 2°00 Ig! 
61 24°5 3 7 20°5 19 2 17°2 16°4 


angle y which the direction of flow over the surface of the cone makes with 
its generators so that y is the deviation of the direction of the surface drag 


from that of the pressure gradient. y is calculated from the formula 


cot x = cu = = LE, (30) 
; Cn! Cn] y=0 fi 
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X 1s the angle which the surface drag makes with the generators of the cone. 


138 G. I. TAYLOR 


Paths of particles on surface of cone 
The paths of particles lying on the surface of the cone can be caleu- 
lated from the expression 


| l "| _1iR 
R, do, surface path ; 
Ry 
n' 9 
so that 6,—8, oa 4 (31) 
1 


where F,,6, are the polar coordinates of a point on the developed surface | 


of the cone, and @, is the value of 6, at the beginning of the path. The 


Direction of current 
1-0 outside boundary layer 
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Fic. 5. Developed surface of cone showing path of particle swept over the 
surface by the current in the boundary layer. 


path obtained in this way is shown in Fig. 5, and since a plane disk is one 
member of the family of cones considered, i.e. the cone for which a = }r, 
the curve of Fig. 5 represents the actual calculated path of a particle lying 
on the bottom of a cylindrical basin when water is fed in at the circum- 
ference with a swirl, and drained off at a hole in the middle. 


Application to swirl atomizers 

The flow in a swirl atomizer differs from that considered in the preceding 
analysis because in the atomizer the velocity outside the boundary layer 
has a component parallel to the generators. In most atomizers the tan- 
gential component of velocity is so large that the component parallel to the 
generators is comparatively small, except quite close to the orifice. Over 


the region where the boundary layer is being built up the approximation 
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here used is probably sufficiently accurate, so that the calculated boun- 
dary-layer thickness 6,, say, at the radius of the orifice should certainly 
provide a sound basis for estimating whether the greater part of the 
outflow enters the orifice from the boundary layer or from the main body 
f the swirling fluid. 

If the radius of the orifice is R,, and the pressure which drives the liquid 
through it is P, then Q is less than R, U where 


1/9 P\ 
U es ) (32) 
N\p 


if R, is the radius of the conical chamber at the point where the liquid 
enters, R,/R, may be taken as the value of R, in Table 1. Thus if 
R, R, =: 1/10 it will be seen from Table 1 that the corresponding value 
f, is 2-0. The value of R, at the orifice, as defined in the analysis, refers 
to distance from the vertex, while R, refers to distance from the axis, so 
that R R, cosec x and 

r) y ee ). (33) 

R, R,sina\ Q 


und in the case where R,/R, 0-1 


° 20-0(si [”\3 
R, 2 (sin a)~# ol] ° 


If, for instance, the vertical angle of the cone is 90°, a 45°, and 


re) ly 5 Vv 
23°8 | so that > 23°8 ——— 5 (34) 
R, AQ R, ° UR, 
For water v = 0-01. If therefore water is driven through a nozzle 2 mm. 
diameter at 10 at mospheres pressure U 45 metres/second so that 
5 0-01 r 
. 23-8 |} a OF (35) 
R, A \(0-1)(4:5 x 108)} 


In this case, therefore, the boundary layer extends inwards to a distance 
one-ninth the radius of the orifice from its surface. When it is remembered 
that in actual swirl atomizers a hollow core runs through the orifice 
occupying a considerable proportion of its area it will be seen that much 
of the water reaches the orifice by way of the boundary layer. If the 
viscosity had been that of, say, olive oil, i.e. v = 1-0 and the operating 
pressure 15 lb. sq. inch the boundary layer would have extended across 


the whole area of the orifice even if there were no core. 


REFERENCES 
1, 8. GoLpsTEIN, Modern Deve lopments in Fluid Dynamics, vol. 1 (Oxford, 1938), 105. 
2. K. POHLHAUSEN, Zeits. f. angew. Math. und Mech. 1 (1921), 257. 





THE STEADY FLOW OF VISCOUS FLUID PAST A 
SPHERE AND CIRCULAR CYLINDER AT SMALL 
REYNOLDS NUMBERS 
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SUMMARY 


In Part I of this paper the flow patterns around a sphere are first computed in 
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detail on the basis of Goldstein’s exact analytical solution of Oseen’s linearized 
equations of motion for the steady flow of an incompressible viscous fluid past the 
sphere, and it is found that, in accordance with the results of observation, a stationary 
vortex-ring is formed behind the sphere. Further, by carrying out numerical calcula- 
tions for small Reynolds numbers, it is found that, even when the Reynolds number 
assumes very small values such as 0-1, the stationary vortex-ring, though of very 
weak strength, is still formed behind the sphere. 

Next, discussions on the drag on the sphere are made, with special reference to the 


pressure drag and the frictional drag separately. It is thus found that, as far as the } 


calculations are made on the basis of Oseen’s equations, the pressure and frictional 
drags contribute to the total drag on the sphere in the ratio 1:2 for any value of 
the Reynolds number. 

Similar discussions are made in Part II for the case of steady flow of a viscous 
fluid past a circular cylinder, and it is thus found that, in accordance with the 
experimental results, the two standing eddies are always formed behind the cylinder. 
Further, it is found that such standing eddies are still formed even when the Reynolds 
number assumes very small values such as 0-05, 0-1 and 0-25, in contradiction to 
the common view that for very small Reynolds numbers, in the neighbourhood of 
unity, the two standing eddies are not formed. 

The drag experienced by the circular cylinder is also discussed, with special 
reference to the pressure drag and the frictional drag separately, and it is found that, 
as far as the calculations are based upon Oseen’s equations of motion, the total drag 
on the circular cylinder is divided, in exactly the same proportion, into the pressur 
and frictional drags for any value of the Reynolds number. 


1. Introduction 

As is well known, the exact analytical solution of Oseen’s linearized 
equations of motion for the case of steady flow of an incompressible 
viscous fluid past a sphere was obtained as early as 1929 by Goldstein 
(1) who also obtained the general expression for the total drag on the 


sphere. For the case of a circular cylinder, on the other hand, Faxén (2) | 


has given the exact analytical solution of the same equations. 


So far as we are aware, however, no detailed theoretical discussions on 
5 


the flow patterns around these solid obstacles seem to have been made on 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1950)] 
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| 
| 
| the basis of Oseen’s equations. In view of the fact that at small Reynolds 
| pumbers the values of the drags on the sphere and the circular cylinder 
ilculated by using the existing approximate formulae? are in good agree- 





ent with those found by experiments, it seems worth while to discuss, 





n the basis of Oseen’s equations of motion, whether the actual flow 
} yatterns observed in experiments can be obtained theoretically or not, 
nd the first object of the present paper lies in such a discussion. 
Detailed discussions on the drags on the bodies have also been made, 
ith special reference to the pressure drag and the frictional drag 
eparately, and interesting theorems concerning these two kinds of drag 
ted in | baVe been found. 
arizeq | 10 Part I the results of our investigation on the flow past the sphere are 
st the | given. For convenience the essential parts of Goldstein’s analysis necessary 
onary | for our investigation are first reproduced as briefly as possible, with, 
— | however, some alterations in calculations as well as in notation. For the 
> omy | purpose of computing the flow patterns around the sphere, the general 
expression for the Stokes current function has been newly obtained. The 
to the | expressions for the pressure and frictional drags experienced by the sphere 
s the | have been derived, and by making use of these expressions various con- 
ig lusions have been drawn. 
| In Part II the results of similar investigations for the case of a circular 
scous | cylinder are given. The exact analytical solution of Oseen’s linearized 
ht | equations of motion for the steady flow past a circular cylinder has been 
a ‘ lerived, by an analysis similar to Goldstein’s for the case of the sphere. 


on to | Although our analysis is quite different from Faxén’s, detailed description 
od 


of the analysis is omitted here for brevity, and only the essential parts 


necessary for our investigation are given. 


that 
drag Part |. CASE OF FLOW PAST A SPHERE 
Sssul . ‘ : 

2. The general expression for the Stokes current function 

Let Ox, Oy, Oz be rectangular coordinate-axes having the origin O at 

the centre of a sphere, which is at rest in a steadily running stream of an 
head incompressible viscous fluid of infinite extent. We assume that the fluid 
sible | 2t infinity flows with constant velocity U in the positive direction of the 
tein | 2X8 of x. Then the components of fluid velocity at any point may be 
the | Wttten as (U+-u, v, w), where w, v, w are the components of the velocity 
1 (2 In the case of the sphere, approximate formulae for the drag have been obtained by 

Useen, by Lamb, and by Goldstein (1), on the basis of Oseen’s linearized equations of 
s on tion, while, in the case of the circular cylinder, similar approximate formulae for the 
j ' ¢ drag have been calculated by Lamb (3) and by Bairstow, Cave and Lang (4), on the basis 
e on 


{the same equations 


. 
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of perturbation produced by the presence of the sphere which become Furtl 
vanishingly small everywhere at a great distance from the body. If terms 
of the second and higher order in u, v, w are omitted, we get, from the| 





Navier-Stokes equations of motion, the well-known linearized equations} 


of Oseen, ee where t 
U (u,v,w) = — . (: ; = 5 : ) proven, v,w), | In w 

Ox p\ox cy cz | velocit: 

where p is the difference of the pressure from the hydrostatic pressure,| 4 de 
p the density of the fluid, v the kinematic coefficient of viscosity, and V?} "= % 
stands for the operator ¢?/éa?+ é?/éy?+ @?/éz?. tg suite 


The fluid being assumed to be incompressible, the equation of continuity | sphere 
is P ‘ 
cu ov. dw 
thn em she ee ame, 


Cx Oy 


If we put k = U/2v, these equations are satisfied by 





— .. 1 6x ee. 2 ] ox ee es l ox (1) 
G2 2k Ox dy ' 2k ey cz | 2k éz where 
,é ® 
and P =e pl a (2) | 
Cx 
provided that V2y—2k X = 0, (3) | 
6x 
and V*6 = 0. (4) t 
The components of vorticity (€, , ¢) are expressed in terms of the function and t] 
x in the forms ‘ 2 
& = 0, n= —, c= o%, (5) 
oz oy 
Introducing the spherical polar coordinates (7, 0, w), the appropriate general 
solution of equation (3) for the case of the sphere is given by 
- | accen 
E 4 l ens > Bn Xm(€) Pre), that t 
m= 
where the B,,’s are constants of integration andt f 
€ = fr, pb = cos8, . 
while the function x,,(€) is expressed in terms of the modified Bessel 
function by n\} 
xn(€) = 2m-+1)(F}* Kyl O} 


and P,,(u) denotes the Legendre function. 


m 


+ € has been already used above to denote a vorticity component, but since it was equal ' 
to zero its use now to denote kr can cause no confusion. 
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Further, the appropriate general solution of equation (4) is given by 


d 1 ~ A, onl) (7) 
— = 


n=0 
svhere the A,,’s are also constants of integration. 
In what follows we shall use, in place of u, v, w, the components of the 
elocity of perturbation along r increasing, @ increasing, and w increasing, 


ind denote them by v,, vg, and v,, respectively. Then it is evident that 


, 
0, and, after some calculations, we get the expressions for v, and 
, suitable for fulfilling the boundary conditions at the surface of the 


sphere as well as for dealing with the Stokes current function, as follows: 


m 1 4 {(n t ua > By, ©, ,(é)! P(e ), ] 





. =e 
; p (8) 
v9 U sin @ > | : 7" : S B,,4 “an(é)| P(e), | 
ae" = ae 
where 
P { mm P m+] | 
Wy nats) | 2m l Xm i(é is Im- g Xm (8) Snnl&) 
oe ee (9) 
Xm \2m = i m-—1,n\%S 2m+ 1° m+1,n i’ r . 
A ¢ | I ¢ I \, 9 ’ 
mn (§) 2 Xm ilé ) 5 ‘ Xm+1(€) { fm.n(€)— 2G m,n(E)}; 
{\2m—1] 2m+3 | 





4 


and g,,,, are the coefficients in the expansions 


m,t 


HP (KH) = & SmnlE)Pr(t)» 
0 


n 


ert P" (us) > Imn(€)Pr(e): 
1 


n 


, accents denote differentiation with respect to ». It can readily be shown 


(2) 
(3) 
(4) 
ction 
(5 
neral 
essel 
(6) 
equal 





that these functions f,,, ,(€) and g 


(€) can be expressed as: 


myn myn 


fh (é) = (Qn1) ¥ (2A) Cm —2))! (2n—2))! 
. font (j!)?(2m+2n--27+1)! ~ 
| (m+n—))! \? } 
\(m—j)! (n—j)!J Pmin-2j(§)> | 


m—1 ° ° mI (10) 
4(2n+-1) % (27+ 1)! (2m—-2j7—1)! (2n—2j7—1)! 


In At . ) 
ni) n(n+1) & (j!)?(2m-+ 2n—2))! “ 
(m+n—)j)! \?-bm+n—25-1(€) 





\(m—j—1)!(n—j—1)!f é 
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~y 
where @,(é) = (2n r(Z)z 4(&). (11 
es 
Using (6), (10), and (11), the functions ®,, ,(€) and ¥’,,,(€) can be 
expressed in terms of the modified Bessel functions of half-integers, 


I,.\(€) and K,, .,(€), and consequently they can ultimately be expressed 


in finite series of their argument €. Some of them, which will be used’ 


later, are given here. Thus, 


p 377 we 1 
Moi(E) = — {f-- oa wi | 
oe ee &} 
(12 
Dar 6.9 6 3 6 
D, o(€) —~ A = of ef _  — 
° a. £ et , #7 
In particular, it is easily proved that 
(2m—-1)z : 
®,,, o(€) ‘ = 2¢2 ° (13) 


Next the constants of integration A, and B,, contained in the above | 


general solution can be determined from the boundary conditions. Since, 
however, the boundary conditions at infinity are automatically satisfied, 
we have only to consider the boundary conditions at the surface of the 
sphere. Denoting the radius of the sphere by a and its diameter by d, the 
said boundary conditions become 


e, LU’ cos @, vg = Usin8, 
atr a. 
Thus, making use of the expressions (8) for v, and vy, these conditions 


vield immediately 


(n- 1)A,, | l N B® (é) (—1l (n 1), | 
es hook 


qnt2 ‘2. ™ 10 (n = 0, 2,% 
m=0 (14 
A lo . l 1), 
3 S BN m n(€o) " 
ant2 2 Z, iii 10 (n ee: Se 


m=0 
where we have written ka == €, which is connected with the Revnolds 
number Rk = Ud/v by the relation 


Elimination of the A,,’s between the above two equations yields 


> B Ann (€o) = — 2 l), (15) 
—_—” 10 (n= 2, 3....), 
W here A... (Ee) ®.n(Eo) T (n _ ] )Fwn.n(Ge)> 


and the B,,’s can be determined by solving these simultaneous equations. 
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The expressions for A,, ,,(€)) corresponding to m = 0, 1, 2 andn = 1, 2,3 
wre given in Goldstein’s paper (1).+ Only the expression for Ay ,(&>), which 
will be used later, will be given here. Thus, 
Noa(Eo) 4 Te “(F + :)} (16) 
4g | £ fy &/) 

Next, by making use of the preceding exact solution of Oseen’s equa- 
tions, we shall proceed to the calculation of the Stokes current function ¢ 
for the flow under consideration which is connected with the components 
f velocity v, and v,, along and perpendicular to the x-axis respectively 
by the relations 

l cd l ob 
ek on wm Cx’ 
where wa = rsin@. Then the radial component of velocity v, is expressed 


. . e 
in terms of & as ; 
l eu 


rou 
If use is made of the expression for v, as given by the first of equations 


(8), we have 


y/ 2 © 

cus . < {(n 1)A,, Fe oe | . 

CU > | r T 2 z= Bn Pnu(8); P(e), 
2 n=0 m=0 


which, by the first of equations (14), is further transformed into 


C us Ua 


Pw tau SS Bylo) (Eo) —P° Onn (€)}P,(u)- 


Cy ? ha ha 1 — | 
m U0 n 0 
Making use of the relation 
2 c 2 < 
a ®,, o\So) / ®,,.0(€) Q, 


which follows immediately from (13), as well as the formula 


] , / 
P(e) Gn J pen (4) PrsalH)}; 


ait 
we integrate the above equation with respect to p. We have then 
‘ Ua P,(2) 
r 3 


iuUS SB — “@,., (9) —7° (e)\** fp alin 
= se m m > m,n on rw 1 


n On 1 
where f(r) is for the moment an arbitrary function of r alone. 
If, as usual, we denote by u& 0 the particular stream-line which is 
It is to be remarked here that there is a difference in the order of suffixes m and n 


tween Goldstein’s A, (é,) and ours. Thus our A, 9, for instance, is equivalent to Gold- 
Stein's A, 1 


be 
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composed of the two radial lines 6 = 0 and 6 = z, the function f(r) caneasily | as the 
be determined as Vai Va3 of the 
= 3r ee 3) F(x). functi 
Thus we have finally Thus, 
; is uns 
| _ Ua? Py(p)— P,(u) 
7 > 3 = terms 
a appro 
qn+2 Fe ust (2) 
] 2 -1(H)—F, 41 (¢ - 
-3l > B,,| 9 ®,,,.n(&)- r ®,, (g)| 9 » (I ‘) | 
\ | 2n-+ ‘ b 
m=0 n=1 | 1 
and this gives the current function for the flow produced by the steady 
translation of the sphere through an infinite mass of viscous fluid at rest | 
at a great distance. 
The corresponding current function for the steady flow past a fixed Fre 
: or . . . . . Dow 
sphere in an infinitely extended fluid can easily be obtained by adding to veyn 
the right-hand side of (17) the current function for the uniform flow, | ct 
which is given by 
b = 4Ur? sin? = 4Ur2{P,(u)—Py(u)}. 
Thus the required general expression for the current function becomes | Hence 
ultimately furth 
|, @\Piin)—Pilpz 
ub Ul r2—— ol) - o(f) 
r 3 
ee: aoe On 
are ja .2 (é) \P, 1 (H) 
-$1 2 D3 Bm) yn ®,,,(&)- U Dn» ( n gf) f - a 2 (18) appre 
m=0 n=1 | 
The first term on the right-hand side represents the current function for 
the steady flow of a perfect fluid past a fixed sphere and consequently the 
second term represents the effect of viscosity of the fluid upon the current 
function. Tk 
If we write further 4, y/(Ua?) and r, r/a, the general expression — expr 
for the non-dimensional current function sb, becomes 
} 2! Po(u)—P,(u) 
By ri— ; — 
ry 3 
and 
x. x 
% ] \\F —P, .4(u) at 
1 __.»2 ~~ n+1\h Q patt 
2 ») p2 B,, re Dnn(Eo) ri ® m, nlEor 1) | “ +] vs ( I It 
m=0n=1 
expr 
3. An approximate expression for the current function Stok 
For fairly small Reynolds numbers, the above expression for y, may be beo 


approximated as follows. By carrying out detailed numerical calculations, 
it has been found that the absolute values of the B,,’s decrease fairly rapidly 
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as the value of m increases and that the value of B, decreases as the value 
of the Reynolds number FR decreases. Also, it is seen that the value of the 
function ®,, ,(€) becomes small when either or both of n and € increase. 
Thus, taking account of the fact that in reality the flow around the sphere 
is unsymmetrical with respect to the plane of yz, we retain only those 
terms corresponding to m = 0, n 1 and m = 0, n = 2, and we have 
ipproximately 


Mo1(o)— r ®o,(€0r:)| sin?@— 


ao Begs l 
wb, ivf sin?@ 1 Bo| 
r, 


l ' _ 
Bo} PoalEo) —17 Dy (Eo rn) sin?@ cos 6. 
1 
From (15) we have B, 6/Ay,(€) approximately. But, when the 


Reynolds number & is sufficiently small, € is also very small, and the 
. . >0 . ? 

function Ap, ,(€) as given by (16) can be expanded as 

aid 2 3 

— 3 2¢% 3). 

Xoa(Eo) ¢ ul aot 30+ O(f5)}- 

$0 

Hence, for sufficiently small Reynolds numbers, the value of By may 


further be approximated to 


2 i 
By — $0 


On the other hand, when both €) and yr, are sufficiently small, we have 
approximately, from (12), 


5) 9 
Dyslé) = —2a(2—1), MaulGors) = —¥e(-= 1}, 
So >0°1 
Do o(Eo) ba, Do (E971) }77. 


Thus, making use of these results altogether, we get an approximate 
expression for ys, in the form 
‘ " ‘ , 
Uh rs } sot ) 4 z(t 3} cos sin2@, (20) 
ry Pe r; 4 ry 
and this may conveniently be used for the computation of the flow 
patterns around the sphere at small Reynolds numbers. 

It seems worth mentioning that in the limit when R > 0, the above 
expression (20) degenerates into the well-known expression (5) for the 
Stokes current function for the viscous flow around the sphere which can 
be obtained on the basis of Stokes’s linearized equations of motion, namely 


l oS ] 

of . % 2 ond 
by -(1 _ 5-3) sin é. 
2\ or, ' 
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4. Computation of the flow patterns 

Using the preceding approximate expression (20) for the current fune- 
tion %,, detailed numerical calculations have been carried out for the case 
when & = | and various stream-lines in a plane containing the axis of x 


have been drawn. The results are shown in Fig. 1, from which it will be 


























Fic. 1. Calculated stream-lines past a sphere at R 1. In the half-field of flow abov: 

the central stream-line yf, 0, the stream-lines, shown by full-line curves, correspond 

respectively to ¥, 0-005, 0-05, 0-1, 0-2, 0-3, 0-4, 0-5, ..., while the stream-line shown 
by a dotted-line curve is at ys, —~0-0001. 


seen that, in accordance with the results of observation, a stationary 
vortex-ring is formed bebind the sphere. 

Further numerical calculations have also been carried out for smaller 
teynolds numbers, and it has been found that even when the Reynolds 
number assumes very small values such as 0-1, the stationary vortex-ring, 
though of very weak strength, is still formed behind the sphere. 


5. The pressure and frictional drags on the sphere 

We shall next proceed to the discussion of the pressure and frictional 
drags experienced by the sphere. As is well known, two different methods 
are commonly used in calculating the drag on a solid body immersed 
in a stream of fluid. In one method the drag is obtained by summing 
up the viscous stresses exerted by the fluid on the surface of the body, 
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while in the other it is calculated by applying the theorem of momentum 
to an infinite mass of fluid surrounding the body. The second indirect 
method has been used by Goldstein to calculate the total drag on the 
sphere. Here we shall use the first direct method and derive the general 
formulae for the pressure and frictional drags as well as that for the total 
drag of the body. 

Denoting the outward normal to the surface S of the body by XN, let 
py, be the x-component of the viscous stress at any point on the surface. 
Then the total drag D on the body is given by 


D } Pr, AS, 
where the double integration is taken over the surface S of the body. If 


we denote the direction-cosines of the normal N by (/, m, n), we have 


D | ( lp dS+- pv ( ( (nn—ml) dS. 


On the other hand, it is well known that the total drag D on any solid body 
is in general composed of the pressure drag D, and the frictional drag D,. 
The pressure drag is evidently given by 

D, }] lp dS, 
and thus, by comparing this with the above expression for D, we have 


D, pv | | (ny -ml) ds, 


If use is made of (2), the pressure drag is expressed as 


waa 
D, pl Pads. (21) 
* * Cx 
On the surface of the body we have v = w = 0 and consequently, by 
(1) and (5), the y- and z-components of vorticity on the surface of the 
re) = c 
C2 cy 


body become 


Thus the formula for the frictional drag is transformed into 


pp eds ab\ .. 
D, pl | (m +n—|d8. (22) 
et 4 ' cy Cz 
Further, by adding (21) and (22) side to side, we get the total drag D in 
the form ao 
ed 
} l IS 9° 
3 | | on (23) 


where ¢/¢N denotes differentiation along the outward normal N to S. 
It seems to be worth noticing that all the general formulae for the total 


drag as well as for the pressure and frictional drags contain the harmonic 


{ 
| 


unction ¢ only. Making use of the above general formulae, we shall now 
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proceed to calculate the pressure and frictional drags on the sphere. In 


this case we have 


1 
D, 27a*pU | (1) dy, 
. COX] pa 
1 
1 ad 
D, = —2na*pU | (mir + n° du. 
: cy "4 a 
—j 


After some reductions by the use of formulae for the Legendre function, 
we have from (7), 


Cd) w= A, n+1 
— —l > " f(n-4 > = 12)P .,(u)!, 
(=) a 0 qrt2 2n-t gi" FAH) ” he AH) 


Ch cd ~wt A, (n+1)(n4+2 
(m +n =——l > a ip (u)—P,, .o(h)}. 
a n+2 9 © ” n+2 
cy Col nan a 2n+3 
and putting these expressions into the above two formulae and carrying 
out the integration, we obtain 
4p fJ2, 8 nf 72 
D, == ¢27pU?A,, D, SarpU2A,. 


p 


But, from (13) and the first of equations (14), we have 
. 7a? 

Ay = —> > (2m+1)B,,. 

0 4€2 7 


Consequently we bave ultimately 


"3 & fe 
; a 0 (24) 
477u ES oe 
Dy = 3 > a L1)B,,, 


where » now denotes the coefficient of viscosity of the fluid. 
Further, adding these results side to side, we obtain the expression for 
the total drag D (= D,+ D,) in the form 
D2 Bi, are 
D== = - (2m+1)B,,. (25) 
m 0 
which is in accordance with Goldstein’s result obtained by applying the 
theorem of momentum, since his B,, is equal to our UB... 
Comparing the above expressions (24) for the pressure and frictional 
drags with the expression (25) for the total drag gives immediately 
D, LD, D, sD. (26) 
Thus we arrive at an interesting result that, if the calculations are based 
upon Oseen’s equations of motion, the pressure and frictional drags con- 
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tribute to the total drag of the sphere on the ratio 1:2, whatever the value 
of the Reynolds number may be. This noteworthy result does not seem 
to have been noticed so far. 


Part II. CASE OF FLOW PAST A CIRCULAR CYLINDER 
6. The general expression for the stream function 
The exact solution of Oseen’s equations for the case of steady flow of 
an incompressible fluid past a circular cylinder can be obtained in exactly 
the same manner as in the case of the sphere. Thus Oseen’s equations 
for the two-dimensional steady flow and the equation of continuity are 
satisfied by 


b l < 
u st + ex xX: 
Ox «2k Ox (27 
27) 
f 
s CD \ ] € Xx | 
cy 2k cy 
od 
and p = pl _ , (28) 
Cr 
provided V?x 2h X — 0), (29) 
Cx 
and V°4 0, (30) 
here i: U’/2v as before, and V? now stands for the operator 6?/¢x?-+ c?/ey*. 
The vorticity is expressed in terms of y in the form 
c=. (31) 


cy 
Using the polar coordinates (r,@), the appropriate general solution of 
equation (29) is given by 


xX Uekroos? S BK, (kr)cos m, (32) 
m=0 


vhere the B,,’s are constants of integration and K,,(kr) is the modified 
Bessel function of order m. Also the appropriate solution of (30) is 

= A,, cos nO (33) 
La nm O® 


n=1 


¢ — UA, logr—U 


here the A,,’s are constants of integration. 
In terms of the functions y and ¢, the radial and transverse components 
of the velocity of perturbation v, and vg are given by 
CO l CX , 
v -~x cos 6, | 


. cr 2k cr 


; 


hd ] 4 ° 
‘ “Xx +xsiné. | 


reco 2kr co 
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Thus, substituting the above expressions for y and ¢, we have, after 
some reductions, 
x 


_< cos né .< » f2,R ee 
UY, == > A, pnt — 4 b3 Bin\ Et >, Pmulkr)eos n6 


n=0 m=0 n 


— 


(34) 
— sin n@ a 

Up | A ii > BY. (kr)sin 8, 

ie n yn 1] 4 3 m ) 


mit ( 


m=0n 1 
where for the sake of simplicity we have written 
® (kr) = (Kunst T is WMLn- nt Riel T | 
+K (I 


m\*m—-n-1T a n+17 igie~a4 manele 
m+1 K,,, WU mn- rn ‘. + a T 
ee en on a 


m PY Dot 1] m~—n mn 


t ‘m.n(K7) (K 


m+n-+ 1)» 
the argument of the modified Bessel functions K,, and J,, being kr. 
These expressions are conveniently used in satisfying the boundary con- 
ditions on the surface of the cylinder as well as in dealing with the stream 
function later. 

Denoting the radius of the circular cylinder by a and its diameter by d, 
the boundary conditions on the surface of the cylinder are 

v, U cos 8, vg = Using, 

at r = a. If use is made of the expressions (34) for v, and vg, these con- 
ditions yield immediately 


Ay, S B 0. 


aT 2" | 

A ly l 1), | 
7 = S By, ®,,,.»(€o) ” ” 6 « (36) 

a” 4 a 10 (n By Syne | 

m | 

/ te : l 1), | 

1, 44% BY Pi { (n ) 


a n+l 4 A m min\> 
m=0 


| O (n ae. oe 


where we have put ka = &, which is connected with the Reynolds 
number R = Ud/v by the relation 


& th. 


Elimination of the A,,’s from the second and third of equations (36) gives 
the following equations for determining the B5,,’s: 


2 4 2 1), 
Ratton OO 


a m*myn 


(37) 
m=0 \0 (n =, a 


where A 


m,n (So) Fas nl€o) A, (Eo) T Fn+n(Eo)Kimiis(€o) + 


Tay n (Eo) Ain (Eo) T Biss + n-1(€o)Kim(Eo)- 
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after In order to compute the flow patterns around the cylinder we shall now 
fnd the general expression for the stream function y defined by u = éyb/ey 
nd v éys/ex. Now the radial component of velocity v 


, is expressed 


n terms of % by v, = éys/r 06. Thus if, as usual, we denote by % = 0 the 





articular stream-line containing the radial line @ 0 as its part, we have 








(34) 
0 
us 7 U,. dé. 
0 
Inserting the expression for v, given in (34) and taking account of the 
frst and second relations in (36), we obtain 
Ua? sin é — = { (a\"+1 : ) rsin 6 
9) y : T tl > S B, (“) Dnn(Eo) ®,,,,(kr) (38) 
i fad 2 hee | 7 } n 
m=0 n l : 
This evidently gives us the stream function for the flow produced when 
» Jy | the circular cylinder moves with constant velocity U in the negative 


con. | direction of the axis of x in an infinite mass of a viscous fluid at rest at 
‘eam | infinity. 

The stream function for the flow past the circular cylinder at rest can 

yd. | easily be obtained by adding to the right-hand side of (38) a term Ursin 0, 

which expresses the uniform flow. Thus we have finally the required ex- 


pression for the stream function in the form 


Tayi a= - . “ “ \n+l1 ‘si 0 
con- ds T (r e : sin 6+ i [ * » By (") Din n(Eo) Dp, (kr) 7 = n : 


m=0 n l 


(39) 

The first term on the right-hand side is the stream function in the case 

of the flow of a perfect fluid past the circular cylinder and therefore the 

(36) | second term represents the effect of viscosity of the fluid upon the flow 
pattern. If we write %, = %/(Ua) and r, = r/a, the non-dimensional 


expression for the stream function becomes 


l ~~ sin n6 
af. ‘ : 1 2 nN +1 . ’ 
(r, , }sin 0 4 S S Bi iP inn (Eo) ry Dn n(Eo71)3 .n 
Ids “ m=0n=1 nry 


(40) 


7, An approximate expression for the stream function 


Assuming the Reynolds number R to be sufficiently small, we shall now 


ives : ; ; ‘ ' . 
deduce from the above general expression for the stream function an 
\pproximate expression suitable for the computation of the flow patterns 

(37) round the eylinder. 


lo do this, we have first investigated numerically the values of the B,,,’s 


ind found that, for a certain definite value of the Reynolds number R, the 


absolute value of B,, decreases fairly rapidly as the value of m increases 
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and also that the value of B, decreases as R decreases. Further, it has 
been found that the value of the function ®,,, , (£,7,) decreases as the value 
of n increases and it also decreases as the variable £,7, increases. Thus if. 
in accordance with the fact that in reality the flow pattern around the 
circular cylinder is unsymmetrical with respect to the axis of y, we retain 
only those terms corresponding to m = 0 and n = 1, 2 in the double 


series in (40), we get an approximate expression for y%, in the form 


he: * in 
by (r ~ Jsin + 1B, (Po.1(€o) —77 Por (Eo7)} : oo 
1 


" 
sin 26 
+-{Do o(&p) Do o(Eo71)} Dp? (41) 
as 
If we retain only the first term corresponding to m = 0 and n = 1 in the 


first of equations (37), we get the first approximate expression for B, in 
the form ps 
By = > 9 (42) 
[\(€) A (So) T 1,(€o) Ky (Eo) 





and further, when the Reynolds number R is sufficiently small, 


9 


By  » (43) 
j—y—log Je, 
where y 57721... isthe Eulerian constant. Further, when the Reynolds 
number R is sie small, the values of € and &)7, are also very small, 
and we may write approximately 


Do (Eo) 2(1—y log 3&5). Do (E571) 2(1—y—log ! £97); 
Do o(Eo) €o(3 —y—log 3€,), Do (E97) £971 (3 —y —log $ £57). 


Thus if we use the value (43) for By we have finally 


: ~« = . : 
sy I 4 ni . Br, log r,|sin LIC () ri—;] — Drj log r,|sin 20, 
| Py r 
(44 


where 4, B, C, and D are all functions of the Reynolds number FP such that 


ww 


bail C k ' 


D BC 


wh 


2(log R—2-00223)’ 16’ 
When the value of R is very small, this approximate expression cat 


represent sufficiently accurately the flow patterns around the cylinder. 


8. Computation of the flow patterns 
Using the above approximate expression (44) for the stream function yj, 
detailed numerical calculations have been carried out for the case R = 0°25 
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and various stream-lines around the cylinder have been drawn. The flow 
pattern thus obtained is shown in Fig. 2. Similar numerical calculations 
have also been worked out for the case R = 4 by using the exact formula 
40) for &, and the flow pattern as shown in Fig. 3 has been obtained. 

From these figures it will be seen that in each case there are two standing 
eddies behind the cylinder. These calculated flow patterns are seen to be 
in very good agreement with those found in experiment.+ 



































Fic. 2. Caleulated stream-lines past a circular cylinder at R = 0-25. In the half-field 
of flow above the central stream-line yw, 0, the stream-lines, shown by full-line 
curves, correspond respectively to us, 0-005, 0-025, 0-1, 0-2, 0-3, 0-4, 0-5,.... while 


stream-lines, shown by dotted-line curves, in the region of eddies are at 
us, 0-001 and — 0-0025 respectively. 


We have further computed the flow patterns for two other cases in 
which the Reynolds number R# assumes very small values 0-05 and 0-1 
respectively, and it has been found that, as in the former cases, the two 
sfanding eddies, though of very weak strength, are still formed behind the 
cylinder. The regions of the eddies in these cases are shown in Fig. 4. 
where similar patterns for the preceding two cases in which R is equal to 
0:25 and 4 respectively are also shown for comparison. 

Fig. 2 should be compared with a photograph for the case R 0-25 given in Tafel 23 

f Prandtl-Tietjens, Hydro- und Aeromechanik, Band 2 (1931), while Fig. 3 should be 


mpared with a photograph for the case R 3-9 given in Plate 31 of Goldstein, Modern 


Jerse opments in Fluid Dynamics, vol 2 (1938). 
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It is worth noticing that, contrary to the common view,+ the two| _— 
standing eddies are always formed behind the cylinder even when the rag ¥ 
Reynolds number R is very small. first, pI 

stresse 
FZ 
where 
and ¢ 
If, « 
} the th 
evlind 
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sreat 
outwe 
ploye 
In: 
a 
Fic. 3. Calculated stream-lines past a circular cylinder at R 4. In the half-field : 
of flow above the central stream-line y, 0, the eoyoy shown by full-line while 
curves, correspond respectively to uJ, 0-005, 0-025, 0-1, 0-2, 0-3, 0-4, 0-5...., while 
the stream-lines, shown by dotted-line curves, in the pr of eddies are at 
wb, 0-001 and —0-0025 respectively. 
— 
d, : 
eylin 
whic 
R=0-05 R=0: R=0:25 - form 
Fic. 4. 
9. The total drag on the circular cylinder wher 

We shall next proceed to the discussion of the drag experienced by the | the « 
circular cylinder. As mentioned in section 5, two different methods are only 
commonly used in calculating the drag on a body placed in a moving fluid, | mat 
and by carrying out calculations similar to those in the case of a three- 

In Goldstein, Modern Developments in Fluid Dynamics, vol. 1, p. 65, we find: ‘For very 
small Reynolds numbers, in the neighbourhood of unity, the two standing eddies are not whic 


formed, and the Stokes—Oseen approximations apply.’ 
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dimensional body, we can obtain the general formulae for computing the 
rag on a cylindrical body with arbitrary cross-section. Thus if, in the 
frst place, use is made of the direct method of summing up the viscous 
stresses exerted by the fluid upon the surface of the cylinder, the general 
frmula for the drag D is obtained in the form 

oF te, 


cn 


D pU 


. 


(45) 


shere integration is taken round the circumferential curve s of the cylinder 
nd @/én is differentiation along the outward normal vn to s. 
If, on the other hand, use is made of another indirect method of applying 
e theorem of momentum to an infinite mass of fluid surrounding the 
linder, we have * od 
D pU . ilo, (46) 


cn 


here integration is taken round a large closed curve o everywhere at a 

sreat distance from the cylinder and ¢/¢n’ means differentiation along the 
utward normal n’ to o. This indirect method has already been em- 
ployed by Filon (6) and others. 


In the case of the circular cylinder of radius a, the formula (45) becomes 


D pU | (r<) dé. 
. \ cr rea 
vhile the formula (46) takes the form 
”~ / ad 
D ol | (r< °) dé, 
. CF } p+ 


o 
if a large circle concentric with the cylinder is taken as the closed curve co. 
Thus, it is easily found that if we use the expression (33) for the function 
é, the above two formulae give one and the same result for the drag on the 


evlinder, namely D onpU2A 
2n Apo, 


vhich, by the aid of the first of equations (36), may also be written in the 


orm 


D 2apU > B,,, (47) 
m=0 
where x. now denotes the coefficient of viscosity of the fluid. This gives 
he exact expression for the total drag on the circular cylinder. If we retain 
only the first term in the series of (47) and use is made of the first approxi- 
mate value for B, as given by (42), we have 
4apl’ 
I,(ka)Ky(ka)+-1,(ka)K,(ka)’ 


which is the approximate formula given by Bairstow (4). 
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Further, if use is made of the approximate value (43) for By, we get 
Lamb’s well-known formula for the drag, namely 
dap 
1 


J —y—log ba’ 


D 








C 
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Fic. 5. Curves of the total drag coefficient Cp plotted against the Reynolds number R. 
Curves I, II, and III correspond respectively to the first, second, and third approxi- 
mations. Small white circles show Relf’s experimental values. 


Defining the drag coefficient Cp) as C, = D/(pUd), the general formula 
for C, becomes | 
Cy - 2 B,,. (48 
m=0 
Numerical values of C,, have been calculated by using this general 
formula. The results are shown in Table I and Fig. 5. In this table the 


PasLeE I. Values of Ch 


R Ist approx. | 2nd approx. | 3rd approx. 
od 6°98 6°78 6°78 
2 4°46 4°04 4°04 
4 3°60 2°90 2°92 
5 3°26 2:16 2°24 


second column gives the first approximate values calculated by taking only 
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espectively the second and third approximate values obtained by taking 


y 


ps 


espectively the first two and three terms in the said series. In Fig. 5, 
xperimental values obtained by Relf (7) are also shown for comparison. 
] : 


\0. The pressure and frictional drags on the circular cylinder 
—— |} As in the case of the sphere dealt with in Part I, the pressure drag D, 
nd the frictional drag D, on the circular cylinder have been computed 
—— | eparately. In the case of a cylindrical body, the general formulae for D, 
nd D, are given respectively by 

D pU 1 ds, D pU mh. (49) 


4 P Cx f J oy 


here (1, m) are the direction-cosines of the outward normal to the cireum- 
rential curve s of the body. 


In the particular case of the circular cylinder of radius a these become 


— S| D pla od cos 6 dé, 
. « Cx r-a 
0 





— . [ [fe ; 
D, pla ( *) sin 6 dé. 
—— 
Seed J \PY) r=a 
0 
Using the expression (33) for the function ¢ and carrying out the integra- 
—— | tions we have ra 
) D, D, mpU?A,. 
ber R, | further, if use is made of the first of equations (36), we finally obtain 
)proxl- oa o" 
D, D, = mU > B,,. (50) 
m=0 
rmula By adding these two kinds of drag we immediately get again the 
| expression (47) for the total drag D. Thus we find 
(48 D, D, = 4D, (51) 
and these results show that whatever the values of the Reynolds number 
eneral | ° ° ° —" ° 
a | may be, the total drag on the circular cylinder is divided, in exactly the 
le the . 3 


same proportion, into the pressure and frictional drags, on the basis of 
Oseen’s equations of motion. This interesting result does not seem to have 
deen noticed so far. 
| i } et ‘a "7 ) fh cie Y 7 Y 72 "a 
| Ifwe define the pressure drag coefficient ( p, bY Cp, = D,/(pU*d), we 


ive 


Cy z > B,,. (52) 
m=0 


g only ¢ approximate values of C, have been calculated numerically, by 
; ' y, 9) 


is give taining the first three terms in the series > B,,. The results are shown 


m=0 
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in Fig. 6, where similar approximate values of the total drag coefficient C,, 
are also shown for comparison. The experimental values of Cp, obtained 
by Thom (8) are also given by small black circles. It is seen that the 











| | 
"i 2 4 6 8 10 20 30 40 50 6070 
R 
Fic. 6. The dotted-line curve shows the calculated total drag coefficient Cp, whil 
the full-line curve represents the calculated pressure drag coefficient Cp, which is 
exactly equal to one-half of the total drag coefficient Cp for any value of the Reynolds 
number. Small white circles give Relf’s experimental values for the tetal drag 
coefficient, while small black circles show Thom’s experimental values of tl 


pressure drag coefficient. 


calculated valves of Cp, are in good agreement with Thom’s experimenta 
values in the .ange of values up to 10 approximately of the Reynold: 
number. Thus it may be concluded that in the case of steady flow past 
a circular cylinder, the linearized equations of Oseen are applicable appro- 
priately provided that the Reynolds number of the stream assumes values 
less than 10. 
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ON THE TORSION OF A SHAFT WITH KEYWAYS 
By H. OKUBO 

(Inst. High Speed Mech., Tohoku University, Sendai, Japan) 

[Received 15 November 1949] 


SUMMARY 

The torsion problem for a circular cylinder with a number of longitudinal notches 
is solved mathematically. In the particular case of one notch, the relations between 
the maximum stress and the radius of curvature at the bottom of the notch or the 
depth of notch are discussed. As an example, the problem for a circular shaft with 
one keyway in a form for practical use is solved. The results are compared with 
those obtained experimentally or numerically by the previous writers. We find that 
there is a considerable discrepancy between the former and the latter in the results 
obtained theoretically and experimentally. 


1. Introduction 
THE torsion problem for a circular cylinder with longitudinal notches has 





been treated by several writers (1). Gronwall and Weber have both given 
a solution for a circular section with one notch. The former treated the 
ease of a notch in the form of an orthogonal circle and the latter the 
case of a semicircle. Some years ago Filon obtained an approximate solu- 
tion for a circular section with one or two hyperbolic notches. Later, 
Shepherd obtained a solution for a circular section with one notch and 
also one with two notches, using elliptic functions. Recently the writer | 
obtained a solution in a simple form for a circular section with a number 
of nearly semicircular notches using elementary functions. 

But in all the work of these writers the relation between maximum 
stress and the radius of curvature at the bottom of the notch or the 
influence of the depth of notch on the maximum stress has not bee! 


elucidated in detail. 

In this paper the problem for a circular section with a number of notche: 
has been treated, and in the particular case of one notch, the relation 
mentioned above have been investigated numerically. In addition to thi 
calculation we solved, as a numerical example, the torsion problem of 
shaft with one keyway in a form for practical use. 


2. The functions used 
We shall consider a circular cylinder of radius unity with m notches, an¢ | 
take the generators of the cylinder to be parallel to the axis of z. 
Then the stresses which do not vanish are given by (2) 


Ye = LT (°—y), Y. — pT 2-2}, (1 


cy Cx 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1959)] 
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where » and 7 represent the modulus of rigidity and the twist respectively. 


1S The function ¢ satisfies the equation 
24 2d, 
C"“@ C 
= —— = 9, (2) 
Ox“ oy* 
) . . . . 7 
it all points within the bounding curve of the cross-section, and the 
condition 9 ‘ 
d—}(a?+-y?) constant (3) 
otches | at the boundary. 
tweet Let us consider a function 
mi | ] zm 
t wit! J (2) - . (+) 
| wit l amzm qa” gm 
where z = x+/y. Then the real part and the imaginary part of the 
re > . . . ’ . . . r r 
| function f satisfy the differential equation (2), respectively. The real 
j 1 . 
part of the function, expressed in polar coordinates, is 
ri 1—a™r™ cos m0 yor _—_qmy™ cos mé (5) 
es nas Q : - ae oe , ° 0 
? 1+aq2"yp2r__Qqaryrrecosmd a*"+r?™—2a™r™ cos md 
given 
“1 the | The function 4, vanishes on the bounding circle r = 1, and becomes 
er the eas , ’ 2n7 . =e 
"© | infinite at m points, viz. r = a, 6 , 2 = 0, I, 2...., m—1, within the 
» aclu- m 
Later, | bounding circle, when a 1. Putting a 1 in equation (5) it becomes 
h al ] p2m 
dy : : (6) 
writel ? 1+ 72"— 2r™ cos mO 
_ which is the solution obtained by the writer (see 1), and the particular 
case when m | was originally obtained by Weber (3). 
ul ‘ . ° / 
aeaee Combining ¢, and ¢, we have 
or tl : 
; l—a"r™ cos mO p2m—__qmy™ eos mO 
t bee ) wy) A 
| a2 pom Pqinym cos mé qzm | rem Iqinym cos mé 
ai ] pom 
Ott h a B > bd (7) 
lations ]+r2"— 2r™ cos m@ 
to tl Equation (7) contains three unknown constants. A, B, and a, which 
m of can be determined by three conditions concerning the form of the notch, 
is will be seen later. 
3. Form of the notch 
sd Ai For convenience in calculation we shall begin by assuming m = 1, since 
the analysis is similar for other cases. Putting m 1 in equation (7), we 
ve 


l+-a?r2—2arecos@ a*+r*?—2arcosé 


1—arcos@ r2—ar cos 8 ]—r? 
" +B 
] -+-7 
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The function ¢ given by equation (8) satisfies the condition (3) on the 


Cons 
bounding circle r 1, and at the same time on the curve represented by solve 
the following equation, follo 

a?(1+-r?)—ar(1+a?)cos 6 14 B l 0 we Cc 
a ———- A+ Ts=% i 
(1+<a?r?— 2ar cos 6)(a? +r? — 2ar cos 6) 1+r*?—2reos0@ 2 


The above equation contains three unknown constants, and consequently 





Fd | s0 0 


] | valu 


b 





we can specify three conditions regarding the form of the notch. We shall 
consider a notch subtending an angle 2% at the centre, and of depth ¢ and | 
radius of curvature p at the bottom of the notch, as is shown in Fig. | 
Then we can determine the constants A, B, and a in equation (9), so as t 
fit the curve with the above specification concerning the form of the notch 





Consequently, we have the equations 
2a?—a(1- a*)cosB 4 x B ra ie | 
(1+a?— 2a cos 3)? 2(1—cos#) © 2 | 
2 As B 4! — 0, \ (10) 
(a—b)(1—ab) * (1—b)? ' 2 | 
pi | 
P= 5 (@r/de),’ | 
where od ae en oe co ae ae ee, 
dé?}, a As,+ Bk,— b’ lie \(1—ab)? ¢ (a- b)2}’ 
41 Lal 2 2b(1-+L 
" ab | a a ry l ab | a " "ign »( 7 
. \(a—b)® © (1—ab)3} (1—b)? ” (1—b)§ 


We 
t+b 
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nthe | Considerable mathematical difficulties are encountered in attempting to 
ed by | solve directly the simultaneous equations (10), which we solve in the 
following way. For certain values of # and b, if we assume the value of a, 
we can readily find the corresponding value of p from equation (10), and 


(9 


ently G7 * 





“0% 0s Wid We th 


Fic. 2. 


so obtain the curve representing the relation between p and a for specified 


values of # and b. Curves representing this relation, when # = 15°, 
h = 0-625, and # = 12°, b = 0-75, are plotted in Fig. 2. From the curves 
2= C750 2=0§00 2 =Wth0 
b- 0625 b- Wis b =G625 
P=097 P=0l80 P=OlH 
O=/5' 10-/f" =/5" 








(10 


Fic. 3. Forms of notches. 


we can then conversely determine the values of unknown constants A, B, 
and a, for any value of p and subsequently, by using equation (9), the curve 
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€ 


representing the specified boundary can be plotted. Fig. 3 shows the curves 
obtained in this way, when # = 15° and b = 0-625, for several values of p. 

[t is troublesome, however, to plot the curve of the notch in this way. 
We can readily plot an approximate curve by a different method. When 
% and b are given, the points s, ¢, and wu are determined. Describe a circle 
of radius st with its centre s, and a circle of radius p passing the point u 
with its centre on the line su, then the bounding curve is described as the 
envelope of these two circles, as is shown in Fig. 3. 


4. The stresses and the torsional rigidity 
When the unknown constants A, B, and a have been obtained, the 

function is determined explicitly, and consequently we can obtain the 
stresses at any point on the cross-section from equation (1). The stress 
on the x-axis (@ = 0) is 

; 

Y. — PT Aa z+ . 9 
: | (l—ar)? © (a—r)? 


and it attains its maximum at the bottom of the notch. Putting r = 6 in 


l 2B 
Je 2B ay 
(lI—r)? J 
equation (11), we obtain the maximum stress as 
aa — —pr7(As, + Bk,—6). (12) 


The torsional rigidity C is given by (4) 


C= 





| (22 +y?—x . r.. y : “) dady 
} dx = Oy 


. 


n(1— $ vce ds), (13) 


where J is the moment of inertia of the cross-section about the axis of 2, 
s is the line element along the bounding curve, ¢ is the function conjugate 
with 4, viz. 


: ] ] 
us Aa™r™ sin mo! 


\1+a2r2"_ 2a™r™ cos mO ' a2m+-72m— 2amy™ cos mb} 


2Br™ sin mé 


1+-r2™_2r™ cos mb’ 
and the integral is taken over the boundary of the cross-section. 
For the cases when # = 15°, t = 0-375, and # = 12°, t = 0-25, the 
maximum stress is calculated from equation (12) for several values of p, 
and the results are plotted in Fig. 4. 


Assuming the form of the notch to be a semicircle when ¢t = p, and 
combining it with the results represented by the curves in Fig. 4, we obtain 
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Fic. 4. The relation between the maximum stress and the 
radius of curvature p, when # and ¢ are specified. 
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i 
the relations between the maximum stress and the depth of the notch for 
several values of p, and they are shown in Fig. 5. 


5. Numerical example 
As an example of the application of equations (5) and (6) we shall 
introduce an approximate solution for a circular section with one keyway 
in a form for practical use. The dimensions of the keyway are chosen as 
follows (see Fig. 6): width of keyway 0-5, depth of keyway 0-2. 
Combining the solutions ¢, and ¢, given by equations (5) and (6), appro- 
priately, we obtain a function ¢, containing seven unknown constants P,, 
P,, Py, Ps, Py, A,, and Ag, given by 
P,(i—r?) 


1+-r2— 2r cos 0 


‘ l 
( l -?) P | ~ J 5 | a 
—, *\1+7r?—2r cos(0@-4 h;) - 1+r?—2rcos(@—h,)} 
‘ 
3. { ]- a; r cos(@ +-@;) : ] “A; r cos(@—6;) 
— *\1+a7 r?—2a;rcos(8+6;) © 1+a?r?—2a,;r cos(6—8;) 
~ r?—a,r cos(6+-6;) 7 r>—a,;rcos(@—6;) | (14) 


2 2s es 2 . ane ? 
a} +r?—2a;rcos(6+6;) a}+r?—2a,;rcos(@—8;)} 


where h, = 4°, h, = 8°, hg = 12°, h, = 14°, a, = 0-93262, 0, = 12° 23’, 
and a,, 6, are tentatively chosen as 
0-87154, 6, — 13° 16’. 


ay; 


Then the function ¢ vanishes on the bounding circle r = 1, and becomes 
infinite at seven points, marked with an asterisk in Fig. 6. In this case 
the equation representing the curve of the notch becomes 


4 
Fy cs > Pp! an I |, 
1+7r?—2rcos6 — 7\1+7? -2rcos(8+h;) © 1+r?—2rcos(6—h,)} 


j=1 
| . 4/ aj(1+-r*)—aj;r(1+a7)cos(8+6;) 
a, *\{l+ajr?—2a;r cos(0+6;)}{a7 +7? — 2a; r cos(0+8;)} 
a(1 n r?) a;7(1+a?)cos(@ 4;) FE. a (15) 


ce aj r?— 2a, r cos(9@—6,)\{a? + r?- 2a; r cos(d—8;)} J "2 

Let us determine the seven unknown constants so that the curve repre- 

sented by equation (15) passes through the seven points on the bounding 
curve shown in Fig. 6. The coordinates of the points are as follows: 








1°) I 2 3 } 5 6 
; oS | 050196 0°80786 | 081787 0*g0695 0°96302 I 
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Substituting the values of ¢ and @ in the above table into equation (15), and 
| solving the simultaneous equations, we obtain the constants as 





PF 0-0079878. P, 0-0019704, P, = —0-0092601, 
P, 0-0003582. Fr. 0-0000054. A, = —0-0012421, 
| re 0-0017322 











The curve represented by equation (15) is plotted in the right half of 
Fig. 6, and the lines of shearing stress around the notch are shown in 
Fig. 7. The radius p of the fillet in the corner of the keyway is assumed 

) be equal to the radius of a circle passing through the following three 


consecutive points on the curve 


Hence p 0-057. 
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The stress along the bounding curve is expressed by If we 
od te. ilues « 

Te perl (r—' cos # + — - ind), (16) 

| er r er | 

where # is the angle which the radial line makes with the normal to the 
», Con 
The s 
yap fil 





2 Sstressfar uy 








2 


Fic. 7. The lines of shearing stress around the keyway, and the stress distribution 


at the fillet in the corner of the keyway. }7 
curve at the point. The stress attains its maximum at the rounded corne! 
(r = 0-828, 6 = 13° 58’), and its magnitude is Ce 
Tmax = 1:84yr. Fic. 8. ' 
Using equation (13), the corresponding torsional rigidity is found by | tained b 
. ° . ° data, by 
means of numerical integration to be * 
Sons, 
~ § c » In a she 
I = 1-495, if Fas — 0-068, — 
os a solid sk 
Orr, for 
and so C 1-427. ‘ 


The torsional rigidity of a shaft of the same size without keyway is | and jat, 


Cy = 1-571, and hence we have C/C, = 0-908. method 
In the preceding calculation a,, 6, were tentatively assumed to be j maxim) 


a, = 0:87154, 6, = 13° 16’, but if we assume a, = 0°8622, 6, = 14° 26. | Omhas 
we can obtain, by a similar calculation, a solution for a keyway of the ye can 
same dimensions except for the radius of the fillet at the corner. In this | the me’ 


case the corresponding values of C, p, and tTyax are greater 


C = 1-424n, p = 0-034, = 2-18pr. can ase 


o~ 
‘max 
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If we choose a,, 6, aS a, = 0°8885, 6, = 11° 48’, then the corresponding 
ilues of C, p, and 7,,,, become 
16) | ” 
( 1-429p, p 0-089, Tmax 1-68p7. 
the " . P ° ° 
», Comparison of the result with those of the previous writers 
The solution for a hollow shaft, obtained experimentally by means of 
ip films, is included in the well-known paper by Griffith and Taylor (5), 
XMT 
Ai afr 
| 27+ 
S 5 
R /6+ 
Nw 
itiol r 
re! r 
0 i ! i 1 4 1 1 1 1 1 f 
0 G02 LH} Gb G06 41 
Fic. 8. The relations between the maximum stress and the radius of a fillet, ob- 
1 by | tained by several writers: (1) Calculated values, by the author. (2) Experimental 
ta, by Griffith and Taylor, for a hollow shaft, with a keyway of the same dimen- 
sions. The maximum stress, however, is given as a multiple of the maximum stress 
, Inashaft of radius unity without a keyway. (3) Experimental data, by Quest, for 
a solid shaft with a keyway of slightly different dimensions. * Calculated value, by 
Orr, for a solid shaft with a keyway of slightly different dimensions. The radius of 
a fillet is not shown in his paper, and it is estimated from his figures. 
'Y ™ | and later several writers treated similar problems (6). This experimental 
method gives quick results, but there are difficulties in determining the 
O be 


maximum stress accurately when the radius of fillet at the corner is small. 
26’. On 


f the 
this 


has solved the problem for a solid shaft by numerical analysis (7), but 
we cannot expect an accurate determination of the maximum stress by 
the method employed by him. The analysis used here gives a result of 


sleater accuracy, since the solution can be obtained explicitly, and we 


{ 


in ascertain the degree of approximation of the result by investigating 
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how the bounding curve precisely plotted will coincide with the given 
form of the notch. 


We shall now compare the result obtained in the preceding article with 


the experimental result worked out by Quest. The dimensions of notches 
in both cases are not exactly the same, being different in their depths; 
viz. the depth is 0-2 in the former while it is 0-15 in the latter. As is 
shown in Fig. 8, the theoretical values of maximum stress are considerably 
larger than those of the latter, and the discrepancy between the two 
results is too large to be explained as merely the effect of the difference of 
depths. The numerical solution obtained by Orr gives a smaller value of 
the maximum stress as compared with ours, though the dimensions of the 
keyway are slightly different in each case. His result agrees well with the 
experimental data obtained by Quest. 


Acknowledgements 
A grant in aid of scientific research has been given for this work by the 
Education Department of Japan. 


The author takes this opportunity to thank Miss T. Mayanagi, his assis- | 


tant, for her help in numerical calculations in this study. 


REFERENCES 
1. L. N. G. Firton, ‘On the resistance to torsion of certain forms of shafting with 
special reference to the effect of keyways’, Phil. Trans. Roy. Soc. A, 193 (1900), 
309. 
T. H. Gronwa Lt, Trans. Amer. Math. Soc. 20 (1919), 234. 
C. WEBER, Forschungsarbeiten, No. 249 (1921), 31. 
W. M. SHEPHERD, ‘The torsion and flexure of shafting with keyways or cracks 
Proc. Roy. Soc. A, 138 (1932), 607. 
H. Oxuso, ‘Torsion of a circular shaft with a number of longitudinal notches, 
(still in press). 
2. A. E. H. Love, The Mathematical Theory of Elasticity, 4th edn. (Cambridge, 1927 
p- 311. 
3. H. Oxvuso, loc. cit., 
C. WEBER, loc. cit., 31. 
4. A. E. H. Love, loc. cit., p. 312. 
5. A. A. GrirritH and G. I. Taytor, R. and M. 392 (1918). 
6. See The Mechanical Properties of Fluids (New York, 1924), p. 245. Also S. Trmo- 
SHENKO, Strength of Materials, vol. ii, p. 622. 
H. Quest, Ing.-Arch. 4 (1933), 517. 
7. J. Orr, R. and M. 1393 (1930). 














For s 
enters @ 
n term: 
stream- 
tion. V 
experin 

It is 
of appl 
lengths 
conditi 


1. Int 
ANU) 
of me; 
the ve 
Emde 
up to 
He e: 
obser’ 
press 
so th 
gradu 
due t 
num 
found 
terms 


press’ 


K va 
migh 
value 

Pr: 


point 


[Qua 








civen 


- with 
tches 
pths; 
As is 
rably } 

two 
nee of 


lue of | 
f the | 
h a 
Vv the 


ASsis- 


z wit 


1900), 
ackKs 


1927 


Tim 








{ NOTE ON PRANDTL’S FORMULA FOR THE 
WAVE-LENGTH OF A SUPERSONIC GAS JET 


By D. C. PACK 


(University College, Dundee, University of St. Andrews) 


[Received 10 May 1949] 


SUMMARY 

For small fluctuations of a uniform supersonic gas jet about the state in which it 
nters a medium with zero excess-pressure it is possible to express the ‘wave-length’ 

terms of the mean velocity of the jet. The mean velocity is that on the bounding 
stream-lines and not the velocity of efflux as assumed by Prandtl in an earlier investiga - 
m. With this correction, better agreement is obtained between theoretical and 
xperimental values for the wave-length of cylindrical jets. 

It is shown that the approximate theory may be expected to have a wider range 
f application to jets with axial symmetry than to plane jets, and to forecast wave- 
ngths for pressure ratios outside the range for which it may be applied to give 


mditions inside a jet. 


1. Introduction 
A NUMBER of experimental investigations have been made with the object 
of measuring the wave-length of a jet of gas issuing from an orifice with 
the velocity of sound. The first real quantitative work was carried out by 
Emden (1), who used an apparatus which permitted an excess-pressure of 
up to 10 atmospheres in the reservoir from which the jet was discharged. 
He experimented with jets of air, carbon dioxide, and hydrogen, and 
observed the way in which they varied with increasing chamber (reservoir) 
pressure. He noticed that the jet had an approximately periodic structure, 
so that one could speak of a ‘wave-length’, and he speculated that the 
gradual extinction of the jet with increasing distance from the orifice was 
due to vortices drawn into the jet from the external medium. Using a 
number of nozzles of different shapes but convergent to the orifice, he 
found that the wave-length A could be expressed with great accuracy in 
terms of the exit-diameter D of the orifice and the ratio N of chamber 
pressure, pp», to pressure in the external medium, II, by means of the formula 
A= KD\(N—1°9); (1) 
K varied from 0-77 to 1-02 according to the form of the nozzle (and, one 
might add, the range of pressure ratios over which it was used). The mean 
value of K was 0°88. 
Prandtl (4) attempted to account for Emden’s results from a theoretical 
point of view. For jets issuing both from circular orifices and from plane 


(Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1950)] 
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slits, he considered slight perturbations about a steady state, and applied 
the results of his work to obtain a value of K which might be compared 
with that found by Emden. His mean value of K was 1-2, which was too 
high. Rayleigh (6) later derived Prandtl’s formula for the wave-length by 


a different method, but he did not discuss the comparison between theory | 


and experiment. 

In 1941 Hartmann and Lazarus (2) gave an account of a large number 
of measurements of the wave-lengths (among other properties) of jets 
issuing from circular orifices, and compared their results with Prandtl’s 
formula. Their values of A, though slightly higher than Emden’s, were 
still below the mean value computed by Prandtl. 

The discrepancies are due to the fact that the assumptions upon which 
Prandtl based his comparison are wrong. He assumed that the mean state 
in the jet was the state upon exit from the orifice; a study of the conditions 
in a supersonic jet soon leads to the conclusion that this assumption is 
incorrect. The velocity increases after issue and reaches a maximum in 
the centre of the first period (see, for example, (3), section 3); thereafter 
it diminishes and, neglecting the occurrence of shock-waves (which is per- 








. . . a . . ° | 
missible for sufficiently small ratios of excess-pressure on exit to pressure 


in the external medium), it only reaches its initial value again at the end 
of the first period. If a periodic disturbance is assumed, it therefore follows 


that the perturbation is about some state which corresponds to a greater | 


velocity than that with which the jet emerges from the orifice. 
The problem of finding theoretically the wave-length of a supersonic jet 


is reconsidered in the present paper, and a new comparison made with 


experimental and theoretical results. 


2. The oscillations of a supersonic gas jet 
The equation which governs the steady flow of a compressible fluid in 
irrotational motion may be written 


¢ «n\ 0? od 2 nee 2 V 
(u?—a?) "4+ 2uv —— + (v?—a?) __ = §a?-, (2) } 
ex? excy ey? y 
, ; é od 
where ¢ is the velocity potential (so that a 4 v = — are the compo- 
xX cy 


nents of velocity) and a is the velocity of sound at the point (x,y). Also, 
3) 0 for two-dimensional (plane) flow, for which Ox, Oy are perpendicular 


axes; these will be taken in a normal cross-section with the origin O at the 


centre of the orifice; 5 = 1 for flow with circular symmetry; Ox is along 
the axis of symmetry, with O in the orifice as before, and y measures the 


perpendicular distance of a point from this axis. Except where otherwise 
stated, the half-width of the orifice (D/2) will be taken as the unit of length. 
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Suppose that a jet of gas issues from an orifice with a velocity U, equal 


to or greater than that of sound; suppose that, in the plane of the orifice, 
the flow is uniform and parallel to the axis of 2, and that the excess- 
pressure is small enough for the assumption to be made that the fluctua- 
tions about the mean motion are only of the first and higher orders of 
smallness. If we write g? = u?-+-v?, so that q is the resultant velocity at 


ny point, the Mach number of the flow at the point is defined by M = q/a. 
Let U denote the mean velocity of the flow, and let M be the corresponding 
Mach number; then, neglecting quantities of the second order, the equation 


9) becomes 24 no a « 
. 26 8 
Te et i ie (3) 
Cx? oy?’ y cy 
Let d Ux+Cf(y)sin Ba, (4) 
where C is small, and write a? for M?—1. Then 
o7 4.99 2p = 0. (5) 


dy? y dy 
For jets with circular symmetry, 
f(y) = AJ)(aBy)+ BY, (aBy), 
where J, and Y, are the Bessel functions of zero order of the first and second 
kinds respectively. Since v = 0 on y = 0 for all xz, B = 0. 
For a two-dimensional jet, 
T(y) A’ cos aBy+ B’ sin aBy, 
and since v 0 on y 0 for all a, B’ must equal zero. 


The form chosen ford ensures that the condition v = 0on a = Ois satisfied. 


3. Jets with axial symmetry 
Since for a jet with axial symmetry 
d = Ux+AJ,(aBy)sin Ba, 
the velocity components are given by 
u = U+ ABJ,(aBy)cos Bx, 
v A vel, ( vBy)sin Be. 

The boundary of the jet is a surface of constant velocity, and the velocity 
is that which results from the free expansion of the gas at the rim of the 
orifice. The boundary velocity may therefore be obtained by an application 
of the Meyer expansion (see 7), since this velocity and the initial direction of 
flow at the boundary are the same as for a plane jet (this may be proved by 

onsideration of the conditions of compatibility for the characteristics of 
the two motions). Let the boundary velocity be U_, and let the boundary 
lepart from y | only by quantities of the first and higher orders of 
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smallness. The condition for constant velocity on the boundary (to the 
first order of small quantities) may therefore be satisfied only if 


U=U, and JA(af) = 0. (6) 


The mean conditions in the jet are thus the conditions on the boundary, 
This fact was realized by Rayleigh (6), but he made no comment upon its 
significance in the present problem. In what follows U and M refer to the 
boundary conditions. 

Let the successive roots of J,(€) be denoted by &, (r = 1, 2,...); from 
(6), 8 may take any of the values f, = &,/a. The solution of equation 
(3) applicable to axially symmetrical jets is thus: 

: 
¢ = Uxt+ > A, Al(E,y)sin(E, v/a). (7) 


a 
r 


Prandtl and Rayleigh considered that the wave-length / which may be 
expected to be measured in experiments is that given by the first term of 


the infinite series in (7). This gives l 2ra/,, or, introducing the diameter 
for comparison with (1), the wave-length, denoted by A in these units, is 
A = 1D/2 = 1:306D,(M?—1). (8 


Prandtl also took the velocity of efflux to be the mean velocity of the 
jet, but with this assumption (8) gives zero wave-length when the exit 
velocity is sonic—the velocity appropriate to the nozzles used by Emden 
and by Hartmann and Lazarus in their experiments. He therefore com- 
puted the exit velocity for various ratios of reservoir pressure to atmo- 
spheric pressure, making the exit pressure equal to atmospheric pressure 
and using the formula for the wave-length of the jet he found values of A 
which would bring Emden’s formula (1) into conformity with the results 
obtained from his own. For values of NV between 2 and 12 he found values 


of K with mean 1-2; nowhere did the value of K fall below 1-19. On| 


theoretical grounds one would expect zero wave-length for a jet issuing at 
atmospheric pressure, since the ‘periodic’ structure of a jet arises from the 
difference of pressure between the jet and the external medium. For a jet 
in which the pressure on exit is equal to that of the external medium, the 
amplitude C in (4) is zero, so that 1 = 0. For higher pressures, within the 
limits of the approximation, the jet variations are given by (7). 

The condition « = U, onx = 0 for 0 < y < 1 has not yet been satisfied. 


x 


On «= 0, u = U+ ¥ (A, €,/a)A(é,y). Multiplying by yJ,(é,y) and inte- 


r=1 
grating between y = 0 and y= 1, it follows from well-known 
properties of the Bessel function that 


A, = —20(0—U,)/€2,(é,). 
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FORMULA 


mn the jet boundary y ~ 1 and the gradient is given by 


dy hy _ a l ; E x 
2k? S — sin 27, (9) 
dx u —— . x 


here k* is written for 1—(U,/U). The roots &, of Jj(€) = 0 tend to (r—})z 


sr tends to infinity, and the series in (9) converges uniformly except in 


e neighbourhood of the points x 0, 2x, 4a, ete. Integrating term-by- 
rm between the limits x 0 and x, 
<i. .€2 
y—1= 4k? - sin® =<. (10) 
= 


his series, being uniformly convergent for all 2, represents the boundary 
ith gradient given by (9). 
It has not proved possible so far to establish exactly the whereabouts 
f the first and subsequent minima of the displacement given by (10). 
However, by means of a computation and the knowledge that > >? = } 
=1 


T 


see below) it was proved that the first minimum occurs for 


A, = «D,(2—1), 


here 1-15 << « < 1*3. A computation over the first forty terms of the series 
10), using tables of zeros of J,(€) compiled by Willson and Peirce (9), appears 
place the minimum between «x = 1-20 and 1-24. We may therefore write as 
ternative to (8) i72 
A, 1-22D,(M? 1). 


The velocity on the boundary is derived from Bernoulli’s equation 


1+ (y—1)M?/2 = NY-UY, 


/ 





where .V po/ Il as before. With y 1-41 and x 1-22, this gives 
A, = 2:695D,(N°291— 1-205). (11) 
Noattempt has been made to compute the lengths of the second and later 
wes. In what follows, by wave-length is to be understood the length ,. 
It is to be noted that, for a nozzle of given diameter, the wave-length 
s independent of the Mach number on exit; it is a function only of the 
ratio of the reservoir pressure to the pressure in the external medium. 
Now A, = 0 when p,/Il 1-9. This is true for the case VM, 1, M, being 
the Mach number at the orifice, since the pressure of the issuing jet is then 
equal to that of the external medium. When the velocity of efflux is 


) 


supersonic po/ Tl 1-9, but unless the ratio is great enough to ensure an 
exit pressure p, Il. a shock-wave will occur at the orifice and the jet 


will contract upon entry into the outer medium. The equation (11) gives 
value of A, corresponding to this intermediate stage, since the linearized 
equation does not make the distinction between compression and expansion 


192.10 N 
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at a corner which is derived from a study of the full non-linear 
equation of motion. Corresponding to the pressure ratio for which p, = I] 
there is no perturbation of the supersonic stream; for greater values, 
equation (11) will express the wave-length of the small perturbations, 
The value of K has been calculated for a range of values of V, by 
equating the formulae (1) and (11). The results are given in the following 
table, + together with the Mach number M corresponding to the pressure II, 
giving the velocity on the boundary and hence the exit velocity for which, 
with the given value of V, there is no perturbation; the theory is strictly 
appropriate only to small oscillations about these particular Mach numbers, 


TABLE (y 1-41) 
1 2 } 6 5 10 12 14 16 1d 
K rr6 | Irs 1700 | 0702 | 0°86 | 0°82 “78 | o-75 | 0°73 | OF! 0°69 
I I*go I°d9 1°05 IS! I*4I I°34 1°25 1°24 1°20 1°10 [33 
M I 1°04 1°56 1°53 2°01 2°16 2°27 2°37 2°46 2°54 2°6¢ 
The mean value of AK in the range V 1-9 to 10 is 0-94; in Emden’s 


experiments over the same range the mean was 0-94 (nozzles A, B, C, D). 





The mean value in the range 2 to 20 is 0-83; for Emden’s nozzles G and H, 
used in this range, the corresponding mean value was 0-78. 

Hartmann and Lazarus found the value K 1-12 for a series of tests | 
in which N was kept constant at 4-5. For values of NV between 3-5 and 
7-5 (approximately), the value 1-06 was found for the second and subse- 
quent periods of the jets; for the first period the value 0-98 was obtained, 
but a term 0:30D was added to the right-hand side of (1). 

It is seen that the corrected values of AK agree much more closely wit! | 
experiment than those originally given by Prandtl, and that the agreement 
with Emden’s results is extremely good. 

4. Plane jets 

For plane jets ¢ = Ux+A’' cosaf’ysin p’x | 
is a typical solution. Taking the half-width of the orifice as unit of lengtl 
the condition for constant velocity on the boundary gives the possibl 


values of : B; (2r—1)r/2a (r 1.2 ) 


dy; 9 Sseecde 


Thus the general solution of the linearized equation is 


@ = Ux+ > A,cos{(2r—1)zy/2}sin{(2r—1)za/20} 
r=1 
and the condition to be satisfied on x 0 is 
U,—O0 = ¥ {(2r—1)nA‘/2aheoss(2r—1)ry/2%. 


how 
r=] 


The comparison given in the table is an indirect one, but the experimental data are 1! 


T 


a form which is not amenable to a direct comparison with equation (11). 
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4 NOTE 





ON PRANDTL’S FORMULA 


linear By the usual methods for Fourier series it follows that 











1 | sx (0—U) 
alues A, ( I : =< 
‘. a” (2r—1)* 
V hy? On the boundary 
owing dy v the. — l iis (2r—1)zx 
ire I] da u wT it OF l 2a ; 
ricth [his expression is the Fourier expansion of the ‘square-topped’ periodic 
jherg | {unetion; it has the sum k*a{ 2/2], (8), where [x/2a]is +1 or —1 according 
s the integer next below x/2« is even or odd. For a/2« integral the sum 
s zero. The jet thus consists of sections of trapezoidal form. 
Since the tangent to the boundary in an exact solution of the jet 
roblem is continuously turning (at any rate until affected by shock- 
aves), this result does not represent physical reality ; it indicates, however, 
| that the jet is periodic with a wave-length / = 4a, or, introducing the 
: : | vidth, D, of the orifice, A= 2D, (2—1). (12) 


nd H,| The wave-length is seen to be greater than that of a circular jet with the 


same initial conditions. The ratio is 2:«, and K in the table must be 


F tests| multiplied in this ratio in order to obtain the corresponding values K’ for 
5 and} a plane jet. The mean value in the range NV 1-9 to 10 is A’ = 1-55. 


| 


subse- | Particular values of A’ are given in the table. 
ained,| There is at present little material available with which to determine the 
range in which (12) is a sufficiently good approximation. Good agreement 


with | results from a measurement from a photograph given by Prandtl (5, 
5 = a 


pment | Fig. 28) for a plane jet with WV, 1-85 and a small under-pressure corre- 
sponding to V 5:8 (ef. table); the measured wave-length is 3-1D, the 
calculated 3-0D. In Fig. 27 of the same paper, for M, 1-85and N = 7-74, 


the corresponding figures are 3-8) and 3-4D. It is possible to compare the 
— vave-length of the jet calculated theoretically in an earlier paper (3), 
“H vhere the reservoir pressure was 7-1 times the pressure in the external 
ssi} medium and the exit velocity was 1-5 times the local velocity of sound. 

rhe calculated wave-length was 4:9), which is to be compared with the 


value 3-3D given by (12). The error here is seen to be very large. 


5. The accuracy of the linearized equation in predicting wave- 
lengths 
Whether or not the linearized partial differential equation gives a 
sufficiently accurate description of the actual motion depends upon the 
are in Jlative magnitude of the terms omitted from the full second-order equa- 
tion, The most important approximation consists in replacing u?/a?—1 by 
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M?—1, and this is equivalent to neglecting 2AM in comparison with J. 
AM being the variation in M within the jet. The expression AM /M is the 
same for two-dimensional and axially-symmetrical jets, according to the 
linearized theory, and the restriction 2AM M o(1) does not permit any 
great amount of excess-pressure on exit. This leads one to ask why the| 
theory gives a reasonably good approximation to the wave-length of a 
supersonic jet over such a wide range of pressure ratios. 

The answer to this question is to be found upon examination of the 
conditions near y 1. The amplitude of the deviations from the mean 
velocity diminishes from the axis of a jet towards the boundary; thus, 
although there may be considerable changes of state in the neighbourhood 


of the axis, there will not be a corresponding deviation from the Mach 





number M on y = 1. In this way the linearized equation may give an 
adequate approximation to the wave-length of the boundary perturbations | 
under conditions where it could not be successfully applied to predict the| 
state within the jet, in particular near to the axis. 

A measure of the accuracy of the linearized equation near the boundary 
may be made by examining the magnitude of the expansion of the 
boundary, assuming the linearized equation to hold. 

For the plane jet, the maximum boundary displacement is obtained at 


once from the considerations of section 4: 


9) ~2].2 
Ymax I aG A ry 
For the circular jet, , 
= Anv2p2 Y ¢£-2 
Ymax I 4ov*h lant * 
r=) ' 


The sum of the series ¥ €>* may be evaluated by consideration of thi 
r=] 


integral p 2~*J,(z)/Jo(z) dz taken round a semicircle of infinitely large radius 
on re(z) = 0 as diameter, there being an indentation at the origin and th 
curved boundary not passing through a pole of the integrand; the sum 0! | 
the series is }. The maximum displacement is therefore less than one-hali 
that of the plane jet with the same reservoir conditions. It follows that the 
boundary approximations will break down at lower excess-pressure 0! 
efflux for the plane jet. 


6. Conclusion ' 

It has been shown that, for small fluctuations about the reservoil 
pressure for which a supersonic jet enters a medium with zero excess 
pressure, the wave-length may be found in terms of the mean velocity 
in the jet. The mean velocity is the velocity on the bounding stream- 
lines, and not the velocity on exit as assumed by Prandtl (4). With this 
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ith W.} wrrection, a much better agreement exists between theoretical and ex- 

‘is the! perimental values for the wave-length of cylindrical jets with sonic exit- 

to the! relocity; there is good agreement for a certain plane jet photographed by 

1it any| Prandtl. The theory may be expected to have a wider range of application 

hy the} to jets with axial symmetry than to plane jets, as the deviation of the 

h of a} boundary for the former is less than one-half the deviation for the latter, 
} calculated from the linearized equation of motion. 





of the In conclusion it should be observed that Emden’s formula arose from 
» mean} lis pioneer work in this field, and that it was the result of fitting a curve to 
: thus.} the results of his experiments. In future experimental work it would be 
urhood} more valuable to test the range of validity of the formula (11) derived from 
> Mach e linearized equation of motion. The considerations of the previous 
rive an| section suggest that the linearized equation will give an approximation to 


vations} the wave-length of a supersonic jet for pressure ratios outside the range 


‘ict the} for which it may be applied to give the conditions inside the jet. 
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ON THE LINEARIZED POTENTIAL THEORY OF 
UNSTEADY SUPERSONIC MOTION 


By K. STEWARTSON (The Royal Fort, The University, Bristol) 





[Received 11 August 1949] 
SUMMARY 

The Laplace transformation is applied to several problems of unsteady supersonic 
flow, assuming that the approximations of the linearized potential theory are valid, 
Simple expressions for the lift distributions on the semi-infinite wing, the swept-back 
wing, and the slender body of revolution at incidence are obtained. A method for 
obtaining a power series for the lift on a ‘ Delta’ wing is also given. 
1. Introduction 
THE theory of unsteady supersonic motion is by no means so well developed 
as that of steady motion, chiefly on account of the greater complexity of 
the governing differential equation. The first writers on the subject were 
C. Possio (1) and 8. v. Borbely (2), who determined independently the | 
pressure distribution on a laminar wing of infinite aspect ratio, executing | 





simple harmonic motion about a line parallel to the leading edge. The 
methods used implied the independence of the flows above and below the 
wing, and consequently it was not possible to extend them to mor 
complicated problems, such as the rectangular wing, where the two flows 
interfere. Several writers, following them, have refined their work, but 
none have overcome this fundamental difficulty. The most successful o! 
these were J. E. Garrick and 8S. J. Rubinow (3), who solved the problem | 
of the symmetrically swept-back wing with leading edges outside the Macl 
cone from the tip (see Fig. 1). Recently, however, J. C. Gunn (4) has 
applied the Laplace transformation to the solution of various stead) 
problems, and in this paper we shall be concerned with its application 
when the flow is unsteady. 


2. The formulation of the problem 

For simplicity it will be supposed that the wing is stationary (apart 
from an oscillatory motion) and that the fluid is moving over it witl 
undisturbed velocity V at infinity. The direction of the stream is taket 
to be the z-axis and a convenient point on the leading edge of the wing 
(such as the tip) to be the origin of coordinates. The wing is laminar and 
inclined to the main stream, but for a first-order theory it may be assumed 
to lie in the plane y = 0. Finally the a-axis is chosen so that Oxyz form 
a right-handed triad. Thé motion is assumed to be irrotational with 
velocity potential V’(z-+-¢), and deviations in velocity components, pressure, 
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nd density to be small so that squares and products may be neglected. 


‘om the theory of sound the equation of motion of small disturbances is 


veg — OF (1) 


ot?’ 
ere c is the velocity of sound, which, in our order of the approximation is 


nstant. Transferring to axes moving with velocity —V in the z direction, 


that z is replaced by z+-V’t, (1) becomes 
| 9 22 
Cc“? »< =4 ra O°D 
eVv2d 2 eee, (2) 
cor czct cz* 


hich is the general equation of motion. 
Ina similar way the excess pressure at any point due to the disturbance 
be shown to be 


od .¢ 
Ap = —p¥(2 +1 *}. (3) 


ct CZ, 


Hence the determination of 4 on the wing is sufficient to determine the 


excess pressure on it and so solve the problem. 


It will be sufficient to assume that the wing is executing simple harmonic 
otion with respect to time, so that ¢ enters only as an exponential 


fact r exp(twt) Let 


hb == y(x, y, z)exp twt; (4) 

Fa . gr Ob — - 

he c27V ds V2 2iw | + (iw )* a5. (5) 
C02* Cz 


This equation may be simplified still further by writing 





lw M | 
u; ae. .ZIex 6 
x (2, y, z)exp - WP | (6) 
where M is the Mach number of the undisturbed stream. Then 
x i x a x - xy] (7) 
5 Fi c y~ \¢ _ 
ere M2 -—-1 and w? C24, (5) 


In practice the wing will have a mean non-zero angle of incidence to the 
in stream and a superimposed oscillation about it, but since it is being 
sumed that the linearized theory is applicable, this motion may be 
ivided into two parts as follows: 
a) a steady motion of the fluid past the wing at constant incidence; 
b) a simple harmonic motion of the wing about the position of zero 
incidence. 
The complete motion is given by the sum of (a) and (b). The solutions 


a) for all the wing shapes in which we are interested are already well 
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known and attention will therefore be restricted to (b). The boundary } This 


condition may then be determined as follows. Let the angle of incidence 
of the wing be Ae’’, and the downward displacement of the leading edge 
from its mean position (taken to be in the plane y = 0) be Be’. Then the 
vertical component of the fluid velocity on the wing is 


[VA+iw( B+ Az) lexpiot, (9) 
which is equal to ¢d/ey. Hence on the wing (part of the plane y = 0) 
uf ; ; 
v8.3 VA —iw(B+ Az). (10) 
cy 


On the rest of this plane, the pressure must be continuous. 


3. The wing of infinite aspect ratio 


This wing occupies the region y = 0, z > 0, so that the boundary 
condition is 
/ 
Cw , ° 
= —VA—iw(B+-Az) (y U,. 2 2 ©). (11) 
cy 
Define us e-P= dz, (12) 


0 
then & satisfies 
oT oT 7 —T 
Cus C“ys e « tw tw\"i of ‘ 
—+— pa? +-2M p— 4 fp = pr, (13) 
6 cy” Cc c ‘ 
since % and éys/éz vanish when z = 0. The boundary condition becomes 
Cus 
cy 
where g(p) is a known function determined from (11). Thus we require a 
solution of (13) so that +0 as y—>oo and such that (14) is satisfied. 
Since ys is independent of a, this is 


g(p), (14 


r. (p} s 
ub _§g P’ exp( wy), for y > 0, (15) 
L 
T g(p) , , 
w= 4 exp(-+py), fory < 0, (16) 
B 
so that (B), a _ 9p )seny (17) 


B 
and hence 


sg f | } v(z—2’ v(z—z2’)' ’ 
Wh « % _ seny | (= exp(—'=26 z ) Jp w( : ) de (18) 
; 3 CY] y=0 eC ee 


applying the inverse transformation formulae and the Faltung theoren. 
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dary | This result agrees with those of Possio, Borbely, etc., and has been 


lence | obtained in a very simple manner. 


edge | 4. The general expression for the velocity potential 





n the In this section we shall find a general expression for the velocity potential 
rising from a given excess pressure on a laminar wing at incidence. This 
( ill enable us to determine an integral equation for the excess pressure 
))  } required for the wing to have the motion prescribed in section 2. To begin 
ith a solution of (13) is required with a singularity at r = 0, tending to 
(10 ero as r-> &, and proportional to sin@é = y/r, where r? = (x—2’)?+ y? 
nd 2’ is constant. 
These conditions are all satisfied by 
ib “ {Ko (rp)} (19) 
PO AoE), 
dary ey 


| where Ay is a Bessel function of the second kind with imaginary argument. 
Consider os T( re’ \bo dz’ , (20) 


(19) | Where f is an integrable function of x’. Then as y — 0, 


; P yf dx’ = = 
us lim —“______. == of f(% , p)sgn y. (21) 
0 (x Cre 


« 


13 . , : 
\ Hence (vs), , is equal to 7 f(2,z)sgny, where f is the inverse Laplace trans- 
form of f. Thus the excess pressure on the wing is simply related to f. 
mes ‘ ; ‘ 


Moreover bv considering ¢/¢y, an integral expression for the normal 
(14) | velocity on the wing in terms of f may easily be obtained. The expression 


20) above is therefore the general formulation of the problem of incidence 


ire ofa laminar wing. To determine the expression for the velocity potential 
sfied in terms of a, y, z we proceed as follows: 
; ey [ 96 
; hi exp(—pr cosh @)cosh 6 dé. (22) 
(15 7 . 
Let tan Acai” 
(16 peo®+iwM’ 
ind Oar w. (23) 
(17 Then 
db, 2 1 na as (p | toe cosh 9 ‘sinh d| x 
r P | cCa* 


9 
ca” 


exp) | (p +12 eosh 9 ‘sinh 9 (24) 


ren. the remainder of the integral being imaginary. 
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Let z yr cosh 3, & yr sinh i, then 
: * de joM\ iwoM\ 
wy — re = (p+ a iAE |exp (ae } rE (25) 
r I ¢ Cae | ca®* 
Combining this result with (20) we have 
fb - | | Y wo dz (: 2")oos( AE— ae C f(x’,2’) \ ) 
t J 4 t Mz — P 
+ (Ms, Jsinfre— |fte’.e IF (26) 
x“C x°C 


The integration is to be taken over that region bounded by z’ <:z, 
€* = (z—2z’)? 


off the wing, and this places certain limitations on f there. If the wing is 


v°72 > 0. In all problems the excess pressure must vanisli 


raked forward it implies that f vanishes off the wing, but if it is raked 


backwards this does not necessarily follow. The method of overcoming 





this difficulty has been found in the steady case (5), but in the unstead) 
problem it has so far proved insuperable. 

From expression (26) we can deduce the integral equation for f(x, 
quite straightforwardly, but it will not be given here as it is only of 
academic interest. 


5. The oscillating slender body of revolution 

This section is concerned with the lift on a slender body of revolution 
pointed at the forward end, which is oscillating in yaw. The boundary 
condition is, as usual, that the normal component of the fluid velocity o1 
the body be equal to the velocity of the body. It will be assumed that the 
body is oscillating about a mean position, of zero yaw, having its axis 
lying along the z-axis. If the mean position is not of zero yaw, it may b 
dealt with in a similar way to that for the wing. Moreover, since we are 
assuming that the linearized theory is applicable, the boundary condition } 
may be written down as though the body were at its mean position. 

Hence wr:ting x = rcos6, y = rsin6@ we find that there are two relevant 
contributing factors to the normal velocity of the fluid at the surface ol 
the body, 

(i) the variation in incidence of the body, Be‘, 

(ii) the displacement of the nose, Ae”. 
Combining these terms we see that if the nose of the body is the origin 
and the body occupies the positive z-axis, then 

Cus 


> —[VA+iw(B+Az)|sin@ as r>0,2> 0. (27 
er 
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\ solution of (13) is required which tends to zero as r > #, and which is 
roportional to sin@. Hence 
: 


ob = ph(p)sin 0K, (rap), (28) 


here h is the Laplace transform of A(z), a function to be determined. 


‘or small 7 


_ Bi ws ] 
um h(p), uy ~ h(z), (29) 
x7 xr 
cus — 
das r— 0, - > h - sin @. (30) 
Cc? _* fi 
Hence h(z) tim WV A+iw( B+ Az)}S(z), (31) 
SIN 0 r—-0 07 


here S(z) is the cross-section of the body at z: the solution is now 


plete. 


6. The symmetrically swept-back wing 

This section is concerned with the symmetrically swept-back wing, 
issuming that the Mach lines from the tip lie on the wing. The boundary 
condition is Gb/cy -e there, and from the solution for this problem, 


t for the boundary condition @s/ey = ez may be deduced by integration. 











*, 
o 
A 
5 
Fic. 1. The symmetrically swept-back wing. (a) Leading 
edges outside the Mach cone. 
OA OA’ isa OA, OA’ are the leading edges of the wing. 
OB OB’ is OB, OB are the Mach lines through O. 


\pplying the Laplace transformation with respect to z, and the Fourier 
transformation with respect to x to the differential equation satisfied by uy, 
e find that 
ey’ — 
c* [ p?(V*—c?) +-2tw Vp +- (tw)? +78? |¥, (32) 


cy* 
re s is the parameter of the Fourier transformation and 


t'(s,p) = w(x, 2). (33) 
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If tan-'o = 7 is the angle of sweep back, then the boundary condition is 


or 2Qeo 
. when y 0. (34) 


cy s°-+o7p* 
Hence, using the method of section 3, the solution is such that when 
y UO 


Zeop ] 


py 5 (35) 


s*+o07p* \ { px” 1 (27m c)pM 1 (iw /c) 


9 


. a"; 
and when y = 0— the sign of the right-hand side is positive. 

There are two ways of tackling the interpretation of (35); we could 
transform 2eop/(s*+-o7p*) and the remainder separately and combine the 
answers by means of the Faltung theorem. In this way we arrive at the 
answer given by Garrick and Rubinow, which suffers from the disadvan- 
tage that the computation of the integrals is laborious. Alternatively we 
can treat (35) as one function and determine its inverse transform as such. 
It will be found that for this specific problem the result is considerably 
simpler and, even taking into account the linear term in the boundary 
condition, the numerical work is lessened a good deal. 


Since the Laplace transform of 








, . / yp s° La M2’ > 
| cos so-l(2—2' Jy) 2 om 1) JeXP ecole F (36) 
I? al Nate? 2) 2c 
0 
is equal to (35), then when y 04 
L , F i eee 
cy € de | dz’ cosso 1(z 2 \Jo) 2! | ae | 
vs 2ro , A | Ni vic? © c?/} 
L 0 
iwMz’ . | s 
<exp| ——; ist). (37 
x“C 
Consider | ds Jy) 2’ (5+5 loxp fer —*@—z'y) ; (38) 
JOO a Nate ca] feP | 2 | 
and let me gax(z-— 2" Veg, (39) 


Then, using Sonine’s integral for J,(¢), (38) is equal to 











} iz y 
l - etdu fC pra wz? «22/28? 
- | ds exp} — St. = =— - (40) 
2m . , 4 4uaitc? 4ua" 
b iz —< 
b+i« 
> = 99 \ 9 19 
X du v2X = wz ~ 
oo | exp u(1 —— | — = b 0) (41) 
Z2aN7 J NU \ - he 4uaxr 
b—i« r 
° = 9, wt? 2Vv2 
~ *¢ L ™~ ti oe 
¥ COS(w/ aC). ( x2X 2 ) (42) 


\(2’2—a2X?2) 


if X72? z'2. and is zero otherwise. 
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Nn is If now we also write X x+(z—z’)/o, then 
5 
(34) Cus * dz cos| (w/arc )a Z x2 X? )] iw Mz’ 
mt : ~ = “eXp - = } 
Cz \ (2 v7X2 ) x“ 
hen | . 
' dz’ cos Ww x7C) (2’2 v2X2 ) iwM2’ . 
- ee lexp end, 
vo } \ (z'2?—a?.X2 ) Ni 
) 
here z,,, 2, are the zeros of z’*—a?X? if they lie in [0,z] and otherwise 
tenes re 0, z respectively. Similarly for z,, z,. Consider now 2’?—a?X2.. This 
the wnishes when 
the 2! (2 —a2a)/o(ato). (44) 
Val > nr ° > > 
Suppose first that a 0. Then if z/¢ x, both of these values of z are 
he egative, so that the integral vanishes. The point here is that throughout 
uch, , f ; ° . : . m.: . 
bie | the interval [0,2] the integrand vanishes in virtue of (42). This result is 
aD | : . ‘ . . e ° ° 
| to be expected as this region is in front of the wing and so the fluid is 
ar’ a 
disturbed. 
If z ro both zeros are positive, so that the lower limit of integration 
; > > ro 
sthe smaller zero of z’” YD, Gall 
If z vv the larger zero is z and is therefore the upper limit of 
(3t ; : . . 
tegration. If z vv the larger zero is z and z is therefore the upper 
t of integration 
Secondly take a 0. If z vv the larger zero is z and z is still 
upper limit of integration. If z xx both zeros are > z and the 
gral vanishes. Again this is an expected result, since it means that 
’ 
form of the leading edge to the right of Oz can only affect that part 
27 : , , oe . . : 
\o ofthe wing to the left of Oz, which is inside the Mach cone from O. Consider 
first integral of (43) (= J, say) in the region where it does not vanish. 
(38 ox), : 
} x >COS C}, zy ~COST—TSMNT, (45) 
“— Go” 
(39) 5 , 
| 
iw Mz, cost . 
c(.M? cos*7— 1) | w2z,sin & 
| 7o COS 
(40 7,/(M? cos?z—1) we, (LM? cos*z— 1) 
{ig Mz, sinz cos € : 
expla i) de (46) 
(4) x ( M? cos?r ] ) 
M? cos?r— 1 ‘ M2. 
ere cos €, (x cos7+2s1n 7T)- sin 7, (47) 
1? XZ, COST x 


i+ and is zero otherwise. 


ifthe upper limit of integration in (46) is 7, then the expression reduces 
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to the value of és cz for the two-dimensional wing with swept leading edge, 


lis COUIa, OF Course, De evs atea airec IsSIng a similar metnoc ( 
Phi uld, of cou I uluated directly using | thod to 


section 3, but it may be deduced from (46) quite easily. In fact 


iw Mz, cost 
cos exp 1) 
J 





I c(.M? cos?r Ww, ‘ 
. ; 8) 
70 ,/(M? cos?7— 1) "\e( M2 cos? — 1 ), a 
iw Mz,’ Dz 
As7—>0 exp _ eat) (49) 
70 x*c J \are 


in agreement with (18). 
A similar expression may be derived for the second integral if we write 


2 cosT-+-xsin7, and change the sign of 7 in (47). Finally writing 


wM A 
Q 7 , (90) 
vc (.M? cos*r 1) 
and M62 M? cos?r— 1, (5] 


we find that 
1 
aO Cds r 


7 


exp(?Qz,) | cos(@Qz, sin €)exp(?Qz, )sin 7 cos € dé 
M cost cz . : = 
0 


of 
¢ 


) 


-exp(tiQz,) | cos(@Qz, sin €)exp(—7?Qz,)sinz cos € d&, (52 


0 
and this expression enables us to evaluate the pressure at any point ¢ 
the wing. If, however, we want to find the total lift and the moment « 
the lift about the a-axis, it is simpler to use (35) instead of carrying out the 
integration of (52). 
Let the chord of the wing be /; then since 


(FP) <6 | da | « Pah(x, 2) dz, (53 
avis 0 
the left-hand side is 
: css ot ; ; ; (54 
ap*,{p>(I 2 ¢?) + 2iwV p+t-(iw)?} 
Hence 
] =/b 1 
F i ‘sony [ eee iMwz' . ai 
dz | us dx pastes at | (1 Z ee be exp Pe: a dz’. (50 
R » XO A x°C XC 
0 2/¢ 0 


This result is comparatively simple and shows that the lift and moment 
coefficients for a swept-back wing are independent of the angle of sweep 
back, which is the same result as for the corresponding steady problem. 


If w is small, (52) and (55) may be expanded in a power series in o, 1 
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lirect contrast to the source solution of Garrick and Rubinow. Alterna- 





ge, 
d to, tively, we may expand (35) and interpret the power series, but there is 
ery little advantage in this. 
The solution satisfying the boundary condition é/éy = « on the wing 
is thus been determined, and from it we may now deduce that satisfying 
(48 cy ez on the wing. 
Let the solution of é/ey = ez on-the wing be ¢,, with Fourier—Laplace 
49 ansform ‘f}, and that of é%/cy =e, be %. Then , ~ the boundary 
\ tt. . 
ue 4eapo*/(s?+-c7p?)? on the wing. But WY ~ (2ec)/(s?+ po?) on the 
ng; hence 
+ J » c ° ~~ 
vrite pt; o?— (V/oc), (56) 
;, co 
) 
| ] 
Cw CW ~=— 
(50 id therefore by Co re. (7) 
C2 Co 
\o r . . J . . 
. The semi-infinite wing 
In this section we shall consider the simplest semi-infinite wing, i.e. that 
ounded by the positive « and z axes. The boundary condition is that 
cy =e on the wing. The solution satisfying e/éy = ez on the wing 
nay then be deduced by integrating with respect to z. Taking the Laplace 
(59 transform with respect to z, the boundary condition becomes @ys/éy = €/p 
hen y = 0,2 > Oand s = Owhen y = 0,2 < 0. Now it has been shown 
nt ( by Lamb (5) that the solution of this problem is such that when y 0, 
nt « V, 
tu 7 e. £4 ‘ 
son yw eT dr. (O38) 
pe pram, 
lhe form of the solution is not very suitable for interpretation but we 
(53 nsee at once that the first term is the ordinary two-dimensional solution 
lid outside the Mach lines from the origin O. On reducing the second 
m by a simple change of variable, we find that 
D4 - 
pibsgn y l zt = a 
exp(—.xpsec?6) dé, (59) 
€ [L 77 ph 
0 
a hich is the Laplace transform of 
wnt # son y cus Ww » F : i w o o@ iMwz 
— Jol ; dO Jol — . (22?— ax? sec4#)]} |exp| ——— 
wee] © <a x“ 7X , °C XC 
U0 
ylem (60) 


w, ll here sec*6, zn. 
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The total contribution to the lift, for a fixed z, from the second term 
of the left-hand side is 


an unexpectedly simple result. For small values of w, the lift from the 
region inside the Mach lines from O is 


. 07 2; 
Pep] (1. Ziel (4 9), in 


20? —3V a2 a 


where / is the chord of the wing. 

In all the problems so far considered, the solutions have been of a 
relatively simple nature, but in the next problem, that of the ‘Delta’ wing, 
not only has it been proved impossible to find any kind of exact solution, 


but even an expansion in a power series in w presents great difficulties. 


8. The ‘Deita’ wing 
This wing has symmetrically swept-back leading edges, and lies entirely 
within the Mach cone from the tip. The method used to solve the previous 


problems breaks down when applied to this wing, because it fails to provide 
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by mul 
eglect 


2 Dz ° . Dz I aya 
. : i= (cos? +iM sin or jexn( —" nso )|: 61) | 
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sould € 
the pro 
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vy taki 
fthe F 
steady 


—- 


a solution in the steady.case. Various methods have been successfully ' 


developed for dealing with this problem, but none offer any chance of: 
solution when the wing is oscillating, and accordingly, hope of finding a 
exact solution was abandoned. The only remaining possibility seemed t 
be to obtain a solution in ascending powers of w and in this section we shal 
carry this out. A general method is developed and used to obtain th 
first three terms of the series. More terms may be worked out if desired 
the only difficulty being of an arithmetical nature, but we can see, eve! 
from these three terms, that the chance of there being a simple solutio! 
is small. The governing equation will be (7) and the boundary conditio 
will be that on the wing 


CX . »9 
x : 5 A, 2". (60 
cy 0 
Hence we know that on the wing |x! < pyz, pp < 1, €x/cy is given, betweet 
2 poz and |x z, it is unknown, and outside |x z it vanishes. O! 


the other hand, x vanishes everywhere on y = 0, except on the wing 
where it is unknown. 
The Fourier-Laplace transformation of the governing differentia 


equation is 


(X) = (p?+A?+s?)X, taking a a (64 


where X(p,8) = x(x, 2). 
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term) Given the value of y on the wing, the value of €X/éy may be found merely 


multiplying X by (p?+A?+-s?)!, and hence the correction due to the 
) neglect of A may easily be found, on inverting the transformation. This 


as ld enable us to obtain a solution correct to order A, and by repeating 
e process, as high an order of accuracy as desired may be obtained. The 
ee | ie problem, therefore, is the solution of the simpler question obtained 
taking A = 0. The method given bélow was deduced from consideration 

(42) \ ofthe Fourier—Laplace transformation of the solution of the corresponding 


steady problem. 


j 
ot a 
‘ N 








ving >X 
tion 
es 
ws \ A B’ J B 4 
vious 2 
vide Fic. 2. The symmetrically swept-back wing. (b) Leading 
ful edges inside the Mach cone. 
OB is OB’ is2 OB, OB’ are the leading edges of the wing 
ol OA is OA’ is OA, OA’ are the Mach lines through O. 
ig a 
ed t If sgn y F(p, 8) X is the solution we require, then it follows that f(x, z), 
shall | the inverse transformation of F', must vanish if |2| > pyz, while the inverse 
n the transformation of F',(p?+-s?) must vanish if |x| > z, and must be equal to 
siret 63) if |x Da2. 
evel Two integrals will now be evaluated, which will together satisfy these 
ution conditions. 
litic Firstly, consider 
(63 I; 4 p” dpe! ~ Js P (6 > 0), (69) 
amt J (p'-r8"po)" 
wee" | where f(x) is a polynomial of degree n, and py < 1. Then the integrand 
‘i m is a regular function of s except at pys Lip, where there are branch 
wit | points. Hence the integrand may be made regular everywhere by inserting 
. ts in the s-plane from pys ip to © arg(—ip) and from p»s Lop 
entia fargip. The integration with respect to s wil! first be recast into a 
re convenient form. 
(64 


. * ds ¢ isrf(g 
mside1 I] Ji P) (66) 


) 2 2 .2\n +3? 
am J (p +8" po)" * 
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and take x > 0. Then the contour of integration may be completed as} } follow 
shown in Fig. 3. The integral round the infinite semicircle vanishes and 
the integral round the circle C, of radius e, is of the form 

















; 
e p A,é, (67) 
0 
’ : mm : _ pepo 
where the A’s are undetermined constants. The integral along the cut is 
equal to 
IAL . ix \ 
. exp(—! me ) aa { 
7 | Po | Po 
. . 7 ss (6S) 
7 Py "aia, aaa | iatide 
l-+e 
writing pgs = —ipa. The sum of these two integrals is equal to (66) and 
s-plane | 
D 
A 0 A 
R ae 'R 
(va lv 
\ ne / 
NN / 
\ ¥ 
\ ce / 
j Re, 
\ , , 
\ / If no 
‘ ee 
Me. - E (pts 
-iR polynor 
ur pro! 
B 
, , ] of the : 
“iG. 3. Contours for J,. . 
oh a , slightly 
OA is real axis; OB is imaginary axis; R, KR the original path of integration ; . 
Dis —ip/p 9; D’ is ip ; DE, DE’ are cuts in the s-plane; Cis a circle centre D case \ 
radius ¢€ [referred to in text]. = s 
(ul) S 
hence must be finite. Therefore the infinite part of (68) must cancel out 
with the infinite part of (67), as « > 0. Since the finite part of (67) is zero 
it follows that 
bh+ig f j Jar\ . da) 
Papen. pee) 
I, — | = i | fo os (69) | Where / 
re er +2 2 v+3 : 
27" pyt p" (a*— 1)"* at s 
2 J : 
Insertin 
where.F is to be read as ‘the finite part of’. | As befo 
Inverting the order of integration and noting that complet 
tha 
: exp! p ~ =] ( yan\ +1 uw 
. y 2- : ' a r | 
; - dp Po / Po 
2a ree | ; 
Ps. 0 otherwise 
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dl as} follows that 


, “\ n+] 
and Po axs(—")(2—™) 
] \ a Po fe 














; _ 5 ' (71) 
mp,(n+1)! | (a*—1)” 
67) ' 
p) > wand J, = Vif zpy < x. Hence J, satisfies our condition on (x),,~». 
ut 1s 
fF | S plane 
68 L A 
R iR 
\ / 
t / 
) and D, 
“ 
~ 
Fic. 4. Contours for J. 
pA ist is: OB imaginary axis R, R the original path of integration ; D is 
E is E’ is ip; DE, D’E’ are cuts in the s-plane. 
If now we can prove that on multiplying the integrand of (65) by 
s*), the new integral vanishes when z < x and is a homogeneous 
olynomial of degree 7 if pyz rv, then we have solved the first part of 
ur problem, i.e. found a general method when A = 0. The determination 
{ the new integral is similar to that for (65), and we shall consider a 
, slightly more general case, which will prove useful, when dealing with the 
use A 0. 
Secondly. consider 
lo 
Ze! . . e—'s© dg f(s/p) ” 
J ¥ prer* dp a ae (72) 
ua J (p?+s*?!4(p?+-8*p5)" 
(89 ere / is an integer > 0. The integrand is a regular function of s except 
it 8 ip, and p,s -ip. Hence it may be made regular everywhere by 
serting cuts from s ip tO pos ip, and from s = ip to pys = tp. 
\s betore, consider the integration with respect to s first, take 2 > 0, and 
mplete the contour as shown in Fig. 4. Then, as with J,, we may show 
” 1 F XP f(x) d 
el ) » ¢ te ) ( 
P ip fla) da 


, ®\_i —" ‘ (73 
al . p* wii p. (1 + x7)! “(1+ x7 pzynri ) 
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; 


Let « = —indw, an elliptic function to modulus ,/(1—p§), and invert the und sin 


order of integration. Then 


uy 


J (—)! ~ { {z—(xnd w/p,)}"*2f(—ind u)(dn w)2"+2! dy 
, 7(n-t 2/)! (1 — pe)" + ; (en u)?n +2(sn u)?! 2 
i 


(74)? Hence 
where 





- ; ee ,' . R (iii 
(a) if zpyp > w, Uy kK, the complete elliptic integral of the first kind, Dn ) 
i > Take 
(6) ifz > x > 2po, dnu, = poz/z, 
(c) ifz<z, J, 0. 
If zp, > x, J, is a homogeneous polynomial of degree n+ 2/ in z and a, and! 
hence, by a suitable choice of / and f, it satisfies the conditions on (@y/éy),_).| as the 
The only problem left is the determination of the appropriate polynomial 
and this will be done below for n = 0, 1, 2. The boundary condition is F 
that (@x/éy) = > A,2" on the wing, and we shall consider the terms of this . 
. 0 ' and si 
series separately. 
(i) Suppose ex /ey 1 on the wing. 
The solution of this problem is well known (see for example 6) but we 
shall give it here for the sake of completeness. Take / == 0, iy | awe 
‘ : . m since « 
in the equation for J,. Then 
(81) th 
K | 
Ae f ; AE we get 
i. F du se*u — =. (73 ‘ 
= mp2 } and B 
0 | 
where £ is the complete elliptic integral of the second kind. Therefore 
’ 
A 7po/ E. and th 
Moreover, 
Zpo/.r 
PY?” ae ae “an 
I, ——,F pa a 3 V(2"pp—*") (76 
7 po z (t#— 1) 7 Pp 
i 
and hence (y),-9. = - Br Po x?) on the wing. 
and hi 
(ii) Suppose €x/ey = 2 on the wing. | The 
Take 1 = 0,” 2, f = A in the equation for J,. Then | Table 
K 
j A is [ dn u(z dn u—a)sn2u du (77)| 
e J 4d)? 
t a(1—p?) entu fe 
0 44 
™ 4; 
——;-- —-[p2(K — F)+(1—p§)E]. (75) f (45 
37 p5(1— p5) 4: 
| 


| 
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rt the} and since J 2 we can determine A. Now 
1 
. 'v (py) 2—2Xa)? Az,/(p§ z2— x?) ” 
" i, a Ff da * - -_—= ; 5 P (79) 
~ AU 27 ps J (a*— 1)? 37 pp 
1 
(1— pj )24/(p5 222?) 
- 1 err-vire* = » 
(74)# Hence (x)y=-0 20K E Bd =" (80) 
po\ 4)-+- (1 — po) 
iS iil) Suppose cx cy 2 on the wing. 
cind, oh a _a. 
Take 7 = 0, n = 3, f = A+ Ba*p?. Then 
K aha 
l ~ [ du sn®u dn?u(z? dn?u+2?)(A — Boz nd*u) ; 
J; _. . , (81) 
277(1— pz)? } en®u 
v, and 0 
y s the term involving xz has no finite part. Now 
oma kK » dn2u d ! 
° ” sn- anu au 7 . e e . ee 2 - 
100 1s a ; pate (AK — E)p§(1+- ps) — E(1— p)(2—p?)} (82) 
en°’u L5p4 
of tl 
nd similar expressions may be derived for 
hk K 
* sn7u du — [ sn®?wdntu du 
, - and F - 
ut J en®u : en®u 
( 0 
Since J, = 22 we must choose A and B so that on the right-hand side of 


81) the coefficient of z* is equal to unity and that of x vanishes. Thus 
_ | we get two simultaneous equations, from which we may determine A 
(io) i 1D - 
ind 6. For example, 


67 pal (K — E)p2(6p2 —4)-+ (8—9p2)(1— p?) E] 


ore As as aan 2 : —3 5 > (83) 
5p3( AK — F)?— 2p? E( K — E)(4—p?)—(4—7p2) E7(1 — ps) 
l there is a similar formula for B. Moreover 
a I, F * (A— Ba)(zppo aa)? dx (84) 
i Gzpéa : (a?—1)8 


VE Po) 2292/44 + B)—22(4+4B)], (85) 


l hence, substituting for A and B, we get x: 
the coefficients of z* and x? in (85) are tabulated as functions of p2 in 
T ible | below. 


TABLE I 


0-1: 0-2: 0-3: 0-4: 0-6; O-S8: 1-0. 
14+] 
{= 2 2-41; $-O1; 5°28; 6°36; 8-16; 9-68 ; ll. 
(78 ' 1B i 
0: 0-35: 0-53: 0-65; 0-73; 0-86; 0-94; —]. 
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In this way we may continue to determine solutions corresponding t; 
the boundary condition éy/éy = 2" on the wing, where n may be as larg 
as we please. It can be seen, however, that the expressions with whicl 
we have to deal are rapidly becoming more complicated, and there does 
not appear to be any way of overcoming this difficulty. The solutions 
obtained here are sufficient to solve éy/éy = expipz, A = 0 on the wing 
as far as terms of order w*. To complete the solution we must determin: 
the error, to the same order occasioned by the neglect of A? in the funda- 
mental differential equation. The correct form of this equation is 


~(X) = (p?+s?+A?)X. (86) 
oy , 
We have deduced that, if @A/e¢y = 1 on the wing and A = 0, 
. 7 pe ] 
(X y 0+ = aan 1° (87) 


E (p?+s8?p?) 
Due to this X the correct value of @X/éy in the plane y = 0 is 

7p (p?+-8*+-A*)! P| (p?+8*)* a 
E (p?+s*p?): E 


ina ol eC ae + O(A‘)|. 
(p?+-8*pj)! ° 2(p?+-8?)!(p?+-s"p5) 


(SS) 


The first term of this expression has been interpreted as unity on the wing, 
while, making use of (74), the second leads to 
€ > K A 

p> d2 ~ { (ednu—x)? du 5 eo uo. me 
——— FP - — ~ > | pa z°(K — EB) +-2?( A p5— E)| 
4(1—pi)E en?u 4H (1— p?) 
0 (Ss! 
on the wing. Off the wing, of course, éy/éy does not assume such simpl 


values, but this region does not interest us. So long as y is zero wher 
|x| > poz the solution is correct. 

Let the appropriate polynomial in J, required to cancel this term (89 
be f(a) C+Dp?.?. Again we can solve for C and D, and obtain the, 
corresponding (y),,9, as in (85). A table of their values is given in Tabk 
II below: 

TABLE I] 


pi 0; O-l; 0-2; 0-3; 0-4; 0-6; O-8: 1-0. 
—(40+-D) 0: 0-66: 0-84: 0-92: 0-96: 1-00: L-Ol: R 

Y p56 N 
—=(4D+C) 2:36; 1:91; 168; 1:52; 1-41; 1:23; 41:10; 1. 

YO Po ' 


Hence combining all these results we find a solution of the problem as) 
a power series in w? of which the first three terms have been calculated. 
It can be seen that it is most unlikely that there is a simple general 
solution owing to the complexity of the first few terms. However, the , 
method given above is clearly capable of supplying as many of the terms | 
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ing to of the series as may be necessary, and superficially at least, appears to be 


; laroe much simpler than any other known method. To determine the integral 


which of (x),-9+ over the wing, from z = 0 to z = 1, the simplest method is to 
e does put s = 0 in J,, omit the integration with respect to s, and multiply the 
27. We find 


utions answer by 


l b+iz 
> Wing rf 1(O) ‘ q eb f(O)in+2 
: 2 hoa dadz . dz dp - J . (90) 
rmine a 271 pnts (n+-2)! 
‘unda wing 0 b—i« 
The principal outstanding problem of wing vibration is that of the 
yeneralized semi-infinite wing, i.e. one which is raked and swept. If the 
(SO 


wing is raked forward, we can obtain a solution using Gunn’s method. 
This solution is complicated, but may be expanded as a power series in w. 
There is nothing new in this and so it will not be included here. If the 
(87) wing is raked backwards, however, the solution is not unique, but it seems 
ikely that the one occurring in practice satisfies the Joukowski hypothesis. 

Taking this to be axiomatic, the procedure is as follows. First rotate the 

1) ], ixes about Oy, until the raked edge becomes the Oz-axis. Determine the 
Fourier-Laplace transformation of the solution of é¢/ey = 1, A = 0, on 
the wing, and from this the solution of é¢/éy = expipz, A = 0 on the 


wing 


g may be deduced. If then we correct for the neglect of A? as in 
section 8, a solution may be obtained to as high an order of accuracy as 
may be desired. This is not, however, wholly satisfactory, and it is hoped 


E) +] 


iat a better method will be evolved. 
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THE FLOW OF TWO ADJACENT PLANE SUPERSONIC 
JETS PAST FLAT-PLATE WINGS 


By M. HOLT 
(Department of Mathematics, The University, Sheffield) 


| Received 20 December 1949] 


SUMMARY 


The linearized method of conical fields (1, see also 2) is applied to determine th 


used. I 


conside 





startin: 
The 

of the 

’ the fiel 


pressure distribution on certain flat-plate wings at incidence to a stream composed 


of two adjacent plane uniform supersonic jets with different velocities and densities, 
and separated by a plane vortex sheet. The leading edges of the wing are straight, 
intersect on the jet boundary, and lie outside the Mach semi-cones of this inter- 
section (see Fig. 1); the edges may be either swept back or swept forward. 

The pressure coefficients on the upper wing surface inside the Mach semi-cone it 
the port and starboard jets are determined. It is shown that there are five typical 
cases for the pressure distribution, and the pressure coefficients on the upper surfa 
have been computed with values of the parameters chosen to represent these fi 
cases. The results are set out in Tables I and II and shown graphically in Figs. 8 to 12 

Two examples of physical application may be cited. Firstly, the solution gives t 
pressure distribution on the part of a control surface of a supersonic aircraft whit 
is subject to interference from a propelling jet. Secondly, the solution may be us 
to give a first approximation to the pressure distribution on a wing in a non-unifor 
parallel stream, by replacing this stream by a series of discrete adjacent unifor 
streams and applying the method of superposition described in (2). 


1. Introduction 
WE consider the flow past a flat-plate wing set at a small incidence to 
stream composed of two adjacent uniform plane jets with different w 


disturbed densities and supersonic velocities. The jets are semi-infinite i!) 


extent and are separated by a vortex sheet, which is plane upstream of thi 
disturbance. The wing is normal to the common boundary of the jets, wit 
a leading edge composed of two straight lines intersecting on the je! 
boundary. The wing (Fig. 1) may be either swept back or swept forward 
but each branch of the leading edge is taken to lie outside the corresponding 
Mach semi-cone. (The corresponding problem for a wing with leading edgt 
inside a Mach semi-cone will form the subject of a further paper.) Onl 
first-order effects of the disturbance are considered, and the distortion 
the vortex sheet downstream of the disturbance is therefore ignored. 
Since the problem is formulated so that no fundamental length i 
involved it may be treated by the linearized conefield method of Buseman! 
(1): the form of this method developed by Goldstein and Ward (2) will be 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1959)] 
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sed. In fact the present problem is a generalization of one of the problems 

NIC. considered by these two authors and the results they give form the 
starting-point of this paper. 

The notation is extended from that of (2). In the jet on the port side 

f the wing the notation used is identical with that employed throughout 

the field in (2). In the jet on the starboard side quantities relating to the 
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surt 
se f YZ 
8 to I? Fic. 1. AOA’ is the wing leading edge, Oz is the trace 
ves t f the vortex sheet, and OB, OB’ are the traces of the 
wl Mach semi-cones. 


undisturbed stream have different values, and variables affected by such 


nifor | differences are distinguished by means of a prime. 

The origin of a rectangular Cartesian coordinate system is taken at the 
apex of the wing, the z-axis is in the direction of the main stream flow and 
the z-axis is in the transverse direction. The angle of incidence of the wing 

. Ys to the plane y 0 is denoted by x. 

: : In the port jet (in x < 0) we denote the values of the undisturbed 
wr velocity, density, and Mach angle by W, p, and p respectively, and the 

- components of the additional velocity by u, v, w. Throughout this region 

de : we employ non-dimensional variables 2,, y,, r, @ defined by 

~wan a/(ztan x) = 2, r cos 0, y/(ztanp) = yy rsin @. 

mding  Qutside the Mach semi-cone of the port jet we also employ the variable o 

gedg! defined by seco = r. In the starboard jet the corresponding variables 
Onl employed are denoted by W’, p’, y’, u’, v’, w’, x}, y}, 7’, 0, o’. The ratios 

‘on tanu)/(tanyu’) and (pW)/(p’'W’) are denoted by AK and c respectively. 
L. Without loss of generality we assume throughout that K > 1. 

gth is We determine the flow in the region y > 0 only. The flow in the region 


mani / < 0 can then be found from the fact that u, w, wu’, w’ are odd functions 


will be of y, and v, v’ are even functions of y. 
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2. Boundary conditions 


On a linearized theory we may take the boundary conditions 


v Wx in the port jet, (1) 
y’ —W’ax in the starboard jet, (2) 
i 


A 








uy 





£6 plane 
Fic. 2. 
to be satisfied on the patt of the plane y = 0 that is the projection of th 
wing. 
At corresponding points on the vortex sheet the disturbed pressure an 
the disturbed normal velocity are continuous so that 
u u | 
cw w’ | 


at corresponding points on the plane x = 0. 


3. Conditions outside the Mach semi-cones 
Port jet 


Figure 2 shows the r,@ plane for the port jet when the leading edge lies 


inz > 0. If f is the angle of sweepback the r coordinate of the leading } 


edge is given by 7, = cotyucotf. From a result established in (2) it follows 
that the disturbance outside the Mach semi-cone is confined to the region 
1, in which v -Wa. 

Figure 3 relates to a wing which is swept forward on the port side. 
The r coordinate of the leading edge is now negative and equal to —7r, 
where 7, = seco, = cotycotB and £ is the angle swept forward. The 


compression shoek starting from the leading edge is reflected at the vortex 


sheet as an expansion wave. In the o,@ plane the compression wave starts , 


from (7—o,, 7) and separates regions 3 and 4. The reflected wave starting 
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In region 4 there is no 


—Wa. 


-Wx is propagated along the 


7) separates regions 2 and 3. 
listurbance: in region 3 the disturbance is constant and given by v 
Hence a discontinuity in v of amount 
haracteristic separating region 4 from region 3; it follows that an equal 
ind opposite discontinuity in v is propagated along the reflected charac- 
and 2. 


teristic separating regions 3 Region 2 is therefore free from 











rd 
/ 
A / 
4 | og 
; a +——__+— > 
. R | 24 














disturbance. In region 1, v lx as before. Thus the conditions on the 
Mach semicircle are the same for both types of wing provided we take f 
as the angle swept forward when the leading edge lies in z < 0 and as the 
angle swept back when the leading edge lies in z > 0. In view of this 
result further consideration is restricted to swept back wings. 

Starboard ye i 

In the starboard jet two cases have to be considered, in each of which 
two regions outside the Mach semi-cone are disturbed, the first adjacent 
to the wing and the second adjacent to the vortex sheet: in the r’, @ plane 
the corresponding first and second regions are bounded by the lines 6 = 0 
and @ = 47, ares of the semicircle r’ 1, and the tangents to the semicircle 
from the points (75,0) and (A, 47) respectively, where r, is the value of r’ 
on the leading edge. In the first case the two regions are independent, in 
the second case they overlap. 

Figure 4 illustrates the first case. In the o’,@ plane the disturbance in 


region 6 is constant along characteristics o’ —@ constant, so that 


nu’. v’, w’, at P’ u’'.v'.w’ at Q’. (4) 


In region 5 the disturbance is constant everywhere with 


Y" W'a. (5) 
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In the second case, to which Fig. 5 relates, conditions in regions 5 and 6} when ¢ 
are the same as those in the first case. In the o’,@ plane, on crossing into | eadins 


region 7 (where regions 5 and 6 overlap), the disturbance on a charae- | 


teristic c’—@ = constant changes discontinuously so that 


v'(R’), = v'(R’),—W'c, ete. 
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Hence, using results given in (2), 
w'(Q') = u'(P’)—W'acoto; | 
v'(Q’) = vo (P’)— Wa 
w'(Q') = w'(P’)+ W'atan p’ cosec a}, | 


where 
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shen Q’ is on the boundary of region 7. Also, o4 is the value of o’ on the 
eading edge, so that , ’ , 
ad 1s SEC Oy cot p’ cot f’, 
here 8’ is the angle swept back by the leading edge in the starboard jet. 
At A, 7’ K: hence at H, 

'e) sin-1(1/K). (7) 


A 4B 

















4, Conditions inside the Mach semi-cones 
Inside the Mach semi-cones we employ complex variables ¢ and ¢’ corre- 
sponding to the port and starboard jets respectively. In (2) the variable ¢ 
is chosen so that the trace of the wing surface should be represented by 
part of the real axis. Here, in order that the boundary conditions on the 
vortex sheet may be expressed in the simplest terms, we choose ft and t’ so 
that points on 6 = 3m correspond to points on the real axes of ¢ and ¢’ 
respectively. 
Port jet. We define ¢ by the relation 
1/t = sin 6/r+-{2(1—r?)! cos 6}/r. (8) 
The transformation from the r,@ plane to the ¢ plane is shown in Fig. 6. 
The trace of the vortex sheet C’'A is transformed into the segment (0, 1) on 
the upper real axis. The trace of the Mach semi-cone A B transforms into 
the part of the upper side of the real axis on which t > 1. The trace of the 
wing C'B transforms into the upper imaginary axis. The coordinate of M 
on the real axis is ¢, r,/(r?—1)!. The components of the additional 
Velocity wu, v, w are the real parts of analytic functions of t, G,(t), G,(¢), 
G,(t) respectively. It can be shown, using results given in (2), that 
dG, /dt: dG,/dt : dG,/dt (1—?é*)?: —2:attanp, (9) 
where the principal value of the square root is taken. 
Starboard jet. The corresponding variable ¢’ is defined by 
1 /t’ K sin 6/r’+- K{i(1—r’?)! cos 6} 7". (10) 


he transformation to the ¢’ plane (in the first case) is shown in Fig. 7. 
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' 
The coordinate of G is | = r,/{K(r2—1)!\}. The components of the addj- | ()) v; ¥ 
r4 2/ | 2 j 
tional velocity are given by 


u’ = RG;(t’), etc., 
where 
dG jdt’: dG,/dt’ :dG;/dt’ (1—A*t’?)}: —7: —it' tan. (11) 


, . hae , ? 
We again take the principal value of the square root. 








0 YK i iz c } In | 
C D oH G aes 


In case 

















and 
Referring to Fig. 4, we now find the relation between the coordinate of 
P in the ¢ plane and the coordinate of the corresponding point Q’ in the ; 
t’ plane. Firstly, We 
: (P) = r(P) = r*, say, basen 
and at the corresponding point P’ these 
rr’) = Cat 

Hence 6(Q’) la7—oa'(P’) 47—sec-!r*’ = cosec-!r*’, j 6.7 
Th 

e ‘ gi fg a Ae ‘ 
o that U(Q') = r*'/k defin 
and '(@") = a P). (12) | has| 
plan 
5. Boundary conditions in the complex planes ditic 


he solutions in the complex planes are found for three sets of boundary | resp 


conditions, corresponding to the following three cases: velo 
I. ¢, 1, tf; > 1 (the first case discussed in section 3) the 

If. 4; > 1, 1/K < t < 1 (the second case discussed in section 3) S 
Ill. ¢, 1, 4, = 1/K, corresponding to a wing with leading edge normal Tun 
to the stream direction. the 


The boundary conditions are 


(a) uvew 0,1 <¢ < t, onthe upper real t-axis (in cases I and II only) 
v Wa, ¢ > t, on the upper real t-axis | 
v Wa, on the right side of the upper) in all three cases. 


imaginary t-axis 








FLOW OF TWO ADJACENT PLANE SUPERSONIC JETS 
e addi- uw’, v',w’ = 0, 1 <t' < t, on the lower real t’-axis 
py’ W's, t’ >t on the lower real t’-axis |. 
: jie ; ; in case I, 
v W’a, on the right side of the — 
imaginary t’-axis 
(I p’ W’a, t’ 1 on the lower real t’-axis 


) in cases II 
v’ W’a, on the right side of the idl 
} and III. 


imaginary ¢’-axis 


4 In ease I, when ft r+Or, t’ r—QOr 


eases II and III, when ¢ r+O2, t’ r—Ot 


“ul u’ . ‘ 
tor Oar < G 
Cw w a 
nd u = w+W'ol K42—1)99) . | 
ate ol ; ‘ ar 8 | for ty << Z. 
| cw w i xf, tan po 
n th 


We now find the solutions in the complex planes in case I. It will then 
he seen that the solutions in case II have the same form and that from 


these the solutions in case III can be deduced. 


| 6. The solutions inside the Mach semi-cones 

The solutions in the complex planes are known when once a function 

defining the velocity components on the disturbed part of the vortex sheet 

(12) | hasbeen found. Accordingly, our first step is to find the solutions in each 
plane which satisfy the boundary conditions 5(a) and 5(b) and the con- 

ditions w P(r), w’ = cP(r), O0O<r<1, when t=r4+0, t' = r—O 

dary | respectively. We then express the condition of continuity of normal 
velocity across the vortex sheet in terms of these solutions and so obtain 


integral equation for the unknown function P. 


Solutions in the t plane. The function P(r) is represented by a step 


‘mal | function, so we firstly find the basic functions dG,/dt, dG/dt satisfying 
the boundary conditions 
uw a) () , S 1) 
us 0 S , , rt On 
nly) ( ' r<t, | when ¢ = r+0 
v Wr 7 t, 
_ v Wax, on the right side of the upper imaginary axis. 


Then d@, dt is purely imaginary on both the real and imaginary axes of f, 
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has simple poles at ¢ = s and ¢t = t,, and dG,/dt = O(1/f*) at infinity. We 


i The so 
therefore take 


dG, i. x 8 that 
dt #—s? ‘ @—# 
Ai t—s' Bi ttt 
and G(t) ‘log prog ("| 
2s t+-s 2t, t,—t 
From the condition w l,t=rt0i, O<r<s we find A = —2/7, and 
From (9) 
dG, Ai Bi The 
dt t(t?—s*)tanu  t(#—é*)tanp L<l 
Ai {?—s? Bi ¢? 
and G,(t) — — log = — — log|—_—-}. nd 
“ 28° tan pu é? 2tj tan pe ti—t* 7 
The condition v = —Wx, t = r+0i, r > t, then gives > the 
B = (—2¢]7Wotan p)/z. 
| 7, Th 
ent OWyit2 +t: 
iia dG, z sr = ni ase (13 In : 
dt a(t°—8*) 77(ty -t*) across 
From (13) we build up the solutions satisfying the general boundar 
conditions 
: Using 
. P(r) vile ti integ 
w= 0, ice <t when ¢=r+0i ' 
v Wr, YS 1] R 
v -Wx, on the right side of the upper imaginary axis. 
1 } 
‘ G. ) [ 2P'(s)\ds 2Whatk? ite 
We find dG, ’ | 2 e de arene (4 
dt a {-—s-* a(t7—t*) 
0 ; 
] 1 T 
G 1—#?)! ( 2sP1(s)ds | 2Wal? 2)! a ( 
and — : , | an e ss U oni. , (15 | 
dt 7mttan p , t-—s* at(t;—t*) ; 
0 
where P!(s) is written for dP/ds. ibis 
} 


Solutions in the t' plane. The corresponding general boundary conditions 
in the ¢’ plane are 


w’ cPi(r), = r= J | 
w' = 0, l<r<t} when t)=r—Oi 
f=; (J 


y’ W'x, on the right side of the lower imaginary axis. 
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Vy. We 1 


he solutions satisfying these conditions are found by the same method 


3s that used in the ? plane. They are 


1 
dG; on) | nit a , 2W'cdy's tan yp (16) 
dt 7. t-—s* a(t.'—t *) 
0 
1 
, dG,  c(l—K*t’?)! [ 2sP'(s)ds  2W'at?(1—K?t'?)! (17) 
9 nt 7 r 1? > — 1/4!9 19 ° ‘ 
~ <8 /77, dt mt tan p t-—<s* mt (to-—t °) 


These two sets of solutions are also valid in cases II and III. For when 
l, (16) gives 
RG; (t’) cP(r), t’ r—0Oi, O<r< i 
nd 
RG,(t') = cP(r)+ Wat, tan p, t’ = r—Mi, i <¢r« i, 


0 that the boundary conditions 5(c) are satisfied. 


7, The integral equation for P(r) 
In all three cases the condition that the normal velocity is continuous 


across the vortex sheet gives 


is t=r+ 
ndat R dG, /dt RdG,/dt'’ when | 


r—Oi } 


Using (15) and (17) it can be shown that (18) is equivalent to the following 


for O<r<l. (18) 


integral equation for P(r), 


1 
* 2s P}(s) ds 2Wat?(1—?*)! tan pw 


R} {(1—#?)! +-e(1 — K#?)}} —e 
(292 PP 


. 
0 


2W’ ot?(1— Kt?) tan w 


(14 t2—#2 
0, when t=r7r4+0:. (O<r<l). (19) 
To solve (19) we firstly write it in the form 
RH(t) = R[ F()—G()] = 0 (O<t< 1) (20) 
it ls understood that ¢ varies along the upper edge of the real axis), with 
aie F(t) = {(1—#*)!+¢(1— K¢?)!} f(), 
tions 
4 2s P1(s) ds 
#292 ’ 


A,(1—#)! | A,(1—K#)! 


1 t2—{2 


A, 2Wat* tan p, A, = 2W’'at? tan p. 
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| 

Evidently (19) is satisfied when ¢ —r+0i (0 < r < 1); when ¢ takes §, The 
real values outside the interval (—1, 1) f(/) is purely real. Equation (19) is The 
therefore satisfied along the whole real axis and as t >, H(t) O(1)t).8 she Mav 
[It follows that H(t) is a rational function of t, purely imaginary when { is 
real. Therefore the only singularities of H(t) are poles, which must be 
poles of F(t) or G(t). a 

In cases I and II the only poles of dG,/dt are at t tt, and the only We de 
poles of dG/dt are at t +t,. It follows from (14) and (16) that F(t) has } magi 
no poles: hence F(#) has no singularities of any kind. G(¢) has simple poles} Port 
at t Lt,, t Lt}. We therefore take 

H(t) = it{A,/(#—t®)+ A, /(t2—#)} 

and choose A, and A, so that the residues of H(t) are equal to the residues! where 
of —G(t) att Lt,, ¢ Lt, respectively. Thus We pu 


A; = A,(#j—1)!/t,, A, A (Kt? 1)/t; 


‘ 
and 
] 


(t : A,{(1—é?)!+ 7t(8 — 1)! /t,}/ (8 —#)- 
= q=pyrcea- Kp): 1{(1—#?)! + at({—1)#/4,}/ (GP) 
A {(1—Kt?)!+ it( Kt? —1)8/t)/(t2—#)]. (21 
In case III we can apply the same argument to show that F(t) has n 
poles. The only singularities in G(t) are now integrable ones at f 
and ¢ +1/K. The function H(t) is therefore regular everywhere, an 
since H(t) >0ast>a, H(t) 0 for all ¢, and 
f(t) = {A,/(1—#)! ++ K24,/(1— K22)4/((1— #2)! +e(1— Ke). (22 
The function P(r) is determined from the relation 
P(r) I f(t)/7, t=r+03, Ox<r<l. 
Equation (21) includes (22) as a special case. Thus the solution of the 
integral equation is defined in all three cases by (21). 
From (13) and (21) 
dG i 
dt ari (1—??)'+-c(1— K#?) 
L Af(1— Kt?) 4+ it(K%t2—1)8/t,\/(t2—0?)]. (23 


9 


) 


, [A.i{—e(1 — K4?)! + ae(tf]—1)4/t,}/(G@—)- 


It can be shown that /(t’) is also defined by (21) (in the lower half of the 


;u here 
und sc 
form 


(" 


p 


Ste 


wher 


i The 


wher 


Thu 


t’ plane) if we replace ¢ by ¢t’ and change the sign of ¢ wherever it appears 


explicitly. Hence 


dG i 


dt’ wf{(1—t’?)§+-c(1 


Kaya Ant 144-0) (G1) 


1 Af(1—t’2)!+ cit’ (K%2—1)!/t}/(t2—t’?)]. (24 


9 
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| 
t takes) § The pressure distribution on the wing surface 
‘(19)ish The pressure due to the disturbance on the upper wing surface inside 


O(1/t).) se Mach semi-cones is given by 


picts p pWw_ in the port jet, 
just b ; ioe ee , 
) p p' W'w’ in the starboard jet. 
~ Ve determine w and w’ by integrating (23) and (24) along the upper 
(t) ha iaginarv axis of ¢ and along the lower imaginary axis of t’ respectively. 
sf . 
e poles! Port jet. If we take the lower limit of integration at t = ix, 
i 
w(t)—wp R | (dG,/dt) dt, 
sidues | where Wp Wx tan px cosec o,. 
Ve put t 1a (1 az} LQ) 
xy k’ cos d/(1—k? cos?¢)!, (25) 
here h (K2—1)!/K, : (1—k?)}, 
ind so obtain the expression for the disturbance pressure in non-dimensional 
form d 
» , . 
” p 1 2h dd ro 
( . 2a tan pe . 4 | 2 24\h. ~> 
1as 1 Lo W? (r7—1) a | {e+(1—k* cos*¢)*} 
‘i 0 
Wr, hk’ cosd+er : ' 
e, a woo, rk cosd+er) \) oO <$<4n). (26) 
\W(r,+cosd) r7 (k’?4 k?r?)cos*} PI 
2 
(9 Starboard jet. With the lower limit of integration at t’ = —ix 
7 t 
w' (t')—w, R | (dG;/dt’) dt’, 
' here wy Wak’ tan p cosec a4. 
of tl 


} The non-dimensional pressure is now found by putting 


t' ik 'Z (1 ary 1-(), 
vhere ty cos ws. 
Thus 
(2 ( p 
7 ) Soll - 
of the 7. =? 
| ws k'W'rg Qk | ds . 
ypea»s és an pu — = 9 9 i) “> 
pea . me W(r—1) a J {c+(1—k* cos*s)'} 
, r,W'[r5(1—k? cos*b)! +e cos p]) 
\7 cos J r,(1 k? cos’)! Wr? cos*y) | 


(24 (0O<p<4n). (27 
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(To avoid a discontinuity across the vortex sheet the divisor }pW®?, rather 
than }p’W”, is used to define C,,.) 

It may be shown that the sign of dC,,/dx, cannot change more than once| 
in the interval (—1,0): the same is true for d¢ 'y" dx, in the interval (0,1). 





Further 

cele the sign of dC,,/dx, is identical with that of gully wpe 
ha j | W’—eW 

at x, = 0 , , oa , , "—eW 

. 1 the sign of dC,,,/dx, is identical with that of oll ' od 
xy 1} (mV cW, 

W here ] Ie’), (r?- ; ] ) r(k’4 cr1)(re+ l ), 


m = k’'r(1+-1,)(e-k'rs)/ry(r?—1). 


and! << 1,1l< m. 

Thus, if dC,,/dx, is negative at x, 1 it may change sign once in the 
interval (—1,0), otherwise dC, /dx, must have the same sign throughout 
this interval. If dU,,/dx, is positive at x, = —1, dC,,/dx, must also lk 
positive at x, = 1. Finally at 2, = 0, x, = 0, dC,,/dx, and dC,,/dzx', mus 
always have identical signs. It follows that the only possible types ¢ 
pressure distribution correspond to the five following cases: 


(a) Ata, -1, de ‘0 dz, < 0;2, 0, d¢ « dz, < 0;2;, 1, dé * dx, <0 
(b) Ata, —1,dC,/dx, < 0;x, = 0,dC,/dx, < 0; 2, = 1,dC,,/dx, > 0 
(c) Ata, 1, d¢ s dz, <.0;2, 0, dt ' dz, > 0;2, = 1, a * dx, > 90 
(d) Ata, -1, de 4 dx, > 032, 0, dé ” dz, > 0; 2; 1, dé ‘' dx, > 0 
(e) Ata, Se “ dx, < 0; 2, 0. dé a dx, . 0: 2 1. de ‘0 dx’, <0 \ 


To illustrate these cases (26) and (27) were computed with different } 
values of 7,, 73, W’/W to cover the five cases and with fixed values of the 
remaining variables. The values chosen were 


k k’ L/v2 Cc > X 0-1 pe lor 
(a) W’/W 0-5 
(b) W'/W 6] ,=2 
; (ec) W’/W = 25, 1, 4 15 20, 
(c) W’/W 5 1, 2 
(d) W’/W 15 


The results of the computations are given in Tables I and II and corre- 
sponding graphs of the variation of C,, with x, and of C,, with a are shown ? 
in Figs. 8-12. 
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TABLE | 


Pre SSUTE COE fficve nts in the port jet 





, A=alg=s Roane G8 
| Wo 5 wu 1s wou 2°5 
( ( ( 
i ld i — 
231 31 0°2065 
[5 25 0°1933 
I 279 0182 
Ig4 211 0 174 
192 348 01694 
194 391 0°1667 
+439 oO 1661 
| 0° 1669 
S65 0°1687 
615 O'I7II 





TABLE I] 


Pressure coefficients in the starboard jet 
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It may be shown that in the limiting case when the two jets become 
identical, i.e. when p = p’, Vi W’, w = p’, (26) and (27) each reduce to 
he equation for the pressure coefficient on a similar wing in a single stream 
given in (2)). On the other hand, it should be noted that the present 
results do not include the case of a single semi-infinite jet discussed in (2) 

, (Le. Wor W’ = 0), since the whole theory is based on the assumption that 
oth jets are supersonic. 

The expressions on the right of (26) and (27) may be rewritten in terms 
lt inverse trigonometrical functions and of incomplete elliptic integrals of 
the first and third kinds. No obvious advantage results from this, since 
the new expressions are long and the elliptic integrals of the third kind 
which arise belong to a class not hitherto tabulated. 
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List of symbols 
Throughout the field 
x,y,2 rectangular Cartesian coordinates with origin at the apex of the 
wing, z measured downstream and z in the transverse direction 
(to starboard). 








Port jet 
W,p,u values of the undisturbed velocity, density, and Mach angle | 
respectively. 
u,v,w components of the additional velocity. 
23, 4,,7,@ non-dimensional variables defined by 
x/(ztanp) = 2, = rcos86, 
y/(ztanp) = y, = rsiné. 
o cos—1(1/r). 
3 angle of sweepback. 
ry cot p cot Pp. 
t r/isin 6+7(1—r?)! eos 6}. 
p additional pressure. 
C', = p/tpW?, additional pressure coefficient. 
¢, defined by the equation 2, = —k’ cos ¢/(1—k? cos*¢)!. 
Starboard jet 
In general, quantities corresponding to those defined in the port jet are 
denoted by the same symbols with a prime added. Exceptions are 
To cot pn’ cot B’. 
t’ r’/K{sin 0+-i(1—r’?)! cos 6}. (N.B. @ has the same meaning in both | 
jets.) 
Cy = p'/2pW?. | 
us cos ly. 


Throughout the field 
K = (tanp)/(tanp’). 
c = (pW)/(p’W’). 
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) thank 4 METHOD FOR THE SMOOTHING 
OF NUMERICAL TABLES 


By A. T. DOODSON 








te? 
(Liverpool Observatory and Tidal Institute, Birkenhead) 
1 Super 
Pd Received 19 July 1949] 
SUMMARY 
The smoothing of numerical tables based on observational data, or obtained by 
ethods which cannot be readily rendered more accurate, is a problem which has 
ceived much attention. There are disadvantages in methods solely numerical or 
of the | ¥lely graphical. The method outlined in the paper is a combination of both. 
' ‘ \ subsidiary function is built up of differences corresponding to those obtained 
ection rom the table to be smoothed and is then compared with the tabulated function. 
The difference between the two is then smoothed graphically. The method is also 
pplied to double-entry tables 
angie 


1. Introduction 

NUMERICAL tables which depend upon observations are liable to show 
irregularities which are a source of inconvenience if allowed to remain. 
In many cases the table cannot be improved by the incorporation of more 
data or by improvements in the processes of reduction from the original 
data. A similar difficulty occurs with mathematical functions of such com- 
plexity that their computation is dependent upon mechanical or graphical 
processes. In either case it may be supposed that the resulting table is 





substantially accurate and that the irregularities in the entries are due to 
random or casual errors. The over-all accuracy cannot be improved by any 
smoothing processes, but such processes may be regarded as advantageous 
in that any smoothed entry is made to conform to the indications of its 
, neighbours, and in this sense it may be considered more accurate. The 
concept of smoothness is based on graphical ideas, and the method 

et are ‘described in this paper ultimately depends upon the use of graphs applied 
| tosmall quantities which show up the lack of smoothness in the original 

table. Its special characteristic is the use of a subsidiary function built up 


both | ‘rom differences. 


2. Smoothness tested by the aid of a subsidiary function 

Suppose that a set of numerical values N, representing fairly well a 
tunction y in terms of the ordinate x, is compared with a function », which 
‘known to be a fairly good approximation to y, so that 7 may be regarded 
‘smooth while y is not smooth. Then if the difference (N —7) is small 
nough to be plotted on such a large scale that errors of a unit in N may 


(Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1950)] 
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be corrected by a simple graph, we can produce a set of smoothed values 
of (V—7n) which may then be added to » to give a smoothed version of V, 
It is not necessary that the function 7 be obtained on theoretical grounds. 


[t may be suggested by the general appearance of the curve for JV, and | 


simple curves such as parabolas, hyperbolas. sine curves, or logarithmic 
i ni 


curves may be used, without any necessity for a theoretical justification, ! 


Whatever the subsidiary function may be, it should be easily calculated 
to a degree of accuracy greater than that required in the table. The sole 
test of its usefulness is whether the differences (.V—7) can be examined 
graphically. 


3. Use of subsidiary function built up from smoothed differences 

A subsidiary function may be simply obtained without recourse to 
mathematical functions by the use of the differences of V. Such a sub- 
sidiary function will more closely represent the characteristics of V than 


any function unrelated to V by any known laws. A smooth representation ‘ 


of NV requires that successive differences also are smooth, apart, of course 
from the accumulation of small errors normally experienced when differ- 
ences are taken, but these are not of much interest in the present problen 
which is to build up a smooth function, so that we may make all differences 
as smooth as we like for-this purpose. 

Our object is to get a good fit to the second differences of .V, the first 


differences of V, and to V. It will be supposed that the second differences | 


may be plotted and a smooth curve drawn ‘through’ the plotted points 


With the smooth values obtained from the curve, a system of smooth first | 


differences may be obtained by successive addition to the first value of A. 
in the table. These could then be used to give a smooth set of values of 
by addition to the first value in the table. 

[t is doubtful, however, whether that straightforward process would be 


useful, for a slight systematic change in the curve for the second differences 


might produce marked results after successive additions through the first 
differences. It is necessary to see that we get a good fit for each set of differ: 
ences and finally a good fit in the final product. The best procedure is as 
follows: 
1. Plot the second differences of V and draw a smooth curve through 
them. Read off the smooth values 7. 


— 


2. Add these successively to the values of AN,, so building up the sub- 
sidiary curve 7,. Plot the differences (AN—7y,) and from a curve 
through the points read off the smoothed values of (AN—»,). Add 
to the smooth system 7,, and so obtain a smooth system for the first 


differences. 
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THE SMOOTHING OF NUMERICAL TABLES 219 
values. 3. Proceed in exactly the same way to obtain », and thence the smoothed 
of V system for .V. Test the final product by successive differences and 
yunds make any small amendments that may seem necessary. 

», and, The modifications necessary if the subsidiary function is built up from 
thmie | \icher differences than the second are obvious, but the difficulties increase 
ation. yith the number of differences used, as successive differences tend to 


ulated | }eeome complex in form and if there are many errors in N then the higher 
e sole | differences may not appear to be susceptible of representation by simple 
mined raphs. 

When used with discretion this method is very powerful, and it .has 


seceeded when the original table of NV’ seemed faulty beyond hope of 


Nces ae —_— an . 
redress. A numerical example is given in Table 1. The curves used in the 
rse to age 
yurse of the work can easily be drawn by hand, as is usual in this method 
1 su . : ; ‘ . , 
+] f smoothing, for the method does not require skilled draughtsmanship. 
} ? 
tation : . . sf , 
4. Application to the extension of tables 
OUTS ° . . . . 
ee [t is sometimes desirable to extend a table, which is smooth to the nearest 
ifte! . . . ° owe 
DI init, in order to give a table which is smooth to one-tenth of a unit. The 
Hien . ° ° . 
nethod of procedure is just the same as described above save that the 
rence . . . 
curves need to be read to a greater nominal accuracy. It is desirable to 
use the original data in such a form that second and third differences are 
e firs ‘ ; nm. : _— 
| appreciable in magnitude. This ensures that the original table has not a 
rences : : ; 
spurious smoothness due to entries at too small an interval. It is not 
OINtTS , ‘ : 
bf necessary to give an example as the process is much the same as in the 
1 first ; 
.,y, example already given. 
of A.\ : 


;of.)| 5, Application to double-entry tables 


The processes described may be applied to double-entry tables by 


ld be} smoothing separately each row and each column, but it is by no means 
ences! easy to obtain smoothness in all directions. A powerful aid, however. is 
e first) by the use of a graphical process in the form of a network, which can be 
liffer applied to the original function or to its differences. 

eC 1S as 


Suppose that the data are plotted on a single sheet so that curves may 
be drawn to express the data in the columns. This gives a family of curves 
rough} and each curve can be used to correct errors in the column to which it 

pertains, but errors in the rows would require a fresh set of curves to be 
»sub- drawn. The double set of curves, however. can be obtained on one sheet 
curve by one plotting of data if the curves for the columns are ‘stepped’; that is, 

Add, after the first column has been plotted the data for the second column 
e first should commence one unit of spacing to the right (or left, if that is prefer- 
ible) of the initial entry for the first column. Then the data for the third 





, 
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; 
column should be plotted so that the first entry begins a further unit of 
spacing to the right, and so on. The commencing points can be connected 
by a curve which represents the data in the first row of the table, and| 
obviously two sets of curves can be drawn, one family representing the 
columns of the table, while the other family represents the rows of the}; 
table. These curves have to intersect one another on the vertical ordinates} —.. 
of the diagram, so that a very accurate system of curves can thus be drawn 
Such a network can be used for the original function or for any differences 
It has been much used for the smoothing of tidal data. 17087 
The differences used in building up the subsidiary function may} 
be entirely in the columns or entirely in the rows, or a combination of) 
these. The example given in Table 2 uses second differences which are) 59° 
the differences along the rows of a table of first differences down the! 42280 
columns. In this example, to save space, the graphical networks are not 
given as they may be easily reproduced by anyone studying the methoi.| 
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unit of 
See 
are TABLE | 
le, an Example of smoothing process 
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Table (111) was then smoothed by a network; after some preliminar 
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NUMERICAL TABLES 
The first column of Table (11) was then smoothed and the rows of the 
lowing table completed by adding successively the values in the corre- 


onding rows of the above table. 
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ON THE CHOICE OF STANDARD SOLUTIONS FOR A 
HOMOGENEOUS LINEAR DIFFERENTIAL EQUATION 
} OF THE SECOND ORDER 

By J. C. P. MILLER 


Scientific Computing Service, 23 Bedford Square, London, W.C.1) 


| [Received 31 August 1949] 
' SUMMARY 

In this paper principles are set up for the choice of a pair of standard solutions 
omogeneous linear differential equation of the second order, suitable both for 
ematical theory and for numerical tabulation. The equation is taken in the 
mal form and, in general, two cases arise: (i) with intervals of 2 called exponential 
n which solutions behave roughly like positive and negative exponentials 
are convex to the w-axis; (ii) with intervals of «2 called oscillatory regions, 
ch solutions cross and recross the w-axis as they oscillate, and are concave 


ywards the axis 


In an exponential region it is shown that solutions should be chosen so that their 


aries as rapidly as possible, and not so that the ratio tends to a constant as 
es indefinitely along the interval; this would indicate that a relatively small 
ion was becoming ‘swamped’ by a large one. In an oscillatory region it is 
vn that suitable solutions can be chosen and from them a modulus and phase 
| by expressing the solutions in the form F'cosy and F’sinx, in which the 
} functions F and y are normally without oscillations. Methods of choosing these 
ons F and x to be as ‘well-behaved’ as possible are discussed. Particular 
nsidered in more detail are (a) the choice of solutions suitable in two oscilla- 
gions separated by an exponential region, or by a region of oscillations of 
vy long wave-length: (6) the case of a finite oscillatory region. A geometrical 
tion, particularly useful in oscillatory regions, is included. 


1. Introductory remarks 
, AXy homogeneous linear differential equation of the second order may be 
transformed to the normal form by change of dependent or independent 
variable or both (see 1, p. 394). Thus the equation to be considered may, 
without loss ot generality, be taken in this form, namely 
) d?y 


iy? I(x)y = 0, (1.1) 


which (2) is assumed to be a real function. 

The purpose of this paper is to set up principles for the choice of 
standard solutions which shall be both useful in mathematical theory 
ind convenient for numerical tabulation. 

Two main cases arise, needing different treatment: (i) Exponential 
‘gions, or intervals of x throughout which /(x) is negative and solutions 
(Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1950)] 

Q 
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convex to the a-axis; (ii) Oscillatory regions, or intervals of x throughout 
which J(x) is positive and solutions are concave to the x-axis. 

In section 2 the general character of the solutions in these two types 
of interval is considered. In section 3 the choice of standard solutions js 
discussed, and in section 4 a geometrical interpretation, particularly useful 
in oscillatory regions, is described. In section 5 a few remarks are made 
on the difficult case of a finite oscillatory region. 


2. General character of the solutions 

2.1. The character of the solutions to this equation depends on the sign 
of I(a) in the interval of arguments concerned. Thus, in an interval where 
I(x) = —J(x) < 0, the solutions are all convex to the x-axis and behave 
roughly like increasing or decreasing exponentials, or like a combination 
of these. In particular, solutions may in general be found which, starting 
with the same value at a particular argument x, become very different in 
size as x increases or decreases. 

In fact, writingt v = y’/y and substituting in (1.1) gives 


ye J (x)—v’. (2.0) 


Thence, by successive approximation or by assuming a solution in the forn 
‘ e Lo) 


y. 3 By 
yew +P- { abs E <a ae eee (2.12 
; 2.F 2.4.2" 244.5°"— 
with J = P?, and assuming P large and its derivatives P’, P”,... of at most 
the same order as P (often they will be successively smaller), the solution 
Pp PrP’ 32PP’—3P" yp} J’ 4J3J°—5J" 
v~ T ee _ ——.e9 Tut 8 . - 
; Ce cas 2.1. 40° 4.3.7 
(2.13 
may be obtained, whence the approximation 
yr~d \a)exp| 4 | Ji(t) in), (2.14 


valid so long as the variation in J(x) is not too rapid in comparison wit! 
its magnitude. Any solution y is a sum of multiples of these two solutions 


ry 


. : : F \ 
which, starting equal at 2 = x», bear the ratio exp|2 | J*(t) dt | to on 


at) 


another at 2 x, A particular solution can cross the axis at most once 


the solutions of (2.14) do not cross it at all. 


+ Here and later, accents denote differentiation. 
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2.2. On the other hand, if J(x) (x) > 0, the relations (2.13) and 
114) become 
, Q’ 200" —3Q"2 —-_ tf. .41f—sr 
1~ +iQ -§— wo tiJt—-_F 4 —— see 
2) 2.4.03 4] 4.8.13 
(2.21) 
teu 999 
| y ~~ ] 412 EXP) rt [*(t) dt (2.22) 
ding to the approximate general real solution 
y ~~ CF '(x)cos) I*(t) dt T ¢ (2.23) 
In this case the solutions are always concave towards the axis and will 
ntinually cross it so long as | /3(t) dt increases or decreases in- 
finitely.t Sturm’s separation theorem (see 1, p. 223) shows that any 





) vo solutions have interlacing zeros, and from (2.23) it is seen that they 
tain the same relative amplitudes throughout. 
2.3. In subsequent discussions an interval in x for which J(x) > 0 will 
» called an oscillatory region for the equation (1.1), whilst an interval 
j vith J(v) < 0 will be called an exponential region. The solutions will like- 
| wise be called oscillatory or exponential type solutions. These terms will be 
haken to apply to the complete intervals considered, even though these 
nay end at zeros of /(x), near which the approximations (2.14) and (2.22) 
nay break down. Modifications necessary near such zeros of J(x), or when 
t issmall, are considered only when they affect the trend of the argument 
see section 4.5). 
|3. Basis of choice of standard solutions 
) 3.1. For a particular differential equation (1.1) the range of x will be 
) livided into a succession of exponential and oscillatory regions by the 
zeros of I(x), or by other points at which J(a) changes sign, if discontinuous. 
| Even minima of | J() may cause significant separation into distinct regions. 
) For mathematical purposes, and particularly for computational purposes 
and for the preparation of tables for general use, it is necessary to choose 
| pair of standard solutions from which others can readily be constructed. 
) The pair of solutions chosen need not be the same in the various regions, 


} Dut it is usually desirable that changes should not be too frequent. 


For tabulation over large ranges of the independent variable solutions 
i the type (2.14) and (2.22) are important. Thus if 
y Le** or y= Fe*'x (3.11) 


; 1 criterion for determining whether or 


10t solutions oscillate indefinitely has been 
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respectively in an exponential or oscillatory region, where L, F, A, x are the solu 
functions of x, it is usually possible to choose these functions so that they} 





|, for v 
may be tabulated at much wider intervals in 2 than are convenient for the niv be | 
solutions y. In an exponential region the alternative auxiliary functions aches 

. | be b 
Iny = InL+A, y' ly L'/E+N’ (3,12\— ¢* = co 
are also suitable for tabulation at wide intervals. and | 
This paper is not concerned with the choice of normalization constants °"Y°™* 
(such as C in (2.23)). The main interest lies in the relative proportions of) the 
two preliminary solutions to be taken when defining standard solutions. a 
To su 
3.2. Exponential region, I(x) J (x) < 0. In an exponential region) wlytions 
as has been remarked, multiples of the two solutions (2.14) that are off fiom one 
comparable size at one argument x = 2, will be widely different at another 
iae aa. ¢ 
x= 2,if | | J*(t)dt is large. Thus in the general solution —— 
To s ; F 
, , F ons Lf € 
y A L(a)e Nwx)_| BL(x)e Ax). (3.21 
) be take 
Zz . 
' . ipproxy 
where L(x) ~ J-*(x), \(w) ~ | J*(t) dt (3.22 PI f 
; pair of : 
Lo 
‘ a , ; .., , | imagina 
and A 0, while A(x) is for the moment supposed to increase indefinitel as 
‘ oe Toe a. These 
the term with coefficient B will be negligible—whatever the value of A el 
ilthoug 


when A(x) is sufficiently large, that is, when a2 (> a») is large enough.f. 
as ; 4 ‘ ; . Pllary fu 
Hence it is desirable, for such an exponential region, that the soluti : 
Big . Pphase ¢ 
L(x)e-“ should be tabulated separately. By a similar argument, traversing} ° ls} 
; analysis 


should | 
As ar 


the region in the reverse sense, that is, taking a») very large and x 
so that A(x) is large and negative, and with B + 0, it follows that #/ 
solution L(x)e+* should also be tabulated separately. These two solution 


—_ 


which form a fundamental pair, are then sufficient to give any others. 

By this choice, when the argument concerned extends to infinity in one? the solu 
both directions, with A(~) increasing or decreasing indefinitely, the tabul'} pair, } 
tion of the solution that does not become infinite is ensured, and this ¢a propert 
be important in physical and other applications. If A(x) does not increa*} the fam 
indefinitely in magnitude, for example because the exponential region? and od 
finite in extent, the above argument may lose some of its force, so that, 1} x > 0, 
adequate tabulation, it is necessary merely to have A or B sufficientl} pair, ar 


greatest value of A(x) that can occu} 








small, the size depending on the 


y ° ° . ° * al + For 
Nevertheless, it is usually found convenient to tabulate the solutio) 
L(x)e**™ in all cases; these always give a satisfactory pair. > Mportan 


As an illustration, consider the equation 
? inthe 
yo AChE not 


y’—y = 0. (3.2 


Nn sectior 
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X al} The solutions sinh» and cosh are useless as a fundamental pair for large 


they} ». for when x is large and positive the solution e-* cosh «—sinhz can 
tthe ly be obtained numerically by keeping an excessive number of figures 


tions i) tables of cosh x and sinh while a similar remark applies to the solution 


(3.12 cosh a+sinh w when 2 is large and negative. In fact, it is the solutions 
and e-~* which must be tabulated in order that all solutions may be 


nveniently obtained to a specified number of significant figures through- 


eR t the range of tabulation. In this particular case, of course, owing to 
ions, | mmetry, only one of the solutions needs tabulation. 

lo sum up, in an exponential region, the choice of a standard pair of 
gion.) wlutions should be such that their ratio alters by « factor as great as possible 
ATE Ol from one end to the other of the re gion. 
other} 

} 3.3. Oscillatory region I(x) > 0. In an oscillatory region, two indepen- 

ent solutions in general retain the same relative importance. The solu- 
3.91 ' tions Fe*'x of (3.11) are conjugate, and the real and imaginary parts may 
""') be taken as a standard pair of solutions. They have the same modulus, 
> ot {proximately /~!, and differ in phase by $7, but they are not the only 
° air of solutions with these properties, which also hold for the real and 
. maginary parts of Fe*"x*©, for any constant e. 
i$ These properties are not particularly useful for table-making purposes, 
" ~ Salthough the possibility of tabulating the relatively slowly varying auxi- 
a | iary functions F(a) and x(x) is important. The choice of solutions with 
7 I phase difference lz is, however, often of considerable importance in 
7 analysis,t and it will be adopted as a desirable principle that solutions 
- | should be chosen to have this property. 
pri As an example. for the equation 
on y’+y = 0 (3.31) 
one’? the solutions cos x and sin x are normally used and are clearly a satisfactory 
abuli} pair. However, the solutions A cos(z-+e) and A sin(x- e) have similar 


118 Cal) properties and are apparently just as satisfactory. In this particular case, 
»TePAY tha f; 4h . a. . a 
creas} the familiar cosx and sina are chosen because they are respectively even 
gion} and odd; this means, incidentally, that they need be tabulated only for 





| 
at, I") « > 0. In not all cases, however, will even and odd solutions form the best 
cient} pair, and suitable alternatives may be used. 
ocel 
lution} For example, the standard solutions of Bessel’s equations, J, (2) and Y,(x), share the 
modulus and differ in phase by $7; the usefulness of this pair depends largely on the 
ance of the Hankel functions 
Hs,,(a J (a ) . Hi,,(a) J,(x)—iY,, (x) 
; t ° ° 
(3.2 ® notation of Jeffreys and Jeffreys (3). Compare also with the geometrical interpretation 
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[f, for example, /(x) is an even function, the criterion may be applied 
that the graphs of the two solutions should be mirror images in the y-axis 
that is, f(7) and f(—x) should be the fundamental pair. This again requires 
tabulation only for positive x. It will be possible to arrange that the phases 
difference is $7, though it may not be possible to have equal moduli as 
well. In the case represented by (3.31), if A N2, € lor, the resulting’ 

air is , ., ' , 
- f(x) cos x— sin a, ji(—2z) cos x-+sin x. 
To sum up, in an oscillatory region the two solutions chosen as standari| 


should have phase difference 4a and equal modulus. The choice of phase} 


constant will depend on the particular equation, but will often be such that 
either (a) the two solutions are respectively even and odd, or (6) when I(x 
is even, the two solutions are of the form f(x) and f(—2), in which case the 
requirement of equal modulus may have to be relaxed. } 


3.4. Further considerations. The criteria for choice of fundamental? 


solutions have been stated in a form that is not as precise as it seems. 


3.41. In the first place, relations (2.14) and (2.23) are only asymptoti 
approximations—good ones for the particular solutions concerned «s 
x — +oo in an infinite jnterval, but needing modification towards zen: 
of J(2), or when |/(x)| has a small minimum value. In fact, the solutions 
represented by these approximations as x — +00 and as x > —© are no 
in general, identical. Thus when /(x) = }x?-+-a, with a positive and smal 


a pair of solutions which have equal modulus and phase difference }: 


=o 


when x, will not have these properties as x > —x, although tl 
deviations are very small when ais positive and large. The phase constant 
(2.23) is then chosen, as it can be, so that the phase difference retains t! 


value 47 as x > —oo. The modulus ratio, on the other hand, cannot, | 
general, be unity both as x > —o and asx -> +00. If, then, it is importar'] 
for the standard solutions to have the form f(2) and f(—2), a normalizatio 
factor is chosen so that the ratio of the moduli of f(—) and f(x) varie 
from, say, k as x > —oo to 1/kas x” —> +00. In order to obtain the modult'| 
and phase from these solutions, only 2 > +00 need be considered, and tl } 
values are determined from the relations 
. l - ° 41] 
F cos x ps) Fsiny = vk. f(—2). (3.411 } 


! 
3.42. Secondly, it is not at once obvious how to define precisely th 


meaning of equality of modulus and a phase difference of $7. The form 
(2.23) is sufficiently definite, but is only an approximation to the exacy 
form CF cos(x-+-e) which may be derived from (3.11); in a finite, possibl 
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short, oscillatory region, this approximation may nowhere be good enough 
to define precisely the solutions needed. 

The difficulty is that the modulus and phase functions F and y are not 
unique, even after allowing for the constants C and e. For if y, and y, 
ire any two independent solutions whatever of (1.1) it is still possible to 
express them in the form 

Y; F cos x, Yo F sin x, (3.421) 
where F? — yi+y3 tan xy = Yo/Yy}- (3.422) 

Nevertheless, the criterion is not meaningless, for it is usually possible 
to pick out satisfactory functions F and y, satisfactory in the sense that 
they are free, or as free as possible, from the oscillatory fluctuations that 
occur in the general F and y, thus rendering them more amenable to 
compact tabulation; these particular functions, being relatively ‘well- 
behaved’, are also useful in analytical investigations. Compare, for 
example, the Madelung transformation (4). Thus, with equation (3.31), in 
which J (a) 1, take 


yy cos 2x, Yo = cos(x+a) (3.423) 
which leads to 
F2 cos2x-+ cos?(a+ «) 1+ cos a cos(2x-+-a) | 
: (3.424) 
tan x cos a—sina tan. J 


It is clear that a | tz yields the least oscillatory (here constant), and 
s0 most satisfactory values for F and y’. 
3.43. In order to study the modulus and phase functions more closely, 
consider the associated third-order equation 
$/2'4+-21'z = 0, (3.431) 


which is satisfied by the product of any two solutions of (1.1), so that its 


ve} 


eral solution is P 2 rT 
on ay; +-by, y.+cys (3.432) 


with the approximate solutions, corresponding to (2.14) or (2.22), 


z~J and J exp) +2 ( Jt) dt} (3.433) 
Or z~l and J exp| +2 ( nw atl. (3.434) 


Clearly L? or F? of (3.11) is a solution of this equation. 
ey . . . . . 
By analogy with second-order equations, or with third-order equations 
having constant coefficients, it may be anticipated that three fundamental 
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solutions of (3.431) can be chosen which will be either (a) all of exponential 
type in an exponential region of (1.1) or (b) two oscillatory solutions of 
equal importance, together with a third solution without oscillations, in 
an oscillatory region of (1.1). 

It now becomes clear that, for the modulus function in an oscillatory 
region, the non-oscillatory solution of (3.431) is needed. In any particular 
application it is usually quite easy to pick out the solution required. 


4. Geometrical interpretation 

It is illuminating to express the results of previous sections in geometri- 
cal terms. 

4.1. The solution Y = 9, +ig, Feix. (4.11) 
in which y, Feosy and y, = Fsiny are two real and independent 
solutions of the equation (1.1) gives the three-dimensional curve (x, F, x) 
in cylindrical coordinates. As an example of similar treatment of the 
Fresnel integrals, see (5) (third and subsequent editions), Fig. 14, p. 36. 

In an oscillatory region this curve is helical in character and winds 
round the solid of revolution obtained by rotation of F F(x) about the 
x-axis. When, as will be assumed throughout section 4, y, and y, are a 
properly chosen pair of standard solutions, F and x will usually be non- 
oscillatory functions, changing slowly and steadily in value when com- 
pared with the oscillatory solutions y, and y,. Similarly the solution 


Y* = y,—ty, = Fe-x (4.12) 
winds round the same solid of revolution, but in the opposite direction. 


4.2. Consider now the solution 
Y, = Ay, +7By, (4.21) 
in which A and B are real, unequal, non-zero constants. This likewise 
winds round a solid, not of revolution, but having similar and similarly 
situated elliptical sections in parallel planes, the semi-axes being AF and 
BF in the real and imaginary Y-directions. 
Another picture is possible and, although clearly less satisfactory for 
most purposes, serves to throw light on the proper choice of standard 


solutions. This is obtained by rotating the curve Y, about the x-axis to | 


give a solid of revolution round which the solution winds as in section 4.1, 


— 


but this time the radius F, fluctuates as x varies, particularly if A and B 


are very unequal, giving an awkward concertina-like solid with transverse { 


corrugations. In such cases z = F? is one of the unsatisfactory oscillatory 
solutions of equation (3.431). 


} 


. 
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ntial | 4,3. The general solution, y, and y, being independent, is given by 
’ « * 
is of }} (4.21), but with arbitrary complex values for A and B. Clearly 


Ss. 1 


, . ‘ P 
f 4(A+ B)(y,+-ty.)+3(A— B)(y,— typ), (4.31) 
ton hich represents the sum of two solutions winding round two solids of 
ule evolution, similar though of different sizes, the x-speeds of rotation of 

erpendiculars to the x-axis being equal but in opposite senses. 
Again, Y, Ce Y Feilx+4 C,e'” Fe (x+€) (4.32) 

ety. y Where C), Cz, y, and € are real constants given by 

2C, elly+ A+B, 2C,e-"y-9) — A—B. (4.33) 
, This represents a solution winding round a solid of elliptical cross-section, 
sinsection 4.2, but with major and minor axes, 2|C,+(C,| F and 2\C,—C, | F 
lent | at an angle y with the real and imaginary Y directions; the sense of 


Fx) # rotation depends on the relative signs of C,+C, and C,—(C. 


the 
4.4, It is now readily seen how any general complex solution of (1.1) 


we | may be made to yield a satisfactory pair of fundamental solutions. First 

the multiply the solution, say that given by (4.32), by the constant e-'’, thus 

pee bringing the major and minor axes of the elliptical sections of the asso- 

sll ciated solid into the real and imaginary Y directions. This gives 

om- # ev, (C\+C,)F cos(x+e)+1(C,—C,)F sin(x+e). (4.41) 
} Next alter the scale of real, or of imaginary parts, or of both in order to 

+12) | convert the associated solid into one of revolution. This gives, say, 

mn. : y. CF cos(x +-€)-+ iC F sin(x t-€). (4.42) 
, The real and imaginary parts then satisfy the requirements proposed in 

91) | Section 2.2, that is, they have equal modulus and phase difference 37 

_ | The normalization constant C and the phase constant « remain to be 

wr decided on the merits of each particular case concerned. 

arty 

and } 4.5. As an example of the application of these geometrical representa- 
} tions, consider the choice of fundamental solutions for an equation having 

for two oscillatory regions R, and R, separated by an exponential region S 

lard | (or by a region where solutions oscillate slowly, so that approximations of 

s to } type (2.22) or (2.23) are not good). 

4.1, | Choose solutions », and », in R, of the form F cos yx and F'siny. Then 

d B i+in, = Fe'x, after passage through S, will emerge in R, with a form 

erse such as (4.31) or (4.32), winding round a solid with sections that are, in 

tory } general, elliptical, possibly with a large ratio of major to minor axis. The 


inclination y of these axes to the real and imaginary axes of Y must then 
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be determined, and a factor e-'Y applied, as in (4.41), so that the elliptical) Als 
section have their axes in the real and imaginary directions. The real and 
imaginary parts are then }z out of phase, while the real part, corresponding 
to the major axis of the elliptical sections, is as large as possible, and the) when 
imaginary part as small as possible for given F and various phase constants | 


in #,. Thus the solution oo th 


tix, = FeX-Y = en, +ing) 

gives a standard pair of solutions y, and y, of the differential equation 
. : ~¢ : . when 
considered which have phase difference $7 in both regions R, and R,, } 
whilst the moduli are equal in R, and have the greatest possible ratio to | 
one another in R,. Fr 


5. The case of a finite oscillatory region are 4 
purp 


In a finite oscillatory region it may be difficult to attach precise meanings ' 
the 1 


to equality of modulus and a phase difference of 4, though this is simpk | 
when asymptotic developments in infinite regions are available. In othe 
words, these properties are, in a sense, asymptotic in character and, a: oie 
with asymptotic series, precision is not easy to obtain when the argument | 

is not large. For purposes of tabulation it is often important to obtain good 


approximations to F and y, in order that tables may be made that are 
interpolable at comparatively wide intervals, particularly when, as with | 
the Weber equation, y”—(}a?—a)y = 0, there is a parameter a as well as 


the variable x, and two-way interpolation is needed. If, on the other hand, | 


wre 


tabulation or theory involves only the solutions of the differential equation 
in the usual form, and a close interval of tabulation is permissible, it may 
suffice to use any two convenient solutions. 


n> 


One case of fairly frequent occurrence seems worth special mention. 


Consider y” +I (a,x2)y = 0, (5.1) 


in which J(a,2?) is symmetrical about 2 = 0, where I(a,0) > 0, so that 
an oscillatory region surrounds the origin, and where a parameter @ occurs 
which, when large, is such that changes in J are not appreciable for several 
wave-lengths of the solutions. It is then possible to determine modulus 
and phase functions F and yx at x = 0 to a good approximation from a 
single solution. Suppose the known solution to be 


y, = Foosy (5.2) 
in which F and y are to be determined. Then 
(ie) = — F(0)x'(0)sin x(0) (5.3) 
dx} .-9 
since, by symmetry, F'(0) = 0. (5.4) | 
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Also, since Fe'x is a solution of (5.1), it follows by substitution that 


P’—Fy?+IF=0, 2F’y’+Fy’ =0, (5.5) 

a oly iy’ a 

whence ** ; (*] v= x =f (5.6) 
e\x/ «6X 


” 


so that, to a first approximation, y’/y’ and y”/y’ being assumed small, 


x/(0) ~ y{L(a, 0)} (5.7) 
whence tan y(0) ~ (y; tn)z-01 (a, 0) (5.8) 
F(0) ~ .{y2(0)+-y32(0)/T(a, 0)}. J 


From (5.4) and these approximations for F(0), (0), and y’(0), which 
are asymptotic in character as a becomes large, but adequate for most 
purposes, the modulus and phase functions may now be determined from 


the relations, derived from (5.5), 
F?y’ K, F’+IF K?2/F8, (5.9) 
where AK is a constant, here approximately 


y3(0)/{L(a, 0)} | y7(0) V{T(a, 0)}. 
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\ SHORT TABLE OF THE CONFLUENT 
HYPERGEOMETRIC FUNCTION (a, y, x) 


) 
By B. W. CONOLLY (Admiralty Research Laboratory, Teddington) whe! 
| 
[Received 15 February 1949] A 
som 
CONFLUENT hypergeometric functions have been extensively studied 
theoretically, but, in comparison with allied functions which arise naturally T 
in mathematical physics, they have received little systematic numerical | 
tabulation. This is perhaps due to the greater number of variables to be | -" 
considered. What seems, however, to be a stronger discouragement to the } _ . 
computer is the multiplicity of ‘fundamental forms’ of the functions in | ili 
existence. 
These fundamental forms arise in two ways, either as most appropriate 
to a particular practical problem, or as giving symmetrical formulae in a 
general theory. To quote only two instances, confluent hypergeometric by 


functions appear in the mathematics of nuclear physics and also in wave- | PY° 
propagation problems where the boundary is a paraboloid of revolution. | "' 
On the other hand, theoreticians concern themselves mainly with certain | © 
canonical forms due to Whittaker (ref. 1, chapter 16). che 


The table appended here is of what is sometimes called the Kummer , 
| Ma 


function M(a,y,2). It is also frequently denoted by ,F,(a;y;x). Of all the 
confluent hypergeometric functions it has been the most extensively tabu- to1 
lated—by J. R. Airey and A. J. Thompson (ref. 2, para. 22.56). 
. . . a ‘ ( > 
Provided that y is not an integer, M(a, y, 2) and #!-YM(a—y-+1, 2—y,2) le 
form a fundamental system of solutions, valid near the regular singularity 
« = 0, of the differential equation 1. 
da J di x 
= 4 EB y 0. (1) 4 2. 
dx?  \x de x 
This equation may be readily obtained by permitting the regular finite 
non-zero singularity of the differential equation 
dy di , 
a(1—a) <4 4 fy—(a + B+ 1)x}F — apy = 0, (2) 

dx? dx 
satisfied by the hypergeometric function of Gauss, ,F\(«, 8; y;x), to become 
confluent with the regular singularity at x = «©. It will be seen that the 
singularities of (1) are a regular one at x 0 and an irregular one at 
x ©. The process of confluence applied to (2) lends its name to the 
solutions of the resulting equation (1). j 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1950)] 
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From (1) the absolutely convergent series for M(«, y,z) may be obtained. 


! Itis J 
M(a,y,2) = > a,x’, (3) 
r=0 
‘ C(y) M(a+r) 1 
where r 7 ; a? 
my D(a) (y+r) r! 


An important formula due to Kummer which has been used in obtaining 


some of the numerical values given is 


died M (a, y, 2) e*M (y—a, y, —2). (4) 
am The tables give 11 decimal values of M(a,y,2) for x = 0-1, 0-2(0-2)1-0, 
ahs ’ 1-0(0-2)1-Oand y = 0-2(0-2)1-0. Since M(a, y, 0) 1, M(0,y,x) = 1, 
0 be 


M(x. «x, 2) e”. certain entries could be made without calculation. 


) the The coefficients a, were provided by Dr. J. C. P. Miller and used to 
Is il is ; 
obtain values of M(a,y,a). Since 
riate et ¥ a,(—x)" = e*M (a, y, —2) M(y—a, y, 2), 
ina ins 
oiete by (4), it was possible to obtain other function values from the same 
ave. | products a,a". As a check, M(y—a,y,x) was calculated without reference 
bias to the a, supplied by Dr. Miller, agreement implying the correctness of the 
hain corresponding .M(a«,y,x). Likewise, all the other function values were 
checked independently. 
— M(x, y, 1-0) has been taken from an unpublished table by Dr. Miller. 
| the Many thanks are due to him for allowing its reproduction here and also 
“_ for permission to use his coefficients a,. 
The values are believed to be correct to within 2 or 3 units of the last 
“) decimal place. 
rity REFERENCES 
1. E. T. Wairraker and G. N. Watson, A Course of Modern Analysis (Cambridge, 
1940), 
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Tables. London, Scientific ( omputing Service, Ltd., 1946. 
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SOME EXTENSIONS TO THE METHOD OF 
INTEGRATION BY STEEPEST DESCENTS 


? 
By P. C. CLEMMOW 

(Department of Electrical Engineering, Imperial College, London) 
' 


[Received | November 1949] 


SUMMARY 

| In practical cases the assumption of small values by certain parameters may 
lidate the usual method of integration by steepest descents when it would be 
therwise applicable. The two ways in which this can happen are treated by using 

} the concept of a partial asymptotic expansion, the terms of which are functions, 
rather than inverse powers, of the relevant variable. It is shown that the theory 
f Watson’s lemma, and a device due to Jeffreys for estimating the ‘remainder’, can 
e extended to these partial expansions. The general terms are expressed in the form 
f Hh, or confluent hypergeometric functions, and the leading term may reduce to 


mplex Fresnel integral. Two examples are given from the theory of wave 


propagation. 


7 


1. Introduction 
THE method of integration by steepest descents is applied to contour 
integrals like . 
| d(z)ek¥ dz, (1) 
; 
and yields an asymptotic expansion in inverse powers of k. By distorting 
C through the appropriate saddle points of #(z), (1) is transformed into 


a sum of integrals of the type 


| f(Aje-** da, (2) 
0 
and the procedure then depends on two steps. First, (2) is replaced by 
| f(A)e-** dr; (3) 
0 
secondly, Watson’s lemma is applied to (3). 
Ifa and f(A) are fixed the asymptotic expansion thus obtained is always 
valid for sufficiently large values of k. But in physical problems the value 
* of bw . : 
ol k will be restricted, and the steps in the above method may break down, 
respectively, if 
(I) | is small, since the passage from (2) to (3) neglects terms con- 
taining the factor exp( ka"): 


Il) the radius of convergence of f(A) is small. 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 2 (1950)] 
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In this paper cases I and II are treated by using the concept of a partial 





asymptotic expansion, defined as follows: if } 1, The 
n , 
fle) = > fale) +R, (2), yf 
m=0 
where , 
where, for all n, z”R,,(z) > 0 as |z| +00 and for values of argz within | 
given interval, then ? 
fol) +hl2)+ fol) +- | 
we Car 
is called a partial asymptotic expansion of f(z) for the given range of arg:. 3) is 
A development such as (5) is clearly not unique; but if each f(z) is | 
replaced by its asymptotic expansion (in inverse powers of z), the unique 2! 
asymptotic expansion of f(z) will be recovered. Alternatively, (5) may le / i 
regarded as formed by a particular grouping of the terms of the asymptotic és = 
expansion of f(z); different groupings give different partial expansiots, | me 
The idea of reaching a useful result by regrouping the terms of a valueles 
expansion is similar to that in Airey’s technique for improving th 
accuracy obtainable from an asymptotic series (1). = 
Case I is treated by obtaining a partial asymptotic expansion of 
: ; and 
| f(Aje*™ da. (6) 
, In (12 
This embraces the development in inverse powers of ka, obtained by| 
continual integration by parts, in the manner just explained, but also’ 
reverts to the normal asymptotic expansion of (6) when a = 0. Case I] 
is more complicated: difficulties arise in showing that the partial expansion) | 
' . ie a hen 
given is always useful, and the general terms cannot be specified for an| 
arbitrary f(A). When the small radius of convergence of f(A) is due to 4 
factor (A?+-z)-” the general terms are confluent hypergeometric functions | 
This was shown by Pauli (2) in considering diffraction by a wedge, where} 
the singularities of f(A) are simple poles on the imaginary axis. Ott (3) has 
also applied the theory to problems in which the singularity is a genera 
complex pole, but there has been no discussion of the nature of thi _ 
expansions used. vee 
In the next section two aspects of the theory of asymptotic expansions | 
are considered, namely Watson’s lemma and a ‘remainder’ estimate due For 
next, 


to Jeffreys. This paper shows that the scope of these may be extended to 
partial expansions. Case I is treated in section 3, case II in sections 4 and 
5; next, both methods are applied to the same integral, and their effect i! 
combination is then noted. Finally, physical applications are given. 


































METHOD OF INTEGRATION BY STEEPEST DESCENTS 


‘os " ul . y > 
) ). The asymptotic expansion of | /(A)« kr dy 


Let fA) = A-* dc, APs, (7) 
$) s=0 
here B is real and positive and rea < 1. Using the result 
phir 
: l/p—l1\, (41) : ' 
, | APe-* da 5 - )!q p+), for rep> —l, (8) 
e can prove by Watson’s lemma (4) that the asymptotic expansion of 
argz.} io) ; 
6” 0) 18 
ae Q ee 
n\2 l 1}. aoe 2}, ¥ ] 1) ” (" X ] Ii: $B) Co 2B x ] tp-B_|. (9) 
cr) a a a A, la 
ay"! In what follows we shall for convenience take « = 0, B = 1, but these 
tot ‘ere —_ , » , 
is restrictions are not necessary. To get an estimate of the remainder we 
S1LONS ° . 
SOFT use the Maclaurin expansion 
ueles n=l 
g tl f(A) 2% TT ys (10) 
* where c, = f(0)/s! (11) 
d 
l e : 
and - , (A- -t)"-1 f(t) dt. (12) 
( (e—i)! . 
In (12), put ¢ A/(1+ p), so that 
ad bi 
als ] An n-l : A 
t usc , ns fm : dy. (13) 
ase I] (n 1) J (1 1 y)e+t l+-p 
nsionp : 
ee hen from (10) and (13), using (8), 
> to . kd? ] ~~ . ‘a—~1) 1 [-—H8+1) R 
tions ] f(Aje™" dn 2 ha “s 2 ) fia , y 
where } ie 
n—1 
3) has > f(k)+R,, say, (14) 
nera s=0 
yf the : | a A" 1 ; A . 
vhere R | fm e-*® dudd. (15) 
(n ty. (1+p)"*?° 1+ 
1s10Ns 0 0 


; ‘ . ° . . 
e due For the A integration in (15) change the variable to z, where A — (1-44 pz; 
led to| Next, for the » integration, change the variable to y, where p = y/z. Then 
4 an 


ook Se 1 “? 
ect | R y” 1f(™)(z)e k(z+y)* dydz. (16) 
. (n—I)' J J ° ° : 
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The y integral in (16) can be expressed in terms of the parabolic o_— Substit 


function (5, p. 349). For our purposes it is convenient to use the form 


if er Rk, 
Hh, (z) yreic+y’ dy, for ren > l, (17) 
nL. 
0 | 
so that (16) may be written 
= ) 
R,, = (2k)-" | f™(z)Hh, _,f2v(2k)} dz. (18) | 
) 
0 


The expansion (14), with R,, given by (18), results directly from successive 
integration by parts (6, p. 593), but the method used here is more easily If we ¢ 
adapted to later requirements. | of zin 

An estimate of (18) now to be given depends on n being large, whichis / 
the case if /-! is small compared to the radius of convergence of f(z). We 
use an approximation to Hh, (z) valid for n large and z not large, whic wsilies 
together with some assumptions about f(z), enables us to evaluate (18 
One weakness of the method lies in the use of an approximation for th 
integral which only holds over the initial part of the range of integration, } 
for though the justification is that the important contribution to the) The v: 
integral comes from that sector of the path near z = 0, this itself is} in (23) 
scarcely evident without a closer examination of the form of 


f(z)Bh,,_,{2v(2k)} by (29 


for n large and z ranging from 0 to «0. Another criticism is that use is} appli 


must 1 


made of methods allied to that which we aim to establish. The argument 

is thus inclined to be circulatory (7, p. 416), but seems reasonable if it is/ 

regarded as similar to that of successive approximations. Next. 
The required approximation (correcting a slip in (6), p. 592) is 


: Qn 1)/2 n ] i re 
Hh, (z) = = ( - Jie evn —2*/4 (19 With 
Initially we shall assume ) 
‘(z z- *pia\—m . | z— mpi 4 9() ; 
f(z) = (2—re'®) or log(z—re'*) ( Putti 
but afterwards the effect of a more general form will be studied. Then 
-\"(m+n—1)! ] 2. new 
J™(2) ( ) ( : ) : f™(0) i — . (21 
(m 1)! (z— petx)m+n ‘ re'® 
In considering (14), note that 
where 
"(] P(s) 0 oie Qi! HEP Sa 8 ' 
Salk) f (0) {(s—1)/2}! 1 Wakes 1 /s | for s large; 
fe-s(k) = f®-2(0) {(s—2)/2}! svk ret 2 svk 
so, if m is moderate, the smallest term is reached when and 





nm == 2kr?. 
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inder | gybstituting from (19) and (21) into (18), we get 


rm ZL 
9 n—2)/2 n 2 - 2 (m+n) : 
R, = (2h)-r/2f C00) ) | (1 = ; e—2 (2k(n—1)} —kz*/2 dz 
7” (n—1)! 2 : re’ 
(li ? 0 
m—~—n ms | z 
f,,-1(*) exp (m+-n)log| 1———} — 
re’ : | ret 
12 0 
—2,/{2k(n—1)}—kez? 2 dz. (23) 
ssive } 
easily | if we expand the logarithm in (23), and ignore second and higher powers 
| of zin the exponential, we have 
ich is ? . 4 m+n—1 : 
R, = = i. i(k); (24) 
\ e'%, {2kr?(n—1)}—(m-+n) 
‘hic } and using (22) this becomes (cf. 6, p. 594) 
(18 
; ] ‘ ™ 
r Ry = aa Iu-all. (25) 
Lt1ol 


o the) The validity of (24) requires that the coefficient of —z in the exponential 
elf is} in (23), namely, 
' /{2k(n—1)}—(m-+n)/(re), 

must not be small; and the method thus fails for « = 0 when n is given 


by (22). But a smaller value of » will indicate that the expansion is still 
Use Is 


applicable, though less accurate. If n = kr? 
‘ment : 
f it is a Aae i fy —y(k). (26) 
Next, consider 
J(z) = (2—1, e™)-™(z—1, e&)-™, where’ 17, > 1. (27) 
(19) With a slight change of notation, write 
} n =% 
; (A—r, e*™)-™ g(A) = > c,A®+r,. (28) 
(20 : s=0 
Putting nN 2h, i, we have 
en 
° : vt 5 . 
(21 | f(Aje** dA = > f,(k)+R,,, (29) 
0 8=0 
) 
w] pte A 
where now — 
: Jatt ' | (A ry (Ot2)me dA, (30) 
0 
‘ : l P n 14(n) 2 , 
and R te ( ) e~ke+u? dydz. (31) 
(22 (n—1)! J J (e+y—ret)m ‘ 
oe 0 0 
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An approximation to (31) can be obtained in a way similar to that for (16), 


In the factor (z+-y—r, e’™)-™, we can put y = ,/{n/(2k)}, z = 0; therefore 


(31) only differs from (16) by the factor 


(,/{n/(2k)}—rq ef%2)-m= = (1, —1, ef02)-m, 


(32) ? 


Furthermore, if 7, is appreciably greater than 7,, the factor (A—r,e!)-™ | 


in the integrands of (30) can be expanded as a power series in \ and tem | 


by term integration carried out, the remainders being of smaller sole 


than #,. Substituting these expansions for f,(/) into (29), an asymptotic 


development is obtained which is evidently arrived at by using the 


| 
| 


ascending series for f(A) and integrating term by term; and we note that 


the smallest term and order of the remainder are determined by thie | 


singularity nearest the origin. If, however, 7, = 7, the situation is compi- } 


cated, since each f,(/) in (30) has a remainder term of the same order xs 


f,,, whose contribution would therefore have to be estimated: for this cas 


aac A = } 
difficulties also occur when a, = 0. Finally, we remark that the abov 


observations are equally applicable when, instead of two factors, as i 


(27), f(z) has an arbitrary number of factors. 


3. A partial asymptotic expansion of e*@ | f(A)e-** dd 


For convenience the integral under discussion will be taken as 
‘. 


eka® [frye ki? dy. 


a 
To obtain an asymptotic expansion of (33), write 
P(k, A) = —f(A)/(2kd) 
and integrate by parts » times. This gives 


x 


> (—)"F, (ke, a)+-(—)me*™* | Fi (k, Ale dd, 
m 1 ’ 


a 


where the dash denotes differentiation with respect to A, and 


F.(k,A) = —F,_,(k,A)/(2kA), for r = 2 to n. 
From (34) and (36) F.(k, A) = (2kA?)-" F(A), 
r—} 
where F(A) = (—)A } (—)8*th,., A8f(A), 
and oe 1 
h,., 2 3r(r—1) 


h,p-3 = §(r+1)r(r—1)(r—2) 


h,o = 1.3.5...(2r—3) 


(34 


(36 
(37 


(39 


Assun 


provi 
inflex 


To 


write 


Bro 
pr ce 
the m 


be wi 


Let 


Then 
for a 


for al 
lor n 


Equ 


asyn 


whe 


for . 
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r (16). Assuming that f(A) is appropriately bounded, 


refore 
kneka® | F’(k,Aje** dA>+0 as k>oo, 
(32) | 7 
provided that boa < arga < 4m and the path of integration has no 
1X2 | . . 9~ or . . "06 . 
) inflexion. From (35) and (37) the asymptotic expansion of (33) is then 
1 tem 
P]-72 7 I)-q2\-2F 9}-72\-3 F 
ee (2ka*)-1F, (a) +-(2ka?)-*F,(a)+ (2ka?)-3 Fy (a)+.... (40) 
ptotie} To obtain a development for $7 < arga < 3m it is only necessary to 
ig the} Write (33) as : 
e that eka® | f(rAje** dA—ek® | f(—A)e™ dd. (41) 
| : “a 
yy the 
ompi-} Broadly speaking (40) cannot be useful unless |ka*| is large. We now 
der | proceed to a partial expansion which is not subject to any limitations on 
seas | themagnitude ofa. Ignoring the possible capture of singularities, (33) may 
abow | be written 
as i | f(A+a)e-2har—ka* dq), (42) 
0 
Let f(A+a) = > c¢,A*° for |A| <r. (43) 
S 0 
Then, if real positive m and d exist such that A~"e-“"f(A+-a) is bounded 
for all real A y, — 
f(A+a)— Yc.A% < KA"e™ 
99 s=0 
5D) 
ralln > mandallA > 0, where K is some fixed positive number. Hence, 
m, 
24 . ls ? 1 
= e~<aKa—l f(A+a) > 4] dx 
K | A” le 2kad | p—(k—d)A* Jy 
(35 
0 
K | Ane-“*-O" dr, for —4n < arga < }z, 
(36 0 
9- 1 -(n 1 . 1—-(n +1) 
(37 }K{——}! {y(k—d)} , (44) 


(38), Equation (44) proves that (42) has, in the sense of (4) and (5), a partial 
asymptotic expansion 


folk) +fi(k) +fo(k)+..., (45) 


(39 where fk) = cg | A8e-2akre-k4* a), (46) 


or —4, : io (otherwise we use (41)). 
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If k is fixed, a may be arbitrarily small in (46) without upsetting its 


? 


validity. For a = 0, f,(k) is the same as in (14). The leading term in (43) | 


may be written as an error integral or a complex Fresnel integral. With [ 


applications to wave propagation in view, we choose the latter, adopting 
the notation : 
F(a) = | eda, (41) 
a 
so that folk) e'4F(a)k beka* F(e-‘4qwvk). (48) 
The Fresnel integral is tabulated for real values of the argument (8, p. 3 
9, p. 571), and for a restricted range of complex values (10). 
From (17) and (46) the general term in (45) is 


f.(k) f(a)(2k) (s+D/2eka" Hh (av2k). (49 

This gives the asymptotic expansion 
| as+1f@(a) {, (s+-1)(s+2) (s+1)(s+2)...(s+2m) | 
f.(k) ~ : _ +. m. +...) 
s\ ) (2ka?)s+1 | 2(2ka?) ( ) 2mm! (2kh« ve | 
(50 


a 


5° | 


Substitute from (50) into (45). Only those f,(k) for which s < r contribute ¢ 


to the k-" term, and it is readily seen that summation over these values 
of s expresses this term in the form 


1) r—1l, * ian 9\! 
(—7F "9 S ( — 2a)* pon wi f(a). (51) 
ai r = s!(r—s—1)! 


s 


A comparison of (51) with (38) with A = a, gives the general relation 


(52) 


r,s 


: ay 1a Cr 


which may be checked against (39). To estimate the remainder in (45) we 
proceed as in section 2. We have 


L 


n—1 
aly 1 ~ 
iP} lq *\ . As tyr 1f(n) at dt. (53) 
mero 2 ‘pm ee 
“5 0 
Then the same substitutions as before give 
n—1 
| f(A+aje-2kar—-K¥* dy = ¥ f,(k)+R,, (54) 
. s=0 
where R. = (2k)-**e** { f™(z+a)Hh, | (2+-@),/(2h)] dz. (59) 


0 


Using (19) and (49), for large s 


fh) = 4f(a) SEE (by e+e avant stare; (56) 
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Ng its f.(k) f(a) ¢ av {k/(2s)} o 
y 0 that = —— ° (57) 

n (45 fe-1(k) ° ff (a) 4/(2ks) 

With If f(z) is given by (20), (57) becomes 

pting } ian — 
f.(k) : m (1 cat ° € avy {k/(2s } (58) 
f.-1(*) rex ,(2ks)\  re® 

(47 


ind since we are only interested in small values of |a|, of the order k-}, 


ind k¢y is large, the numerically smallest term is again given by (22). 





(4 pe ‘ as ‘ 
Then from (55), using (19), 
ee ((n—2)/2}! , si casita and 
R : 1(}:) n f™(a)e ka*/2—av {2k(n—1)} x 
is (n—1)! ‘ 
bd ~ (m+n) 
(49) } (1 @ e-2s (2k(n—1)}—akz—kz*/2 qs 
f re'*__q 
0 
m——-n l a 1 i { , z 
' = f._,(k) (1. exp| —(m-+n)log|1—— ete 
", hae re 3 re'*_q 
0 F 
(50 ‘ 
—2x(./{2k(n—1)}+-ak)—ke? 2| dz, 
ibute ¢ : ro . . — e . 
; by using (56). Evaluating this in the same way as (23), we find 
alues * 
i-th ] > 
R ( fn-ik), 
(re'*—a)(, {2k(n—1)}+-ak)—(m-+-n) 
(51 which, with x 2kr?, and a/r small, again gives (cf. (25)) 
R tf _slk) (59) 

a “n pia Sar =—s ‘ia 
(359 The concluding remarks of section 2 also apply in this more general case. 


») Wwe 4. A partial asymptotic expansion of | f(A)e-**" dA 


0 
When the radius of convergence of f(A) is small, write 
(53 f(A) = g(A)h(), (60) 
where h(A) contains those singularities of f(A) which are near the origin. 
Let x 
g(A) = ¥ eX. (61) 
(54 ‘ne sis 
Then the required partial expansion is 
(55 fl) +Alh)+fe(k)+...; (62) 
where B C, ( ASh(A)e kX? dx. (63) 
(56 Ifh(A) is unspecified, it does not seem possible to give a rigorous justifica- 


tion of (62) by a treatment similar to that of Watson’s lemma (cf. section 3). 
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In the next section a particular form of h(A) will be examined, but first it} Then 
' 


is worth while considering the remainder estimate for the general case. 
By the same method as before 


r e-3 
| f(Aje** dA = > f.(k)+R,, (64 
0 s=0 


where now 


R,, | | y” 1g (z)h(z | y)e k(z+y)* dydz. (65) \ 


(n—1)! 


0 0 


If (65) is evaluated in the same way as (31) we have 


h|/{n/(2 


| 
(n—1 


ky pe 
; | | y"19(z)e-Ke+w* dydz. (64 
0 0 


Furthermore, normal saddle-point technique applied to (63) gives 


fa(k) = ceva s8?(2k)-+De-sl2H/ /{s/(2k)}] 


e—T] " ' 
i, = Jb) 6+D/2h/ ,/{s/(2k)}], (67 
when s is large. This shows that the smallest term is again reached when 
n = 2kr?; thus (66) again gives the result 
R = I ' 
a = ela <a Su i(h). (68 


If in each f,(/), h(A) were expanded as a power series in A, term by term 
integration effected, and corresponding inverse powers of k collected 


together, the asymptotic development (9) would be recovered. This repre- 


sentation breaks down because of the extremely rapid increase of h'*)(0 
with s, and it is worth noting that in our treatment of the partial asymp- 
totic expansion this difficulty is avoided by withdrawing h(A) from under 
the integral sign at A = ,/{n/(2k)}, which is not small, instead of at A = 0: 
in fact, the greater is n in R,, and f,,_,(k), the further does the saddle-point 
retreat down the real axis, away from the origin and the singularities of h(A). 


5. A partial asymptotic expansion of | g(A)(A2-+-2)-"e-*" da 


Some results for the confluent hypergeometric function (5, chapter 16) 


will be required. The most convenient form for our purposes is 


U(x, y, 2) ; ya" dt, for re(a—y) > 1. (69) 
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Then U(x, y,2) = 2!-7U(a—y+1, 2—y, 2); (70) 
Lhp- ne(” it (p p “Sk for ren ] (71) 

ind e*/2Hh,(z) = 2 vey (* sae 9 5) (72) 
The asy:nptotic expansion for —a < argz < vis 

U(o dues =) 


z 
x(a+1)...(a+m—1)(y—a—1)(y—a—2)...(y—a—m) — (73) 
m! 2™ oe i 

Series in ascending powers of z are more complicated, particularly when 
y isan integer (11). A unified treatment can be based on Barnes’s repre- 
sentation (12) 

. l l P 

U(x, y, 2) (st+-a—1)!(—s—1)! (—s—y)! 2° ds; 


2m (x 1)! (a yy: . 


the path is completed by an infinite semicircle on the right, and the 
residues of the enclosed poles evaluated. When y is not an integer these 
are all simple poles; otherwise double poles occur. In the following series 
the factor preceding 2” under the summation signs is to be taken as unity 
for = 0, and F(r) is the digamma function. 


Cas 


when y not a positive or negative integer or zero: 


ene = St > et eet 5, 
(a—y)! r!y(y+1)...(y+r—1) 


(y 2)" 2 > (a EBD 2)..(a—y- r) or. (74) 
(a—1) a 7 (2—y)(3—y)...(r+-l—y 
Case whe ny | 
U(c,,2) aN DV ala re t—)) tog 2+-F(a- r—1)—2F(r)]; 
{a . — et ag 


Case when ya positive ante ger gre ater than 1: 


U(a,y,2) = Y= 2) ay ¥ ye (a—y + Ia—yt2)..(a-y tr)» 
(a—1)! os r! (y—2)(y—3)...(y—r—1) 
(—)” S o(a+1)...(a+r—1) .,. 
La, y+l1)...(y+r—1)~ 
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Case when y a negative integer or zero: 


t <4 & 7 
Diese) = Ot SF let Mletr—D) 


(a—y I ‘ r!y(yt+1)...(y+r—1) 
: (—} al-y — (a—y+ 1)(a—y4+2)...(a r+)», 
(a—1)!(1—y)! L  ri(2—y)(3—y)..(I+r—y) 
<[log 2+ F(a—y+r)—F(1—y+r)—F(r)]. (7) 
Now in (60) take h(A) = (A?+2)-?, (78) 


where, for convenience, we suppose p (but not z) to be real. Then, fron 
(63) and (71), 


fAk) = $c,k” won? ‘}! U(r.r—"5 . : kz) (79) 


It is required to prove that 


krl2 5c MdA>0 as k>oo. (80) 
0 x 


Three cases must be treated separately. 


(a) p< 0. Here the integral in (80) is less than or equal to 


es n 
AY _ 4-20 
J (A?+ |2|)? 
0 
1 /n—1 ; —] . 
jpn ni (S jie (p.p— "5 sk), from (71), (81) 
l \ . 
os 5(" 9 ! z|-Pk-@+DI2 as k-+>oo, from (73). (82) 


This proves (80); we note that we may take z = 0 without invalidating 
the result, for then, using (74) or (77), (81) becomes 


1 /n—1 
("5 p)! ke (n+1)/2_ 


(6) p > 0, rez > 0. Writing argz = 4, the integral in (80) is less than 
or equal to 
¥ An 


2 -— e-kX* qj, (83) 
(A2+- |z|!cos 8)” 


so that this case is similar to p < 0. We can only take z Oifn > 2p—l; 
otherwise, note that, from (74), (75), and (76), the leading terms of the 
series for U(a,y,z) are 


4 = 
- -logz, for y i. "tie zi-y, for y>1; 
(a—1)! )! 
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so that, if k\z| is very small, the f,(k) given by (79) decrease in order of 
nagnitude as s increases. 

) p 0, rez 0. Let 6 ta+2e, where 0 < € < }x. Distorting 


the path through the angle e, the integral in (80) becomes 


in+Leh\n P 
A p—ke*€A? dr vn 2Pe kA® cos 2e dX 
Nezie + 2\ 
(77) 0 0 
- } Dy - 
(7/4) l (nm <p l ! (k cos 2e)” (n+1)/2 
9 9 ‘ 7 
, tror 


While this proves (80), it does not indicate that the expansion will be of 
any use, for a fixed value of k, if « is too nearly 47. For this case it seems 
that we must rely on some such approach as that given in section 4. 


A case with important applications is that for which p = 1. Then 


s—l1\,,, 3—s , 
al f,(he) = 4c, i! + 5 U1, ke 


—  ' A 2 . “7 
c( 1 (s—1)/2)41—s)/2 -e-k dd, from (71), 
2 (A2+-z)@+b2 
r 
(81) . 2 © 
; where \z is taken with a positive real part. For s 0, the integral in 
(82 . ‘ , 
'' (84) is a Fresnel integral. 
lating 6, Further remarks 
In (33) the substitution A? = p?+-a? gives 
* uf 21 @?)! 
Hw (H - My kp? dy, (85) 
(u2+-a2 
than 0 vere) 
to which the method of section 4 could be applied. The result would, in 
(33) general, be different from that obtained by the method of section 3. For 
example, conside1 
oki 
p—1; al sa ™: (86) 
7 T b? 
f the ‘1 


where byk is large. From (49), the general term of the partial asymptotic 
expansion of section 3 is 


fk) = f(a)(2k)-6+P eke" Hh fa,/(2k)}, (87) 
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where f(A) = (A?-+-b?)-!. But (86) may be written The 
2 the sh 

| — eT el (83) } region 

. (A?-+-a?+-b?),/(A?+-a?) the aj 





and the method of section 5 gives a partial asymptotic expansion with the result 
general term 
Y/(]. ( )° +\—s-1/2 le! J7(1 1 @ Jyy2 , 
E(k) (a2 pays *) S12 38! U(3,3—8, ka?). (5!) | hick 
itt 


For a ~ 0, (87) and (89) are different, except for s = 0. When a a 

however, f,(k) = F,(k) for all s. _— 
In section 3 the radius of convergence of f(A) was assumed not to le 

small. Now suppose otherwise, and consider aa 


a 


eka® [ g(A)h(A)je-*** da, (9 


a 
where both a and the radius of convergence of h(A) (but not of g(A)) may be | 
small. A little consideration shows that a combination of the methods of) The 
sections 3 and 4 will be formally applicable. In fact, if ' 

g(A+a Yc, As 
g(a ) a glves 
then (90) has a partial expansion with the general term 


ra 


¢. ( h(A+-a)Ase—2kad—-k% dy, (91) 


z whe 


7. Physical applications N 
The application of Fourier or Laplace transform theory to problems of a = 

wave propagation gives rise to the type of integral discussed in this paper app 

(see e.g. 7), which is therefore frequently encountered. Two examples will 

be given. 


First, consider diffraction by a ‘black’ screen. We shall show how the 


type of integral given by the usual Kirchhoff assumptions or their equiva- | i 6 
lent can be solved by the method of section 3. The cumbersome approach whi 
commonly given in text-books, whose validity is suggested rather than 
proved, is thus avoided. ve 
Ce 


Let the incident plane wave be specified by 


e—ikr cos(6—a) (92) 5 


where the ‘black’ screen occupies y = 0, 2 < 0 and r,6,z are cylindrical 
polar coordinates; polar coordinates p,¢ referred to an arbitrary point wh 
(€,0,0), so that x—é = pcos¢, y = psin¢d will also be used. 
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The assumption is made that some component of the field vanishes on 
the shadow side of the screen and maintains its undisturbed value in the 
0, so that the Kirchhoff integral is only evaluated over 
the aperture. The z integration yields a Hankel function, and the sort of 
result obtained is 

k | HP (kp)e-*é cos « dé, 


0 


(93) 


vhich mathematically puts Huygens’s principle in evidence. 
If the point of observation in y > 0 is a large number of wave-lengths 


from the aperture, we may write (93) as 


kf etki 


pinld e-iké cosa gg 
Num J Wp : 
0 
2k 3 
gn j—¢-Groowt—-a) p-texp| —thk{p—(x—£)cosa—ysina}| dé. (94) 
Al 7 ? 
0 
The change of variable from € to A, where 
. »(d—a , , — 
2p sin? m p—(x—&)cosa—ysina = —ird?, 
gives for (94) 
2 i dta ' 
oN (kre ikrekra cosee| 9 al dx, (95) 
. (d- 
where a city 2sin( 4 
2? 
Now (95) is of the form (33). At a point on the shadow edge, 0 = « and 


a=(. Also, at A ad G. 


approximation is given by (48); in the shadow region 6 


Hence, provided 6+-« is not small, the first 
x it is 


2 G+a . (O@—a« 
e714 accec p—ikr cos (6-a) F Ibvr)s 
om ( sec - P| y(2krsin( 5 |: 


x we must use (41), thus introducing the ‘geometrical optics’ term, 
which is the undisturbed plane wave. 


(96) 
If @ . 
For the second example, consider radio propagation over a finitely 


conducting, plane earth. For a two-dimensional source, the reflected wave 
an be put in the form 


») 
e~izl4 |" H’2)(hr) 
| 7 


f(sin x jetkr cosG +a) dx, 


| (97) 
Vie7) 


1x 


Where r, 


6 are measured from the transmitter’s image in the earth’s 
surface, and | 


/(sin «) is the Fresnel reflection coefficient for a plane wave 
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incident at an angle a. Ignoring contributions arising from the branch 
points of f(sin x), the second term in (97) can be written 


tor+in 
im/4 4 


| f[sin( x4 6) |e ikr cosa dx. 


é 


a (277) ‘ 


—\knr-—ix 


and the substitution A = e-'7/4y2 sin(4a) thus gives 
sec (4a) f[sin(a+6)]e*™ da. (98) 


If the waves are vertically polarized and @ is small, f[sin(«+4)] cannot 
be removed from under the integral sign in (98) at « = 0 (A = 0), becaug 
the fact that the earth’s (complex) dielectric constant is large means that 
f(sinax) has a pole near « = 0. However, the method of section 5 3 
applicable, leading to a first approximation in terms of the Fresnel integral 
with, in general, a complex argument (see 3). 
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